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IRE TRANSACTIONS ON ANTENNAS AND PROPAGATION 


Comparison of Experimental with Computed 
Tropospheric Refraction* 


W. L. ANDERSON}, MEMBER, IRE, N. J. BEYERS?, anp R. J. RAINEY}, sTUDENT MEMBER, IRE 


Summary—Limits of applicability of ray tracing in computing 
tropospheric refraction at White Sands Missile Range have been 
further explored. 286 comparisons were made, all for a path from 
radar to fixed beacon of about 45 miles and an elevation angle of 
17.99 milliradians. A horizontally stratified atmosphere was assumed. 
Refractive index profiles were prepared from a variety of weather 
data, and classified A, B, C, or R, in descending order of reliability 
prior to ray-tracing calculations. 

Angle observations were made with an FPS-16 C-band radar 
having a quoted instrumental accuracy of 0.14 milliradian- rms. 
Angle readings varied from 18.36 to 20.54 milliradians, with mean 
of 19.01 milliradians and standard deviation of 0.41 milliradians. 

The rms deviation of computed from experimental angles ranges 
from 0.28 to 0.41 milliradian for different classes of data. The ratio 
of this deviation to the deviation from over-all mean varies from 0.68 
for Class A to 1.00 for Class R. Thus the improvement over a “stand- 
ard atmosphere” varies from 32 per cent to 0, and correlates directly 
with quality of weather information. 

For this experiment it is concluded that most of the rms eleva- 
tion angle error is contributed by atmospheric conditions. Although 
ray tracing methods provide a significant correction when sufficiently 
good weather information is available, there still remains a large 
uncertainty not accounted for by equipment. 


* Received by the PGAP, December 4, 1959; revised manuscript 
received, March 12, 1960. This work was performed under Contract 
No. DA-29-040-ORD-1238. 


} Elec. Engrg. Dept., University of New Mexico, Albuquerque, 
aM 
t Missile Geophysics Division, White Sands Missile Range, N. M. 


INTRODUCTION 
IL) siseras the period July, 1958—July, 1959, a con- 


siderable amount of data was obtained at 

WSMR for the purpose of comparing computed 
with observed tropospheric refraction. Some results ob- 
tained prior to July, 1958 and presented in a previous 
report! are not presented here, because the recent data 
involves a different ray path and a more accurate radar. 
However, conclusions of the two reports are compared. 
On the basis of present information, it is possible to de- 
fine more definitely the magnitude of the elevation 


angle errors to be expected from meteorological uncer- 
tainties in the WSMR area. 


RADAR PATH AND TOPOGRAPHY 


The White Sands Missile Range is situated in the 
Tularosa Basin in southern New Mexico. The basin is 
oriented with its major axis north-south and is about 140 
miles long and 40 miles wide (Fig. 1). It is sharply 
bounded by mountain ranges which average 7500 to 
9000 feet MSL while the basin floor is about 4000 feet 


*'W.L. Anderson, N. J. Beyers, and B. M. Fannin, “Comparison 
of computed with observed atmospheric refraction,” IRE TRANs. 
ON ANTENNAS AND PROPAGATION, vol. AP-7, pp. 258-260; July, 1959, 
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1960 


MSL. The sandy soil on the basin floor supports desert 


vegetation such as mesquite, greasewood, cactus, and so 
forth. 


The radar unit used in this experiment is situated on 


3 the basin floor as shown in Fig. 1. All radar observa- 


_ tions were taken over a path extending from the radar 


site to a beacon target, equipped with a horn antenna, 
located in the Sacramento Mountains northeast of the 


radar. The slant range is 79,060 yards with the target 


es LS Sr 


Fy Le ee bk aK 


altitude 5610.5 feet above that of the transmitter, and 
the resulting elevation angle is 17.99 milliradians. Ex- 


cept for the mountains at the beacon end of the path, 


most of the terrain between transmitter and target is 


- quite flat. The topography is shown in cross section in 
_ Fig. 2, where the altitude scale has been multiplied by a 


factor of 7.5. 


CLIMATOLOGY 


The high plains region of southern New Mexico is 
characterized by short, mild winters and long hot sum- 
mers. Average temperatures range from 23°F low to 
61°F high in the winter and from 59°F low to 98°F high 
in the summer. The arid climate of the Tularosa Basin 


TULAROSA 


TULAROSA BASIN 


ALAMOGORDO 


AFB 


AFMTCC) 


HOLLOMAN 
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yields an average annual rainfall of about 8 inches, 
much of which is the result of the numerous thunder- 
storms occurring in late summer and early fall. These 
thunderstorms cause the formation and dissipation of 
convective cells which disrupt the horizontal stratifica- 
tion picture during this time of year. Clear skies gen- 
erally prevail during the winter months, while the spring 
is characterized by frequent winds averaging 20 to 30 
miles per hour. 


METEOROLOGICAL DATA 


The refractive index profiles used in computing ex- 
pected bending were prepared by personnel of the 
WSMR Missile Geophysics Division, U. S. Army Sig- 
nal Missile Support Agency. The data used in making up 
the profiles were obtained from radiosonde, wiresonde, 
surface observations, and helicopter mounted refractom- 
eter. At the time of preparing each profile, and before 
bending calculations were performed, a classification 
was assigned. The criteria for classification involved the 
following: ; 


Class A 


1) 
2) 


use of at least two observational methods, 

close agreement between observations made by 
different methods, 

obtaining all data within two hours of radar ob- 
servation, 

a stable synoptic situation, 

sounding made within 5 miles of propagation path. 


3) 
4) 
5) 
Class B 
1) 


one sounding available, 


2) 


3) 
4) 
5) 


1) 


RADAR 


CANTONMENT SITE MILES 


AREA 


1) 


WHITE SANDS MISSILE RANGE 
Migeads 


KILOFEET 


MILES 


Fig. 2—Terrain profile along ray path. 


obtaining data within two hours of radar observa- 
tion, 

a stable synoptic situation, 

soundings within 5 miles of propagation path, 
use of supplementary data more than two hours 
old if in good agreement. 


ClhassnG 


data which did not qualify for Class A or B, but 
were still considered to be a fair representation of 
existing conditions. 


Class R 


This category includes all profiles prepared solely 
from radiosonde, although in every case at least 
part of a radiosonde profile was also used in mak- 
ing up the Class A, B, or C profiles. The reason for 
this type of classification is to give a comparison 
between the Class A, B, or C results and those 
which would be obtained if the only available in- 
formation were radiosonde data, sometimes taken 
at times differing up to several hours from times of 
radar observations, and taken with no factor of 
meteorologist’s judgment entering the picture. 
Thus there are cases, for example, when the Class 
R profiles are the same as the B or C profiles. 
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TABLE I 
rms Per Cent 
Time of | Number of €rms 
Ce Day |Comparisons ee cae ay Improvement 
0.663 milliradian 0.426 milliradian 0.643 35.7 per cent 
A oak ‘6 0.288 milliradian 0.273 milliradian 0.948 5.2 per cent 
Both 37 0.534 milliradian 0.368 milliradian 0.689 31.1 per cent 
5 0.408 milliradian 0.307 milliradian 0.752 24.8 per cent 
B ea 2 0.292 milliradian 0.261 milliradian 0.894 10.6 per cent 
Both 111 0.354 milliradian 0.285 milliradian 0.805 19,5 per cent 
M. 8 0.482 milliradian 0.339 milliradian 0.703 29.7 per cent 
(Cc eat 9 0.319 milliradian 0.333 milliradian 1.044 —4.4 per cent 
Both 17 0.404 milliradian 0.336 milliradian 0.832 16.8 per cent 
A.M. 84 0.491 milliradian Wi trie me i: per ou 
i M. 81 0.294 milliradian 0.272 milliradian : .5 per cen 
cae Both 165 0.406 milliradian 0.311 milliradian 0.766 23.4 per cent 
A.M. 75 0.478 milliradian 0.470 milliradian 0.983 1.7 per cent 
R P.M. 68 0.272 milliradian 0.330 milliradian An 213 —21.3 per cent 
Both 143 0.394 milliradian 0.409 milliradian 1.038 —3.8 per cent 


Ray-Tracing Method 


The procedure used in ray-tracing has been fully 
described in previous reports,” and will not be further 
discussed except for the remark that the computational 
accuracy is conservatively estimated to be within about 
0.03 milliradian. 


RADAR READINGS 


The internal precision of the FPS-16 radar is quoted 
at 0.10 milliradian rms, but since additional fixed error 
or “bias” can occur, the over-all internal accuracy of the 
radar is quoted at 0.14 milliradian rms. This is consid- 
ered a conservative figure. 

A total of 165 observations were made during the pe- 
riod July, 1958 to July, 1959. Of these, 84 were taken 
at around 0900 and 81 around 1400 local time. Mean ob- 
served angle was 19.01 milliradians, representing an 
average bending of 1.02 milliradian. Extremes in ob- 
served angles were 18.36 and 20.54 milliradians. Stand- 
ard deviation was 0.406 milliradian. 

In nearly all cases the so-called “observed angle” is 
the mean of 5 readings taken at 1-minute intervals. The 
maximum difference between extremes during ‘any 
single 4-minute period was 0.58 milliradian, but this 
was quite unusual. There were no other cases in which 
the difference exceeded 0.25 milliradian. The mean dif- 
ference between extremes was 0.07 milliradian, and the 
rms was 0.16 milliradian. Since bias is probably con- 
stant during a 4-minute interval, the internal radar 
precision figure, 0.10 milliradian, can be used to arrive 
at a value for the rms contribution of random atmos- 
pheric fluctuations, or scintillations, 


+/0.16? — 0.102 = 0.125 milliradian rms. 


* Ibid. ; and B. M. Fannin and K. H. Jehn. “A study of radar 
elevation-angle errors due to atmospheric refraction,” IRE TRANS. 
ae eee AND PropaGaTion, vol. AP-5, pp. 71-77; January, 


This should be considered a very tentative figure only, 
since at present there is no way of knowing how much of 
the change during an observation period is contributed 
by random fluctuations about a mean, and how much 
by sudden shifts in the refractive index profile. 


Comparison, Observed with Computed Values 


An evaluation was made of 165 Class A, B, and C 
profiles. In addition, 143 Class R comparison profiles 
were considered. Table I gives a breakdown of the re- 
sults, subdivided in various ways. The quantity € is de- 
fined as (X — Y), where X is the computed angle and Y 
the observed angle. The quantity ay is introduced and 
defined as 


N 
>> (Y; — 19.01)? 
i=1 


ay = ) 


N 


where 19.01 is the mean observed angle for all 165 cases 
studied. ay is used primarily as a measure of departure 
of any group of recorded angles from the over-all mean, 
and in the special case where all 165 readings are con- 
sidered, ay is the standard deviation. 

The ratio €:ms/ay has been adopted as a criterion for 
judging results for two reasons. First, it tells at a glance 
whether any improvement was gained over a “standard 
atmosphere” correction, and second, it makes different 
samples of data directly comparable, whether or not the 
samples differed among themselves in their statistical 
properties. Thus, although in footnote 1 the value of 
€rms for Class A data was 0.17, and in this report it is 
0.37, the ratio to ay in each case shows that approxi- 
mately the same degree of improvement was effected, 
namely, a reduction of rms error to about 30 per cent 
to 40 per cent below ay. Considering both reports, the 
reduction for Class B lies in the neighborhood of 20 per 
cent to 40 per cent, and for Class C, —40 per cent to 20 


) 


1960 


per cent. Since the number of Class C cases is quite 
small in both reports, probably the safest conclusion is 
that Class C gave no improvement. Nor was there im- 
provement from Class R data. 

It will be noted from Table I that the rms values for 
€ms are considerably higher in the morning than in the 
afternoon. This fits the pattern of the diurnal meteoro- 
logical cycle, as discussed by various investigators. 
Day and Trolese conducted experiments in the Arizona 
desert during winter months, and obtained averaged 
profiles of refractive index hourly over 24 hours. At 
1400, the atmosphere was nearly standard with time 
rate of change near minimum. On the other hand, when 
observed at 0900, the refractive index near the surface 
was high, owing to existence of a temperature inversion 
formed during the night, but it was beginning to de- 
crease quite rapidly due to heating of the ground by the 
sun. In the desert type terrain of the Tularosa Basin, in 
the WSMR area, the same general conditions prevail. 
During periods of time instability, there is almost cer- 
tainly a concomitant space instability, with the result 
that the horizontal stratification assumption is less valid 
than at other times. This may seem to contradict the 
requirement of “stable synoptic situation” imposed in 
classifying the A or B profiles; however, the judgment 
of stability is based on an over-all profile which ‘is of 
necessity local to the weather sounding equipment. An 
exception to this limitation occurs when data are ob- 
tained by helicopter mounted refractometer, but here 
the data collecting time may introduce some uncer- 
tainties. 

As an indication of the seasonal behavior of the di- 
urnal cycle, it is instructive to compare a@y’s for various 
periods of the year. The trend is illustrated in Table II, 
where the 3-month groups have been chosen in such a 
way as to isolate the conditions in the period centered, 
more or less, on the winter solstice. During this time, 
surface heating in the hours before 0900 should be con- 
siderably less than for other times of the year. 

These figures seem to fit the general weather picture 
very well. In all cases ay(A.M.) is higher than ay(P.M.), 
and the tendency is most striking for the November— 
January period. During August—October they are nearly 
the same, and both are quite high, probably because of 
occasional rains during this time. In the February—April 
period, ay for both morning and afternoon is fairly low, 
due perhaps to winds and consequent good mixing of 
the atmosphere. The minimum is reached for the May— 
July period, due to ground heating well before 0900 and 
the general dryness. A contributing factor may also be 
May winds, which are not unknown in this area. 

The tendency for €ms to be high when the data have 
high ay is therefore probably related more to the insta- 


3 J. P. Day and L. G. Trolese, “Propagation of short radio waves 
over desert terrain,” Proc. IRE, vol. 38, pp. 165-175; February, 
1950. sit ’ 
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TABLE II 
: : No. of Cases 
Time Period (A, B, C) ay 
August through October A.M. 11 0.54 
P.M. 9 0.53 
Combined 20 0.53 
November through January A.M. 24 Ond2 
P.M. 25 0.25 
Combined 49 0.53 
February through April AM 32 0.39 
P.M. 27 0.33 
Combined 59 0.36 
May through July A.M. 17 0:21 
P.M. - 20 Oon3 
Combined Ryi 0.17 


bility of the meteorological situation near the ground 
than it is to a failure of the ray-tracing procedure when 
applied to extreme cases. Unfortunately, the data pres- 
ently on hand do not include enough cases of high bend- 
ing during times of supposed high stability to give any 
valid statistical weight either for or against the hy- 
pothesis. Nevertheless, the low figures for ém; at 1400 
lead to a strong presumption that, during the other 
stable period, around 0200, it should also be low, even 
though bending at this time would be expected to be 
high. 

Another point that deserves special comment is the 
very high ay for Class A (A.M.) cases, as contrasted with 
B and C. This arises because the classification A gen- 
erally necessitates the inclusion of wiresonde data. The 
use of wiresonde is restricted to days on which winds are 
less than 10 knots; consequently, the atmosphere is not 
as likely to be well mixed at such times as it is otherwise, 
and a variety of profiles may occur, which have sharply 
stratified moisture layers, temperature inversions, etc. 
Class A (pP.M.) cases, on the other hand, have the lowest 
ay of any of the groups because on clear, calm after- 
noons the morning inversion has generally disappeared, 
and the temperature distribution is near adiabatic. 

In the scatter diagrams (Figs. 3 and 4), the main 
diagonal represents perfect agreement between ob- 
served and computed values, while the two unbroken 
lines on either side show the radar accuracy. The dotted 
lines show the standard deviation of all observed angles, 
and are intended merely as a rough guide in evaluating 
the over-all results. 


Surface Refractive Index 


In the course of this study, the following question has 
frequently arisen: might there not be an adequate angle 
correction procedure based only on a knowledge of sur- 
face refractive index? Such a procedure, as compared 
with ray-tracing from a complete profile, would greatly 
mitigate the amount of effort required in making correc- 
tions. 
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One such procedure that has been proposed! involves 
the assumption that a fair representation of weather 
conditions is given by a refractive index profile which is 
exponential between Np (ground value of V) and a fixed 
value V;=128 N-units at 20,000 feet above the surface. 
Specifically, the profile is described by: N= Noe-8’, 
where Y is height above surface and N the correspond- 


4K. E. Pearson, D. D. Kasparek, and L. N. Tarrant, “The Re- 
fraction Correction Developed for the AN/FPS-16 Radar at White 
Sands Missile Range,” U.S. Army Signal Missile Support Agency, 
White Sands Missile Range, N. M., Tech. Memo. 577; November, 
1958. 
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ing N-value. For given No, B is obtained by evaluating 
the expression at 20,000 feet. Above this height, an aver- 
age profile is used. 

Profiles made up in this way can be ray-traced to give 
a table of correction factors for the pertinent range of 
N>. These factors are used in a simple formula which 
gives the angle correction for a given range and angle. 

As is shown in Table III, the method gives correc- 
tions which, statistically, are a slight but probably not 
worthwhile improvement over those that would result 
from a fixed standard atmosphere. The average im- 
provement is 7.5 per cent, or about 0.03 milliradian. 
Since the method does not involve a classification of the 
weather situation, the results are given when this 
method is applied to all the cases which were used in 
Classes A, B, and C. It is perhaps worth mentioning, 
however, that €ms/ay varies only slightly for different 
samples of data, provided the samples contain a reason- 
able number of readings. 


TABLE III 
CORRECTIONS BASED ON SURFACE REFRACTIVE INDEX ONLY 


Time Num- 

Group ber ae bs érms/ay 
A.M, 84 0.450 milliradian 0.491 milliradian 0.916 
P.M. 81 0.280 milliradian 0.294 milliradian 0.952 
Combined | 165 0.376 milliradian | 0.406 milliradian 0.926 


CONCLUSIONS 


The following main conclusions are drawn. 

1) The high accuracy of the radar used in this ex- 
periment allows a fairly definite conclusion as to the 
main source of the discrepancies between computed and 
observed elevation angles. The rms error for combined 
A, B, C data amounts to about 0.31 milliradian, and 
taking radar accuracy to be 0.14 milliradian, the error 
contribution from weather is evidently about 0.28 milli- 
radian. This is of nearly the same magnitude as the 
total error. 

If this weather contribution is assumed to include an 
rms error of about 0.13 milliradian arising from scintil- 
lations or other sudden, unpredictable changes, and if 
this is subtracted quadratically, then the remaining 
meteorological uncertainties are given by the figure 0.25 


milliradian rms. This quantity should represent two 
main factors: 


a) limitations in determining the true refractive index 
profile, 


b) departures from the horizontal stratification pic- 
ture. 


2) In evaluating the extent to which bending can be 
predicted, it is apparent that when weather deviates 
considerably from the average condition, the probabil- 
ity of error in computing the bending is increased. This 
probability is minimized, however, when best quality 
meteorological information is used. The correlation be- 


a ie 
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1960 


tween quality of weather data and the improvement 
gained over a standard atmosphere is direct and definite, 
both regarding Anderson, et al.,! and this paper. 

3) The Class R data, consisting entirely of radiosonde 
profiles, give no reduction of error over that which would 
have resulted from using a standard atmosphere. This 
and the Class C results indicate quite strongly that radio- 
sonde information alone, under the conditions of this 


_ test, is virtually worthless for ray-tracing purposes. 


4) The corrections deduced from surface refractive 
index, by the method described in the previous section, 


_ gave an improvement of about 0.03 milliradian rms 
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over the error which would have resulted from a 
standard atmosphere assumption. This is considered 
insignificant for most practical purposes. 

5) In the WSMR area, particularly, the time of day 
and season of year constitute valuable information on 
a) what type of average profiles can be expected, and b) 
what type and degree of variation from these averages 
can be expected. Possibly this knowledge can be guided 
by surface refractive index and initial gradient readings 
taken near the radar site to obtain corrections rivaling 
those obtained by ray-tracing from complete Class A 
profiles. 


The Flared Slot: A Moderately Directive 
Flush-Mounted Broad-Band Antenna* 


J. W. EBERLE}, C. A. LEVIS}, seNIok MEMBER, IRE, AND D. McCOYt, MEMBER, IRE 


Summary—tThe flared slot is a simply constructed antenna de- 
signed to be flush-mounted in a horizontal surface. Its pattern is a 
moderately directive pencil beam over a 4:1 frequency band, and 
its impedance bandwidth is essentially that of a coaxial line-to- 
ridged waveguide transition. Both vertically and horizontally polar- 
ized models are discussed. 


I. INTRODUCTION 
Oat ac having a single, smooth beam of mod- 


erate directivity are finding increasing use in air- 

craft and airport installations in connection with 
homing or ground-directed approach systems. The Yagi 
arrays and horns commonly used to obtain such patterns 
are frequently inapplicable to high-performance air- 
craft or to airport installations where the only available 
locations often are flat surfaces, such as the underside of 
the fuselage, wings, or horizontal stabilizer, or the 
ground near or even in the runways. A class of traveling- 
wave antennas suitable for installation in such locations 
was briefly described by Stephenson and Walter [1]. 
The further development of one of these antennas, 
named a flared slot antenna, is discussed below. Its con- 
struction, far field patterns and impedance over a wide 
range of frequencies, bandwidth, and near-field struc- 
ture are discussed. The effect of certain modifications 
which might be necessary in typical installations has 
also been investigated. 


* Received by the PGAP, August 3, 1959; revised manuscript 
received, March 3, 1960. This work was performed at The Ohio State 
University, Columbus, and was sponsored by the Air Research and 
Development Command, Wright Air Development Center, Wright- 
Patterson Air Force Base, Ohio. : 4 

+ Antenna Lab., Dept. of Elec. Engrg., The Ohio State Univer- 
sity, Columbus, Ohio. 

t General Electric Co., Syracuse, N. Y. Formerly at Antenna 
Lab., Dept. of Elec. Engrg., The Ohio State University, Columbus, 
Ohio. : 


II. VERTICALLY POLARIZED FLARED-SLOT ANTENNA 


A vertically polarized flared-slot antenna is pictured 
in Fig. 1, and the basic construction appears in Fig. 2. 


Fig. 1—Flared-slot antenna for vertical polarization. 
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Fig. 2—Basic construction of vertically polarized antenna. 
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Energy entering at the coaxial feed is transferred 
through a ridged dielectric-filled waveguide section to 
the radiating aperture. The coax-to-waveguide transi- 
tion was designed with the use of ridged-guide data by 
Cohn [2], [3] to give a wide-band transition to 50-ohm 
cable. The antenna itself is based on a traveling-wave 
design by C. H. Walter. It employs a phase velocity 
somewhat slower than that of free space and exhibits a 
clean endfire pattern over a 4:1 frequency band. 

The impedance locus of the antenna is shown in Fig. 3 
as measured on a Hewlett-Packard slotted coaxial sec- 
tion through UG 491 A/U and UG 201/U adaptors. The 
impedance bandwidth is limited by cutoff at the low-fre- 
quency end, and by the waveguide back plate at the 
high-frequency end, where it appears as an effective 
short-circuit at the feedpoint. The bandwidth for which 
the VSWR does not exceed a given value is plotted in 
Fig. 4. 


240, 
I 


180}— 


fo-f 
% Bandwidth = + x 100 
j ‘ 
Where fz Is The Highest And 
The Lowest Frequency For Which 
The Given VSWR Is Not Exceeded. 


Z, = 50 Ohms =} 


| 
} l 


7 8 9 10 


5 
VSWR 


Fig. 4—Impedance bandwidth of vertically polarized antenna, 
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The far-field strength (voltage) patterns of the an- 
tenna are shown in Fig. 5, which pertains to the resin- 
filled model of Fig. 1, and in Fig. 6, which gives the pat- 
terns of an earlier paraffin-wax-filled model. The ~ 
coordinates used appear in Fig. 7 (p. 464). The resin 
used was Stycast 35 and had a relative dielectric con- — 
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Fig. 5—Field-strength (voltage) patterns of the vertically polarized 
flared-slot antenna of Fig. 1. 
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Fig. 6—Field-strength (voltage) patterns of a vertically polarized flared-slot prototype. Horizontal plane patterns at 6 
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p= 0° 
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Fig. 7—Coordinate system for field-strength patterns. 


stant of 2.58. The dielectric constant of paraffin is 2.25. 
The extremely wide range of endfire operation (over a 
4:1 frequency band) is clearly evident in Figs. 5 and 6. 
The limitation at the high-frequency end is due to 
higher-order modes in the waveguide which are launched 
into the aperture and cause beam deterioration. At the 
low-frequency end, the discontinuity at the beginning 
of the aperture can cause spurious radiation. A “discon- 
tinuity minimizer” consisting of transverse Aquadag 
strips imbedded in the dielectric (Fig. 2) is an effective 
means of reducing this effect with some sacrifice of 
efficiency. Fig. 8 (p. 465), shows the effect of two differ- 
ent sizes of the discontinuity minimizer. The efficiency 
of this particular antenna has not been tested, but 
measurements on a paraffin-filled antenna without ridge 
feed showed an efficiency reduction from 80 per cent to 
to 55 per cent when the minimizer was used. For receiv- 
ing applications this is not serious, since the directivity 
increases, so that gain is not affected adversely. 

The structure of the antenna near field at midband 
(Xo =5.017 cm) is shown in the phase and amplitude 
plots of Figs. 9 and 10, respectively. The quantity 
plotted is the electric field component tangential to the 
dielectric-air interface and parallel to its centerline at a 
distance of 0.005 cm from the interface. The orthogonal 
tangential component was lower by 25.1 db over the 
entire surface. The measurements were made with a 
phase-sensitive coherent detection system [4], [5] de- 
veloped at this Laboratory. Contours are not plotted 
near the discontinuity minimizer because the probe was 
too large to resolve the rapid variations of the field in 
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8=0° 


this region. No data were taken at the far edges of the 
antenna. 

The polarization of the antenna far field is pure verti- 
cal; over the entire frequency band, the cross-polarized 
component is 15 db below the maximum of the desired 
component. This can easily be understood in terms of 
the measured near field. Referring to Fig. 11 (p. 466), 
the longitudinal electric field in the aperture is equiva- 
lent to a set of magnetic currents flowing in the surface 
at right angles to the indicated field lines. The equiva- 
lent magnetic currents are related to the electric far 
field as ordinary electric currents are to the magnetic 
field, and thus they produce the far field shown, which 
is vertically polarized. 


III. HortzoNTALLyY POLARIZED 
FLARED-SLOT ANTENNA 


A horizontally polarized flared-slot antenna is shown 
in Fig. 12. The basic construction appears in Fig. 13, 
and a method of feeding is sketched in Fig. 14. The sim- 
ilarity of this antenna to the vertically polarized model 
is immediately apparent. The major differences are in 
the feed probe, which is parallel to the aperture for hori- 
zontal polarization, and the greater feed section and 
aperture depths required to propagate the resulting 
waveguide mode. 

The direction of the electric field in the aperture is 
sketched in Fig. 15. The corresponding equivalent cur- 
rents are perpendicular to the electric field and are also 
indicated in the figure, and the resulting electric far field 


, | 
September — 


4 


1960 Eberle, Levis, and McCoy: The Flared Slot 465 


e' 
va 
‘ 
i 
a 


A SSRN 
Lz bs: Sy Nyy, 
ps SOY 
ED ] 
==) 
ZOMG 
No Lossy Strips ” {2 
(a) 


(d) 


Fig. 8—Discontinuity minimizer effect on ‘pattern of vertically 
polarized antenna. Dimensions are given in inches. The shaded 
strips represent Aquadag paint. 
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Fig. 11—Polarization relations of vertically polarized flared-slot 
antenna. Solid arrows denote electric field lines; dashed arrows, 
equivalent magnetic currents. 
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Fig. 12—Horizontally polarized flared-slot antenna. 


Fig. 13—Basic construction of horizontally polarized antenna. 
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Fig. 14—Feed section of horizontally polarized flared-slot antenna. Fig. 15—Polarization relations of horizontally polarized antenna. 


Solid arrows represent electric field lines; dashed arrows re- 
present equivalent magnetic currents. 
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Fig. 16—(Continued.) 


is shown. It is apparent that a substantial portion of the 
main beam is horizontally polarized, but there are two 
vertically polarized regions at the beam edges. These ap- 
pear as distinct lobes in most of the horizontal-plane 


WV, SH 
eS LOSE e 7 patterns. 
LR CR The far-field patterns of the antenna of Fig. 12 are 
Ee Ah Sl M shown in Fig. 16. The dielectric is paraffin wax, and no 


discontinuity minimizer was used. If a discontinuity 
minimizer were to be employed to extend the antenna 
bandwidth, its lossy strips should run longitudinally, 
since this is the direction of the unwanted electric field 
component. The effect of such a discontinuity minimizer 
is shown in Fig. 17. 

The horizontally polarized model received less de- 
velopment effort than the vertically polarized version. 
All models employed paraffin wax as dielectric, and no 
attempt was made to broadband the impedance char- 
acteristics. 
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Mathematical and Experimental Studies of a 
Wide-Band Vertically Polarized Antenna* 


P. FOLDES{, MEMBER, IRE 


Summary—Theory and experimental results are presented for 
an omnidirectional, vertically polarized antenna. One of the de- 
veloped models has 7 db gain, VSWR <1.7, and 50 kw power han- 
dling capacity in a 2:1 frequency band. The new antenna is simple in 
construction and has a single coaxial feed point. 


INTRODUCTION 


N recent years considerable effort was exerted in de- 
| veloping very wide frequency banded antenna sys- 

tems. The definition of “wide-band” was continu- 
ously extended and at present this term usually refers to 
a band which has an end frequency ratio of about 1.5:1 
or more. 

In the present paper, an antenna type which belongs 
to this group will be discussed. It will be shown that a 
circularly symmetrical, vertically polarized antenna sys- 
tem can be obtained by means of a coaxial feed horn and 
a circularly symmetrical reflector. The antenna has 
omnidirectional pattern in the horizontal plane and 
highly directive or shaped beam in the vertical plane. 
The feed system has a coaxial input, with a very low 
VSWR. The antenna system does not contain any res- 
onant element and therefore its impedance characteris- 
tics remain practically constant in a very wide frequency 
band. At the same time, the reflector configuration as- 
sures that the pattern changes only slightly in a 2:1 fre- 
quency band. The novel antenna structure can be re- 
lated to a discone-type antenna.' However, its feed sys- 
tem can be considered as an improvement because of the 
mechanically favorable arrangement. The system is 
particularly useful for very-high-power application in the 
UHF band. Theoretical explanations will be given for 
the feed system: and experimental verification for the 
performance of the new antenna. 


GENERAL DESCRIPTION OF OPERATION 


The basic electrical arrangement of the proposed an- 
tenna is shown in Fig. 1. This represents the principal 
parts of the antenna, namely: 


1) a rigid coaxial transmission line, propagating the 
TEM mode; 

2)-a toroidal transition section, which inverts the 
inner and outer conductors of the coaxial line 1); 

3) a coaxial horn section, which feeds the circular 
aperture with a TEM mode; 

4) aniris, to cancel the aperture reflection of the horn 
‘in a wide band; and 


* Received by the PGAP, December 1, 1959; revised manuscript 
received, March 18, 1960. 

+ RCA Victor Co., Ltd., Montreal, Quebec, Can. 

1J. J. Nail, “Designing discone antennas,” Electronics, vol. 26, 
pp. 167-169; August, 1953. 
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Fig. 1—Electrical schematic of UHF antenna 
for vertical polarization. 


5) a reflector, having particular curvature for radia- 
tion of energy into space. 


The operation of this antenna can be summarized as . 
follows: the antenna consists of two main parts, the 
feed system and the reflector. The feed system is basi- 
cally a transition between the small cross section of the 
coaxial line and the considerably larger cross section of 
the coaxial horn aperture. The reflector transforms the 
primary spherical wave, radiated by the horn, into the 
far field distribution. ; 

Fig. 1 indicates some of the details which are neces- 
sary in order to obtain the desired operation. A tapered 
coaxial transmission line section serves at the input of 
the antenna. This line section transforms the charac- 
teristic impedance of the 50-ohm feed line to a smaller 
value. This lower impedance is necessary to minimize 
the dimensions of the toroidal line section. For low 
impedance lines, the diameter of the outer conductor 
does not differ much from that of the inner. Therefore, 
the decrease of impedance must be accompanied by an 
over-all increase in the diameters to minimize losses. 

The toroidal section 2) can be designed for constant 
characteristic impedance. This is achieved by varying 
the distance between the two conductors along the 
direction of propagation. At the same time, the radius 
of curvature of this section must be kept large compared 
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to the difference in radii of the conductors. Optimum 
dimensions will result if the ratio of the distance be- 
tween the conductors to the radius of curvature re- 
mains constant. This constant ratio will result in a 
“quasi-toroidal” transmission line having a spiral cross 
section through the axis of the conductors. The toroidal 
line is terminated by a coaxial horn 3). This section must 
provide a matched termination in the required fre- 
quency band and, at the same time, it must illuminate 
the reflector in a specific way. The impedance character- 
istic of this horn is principally determined by its length 
and the angle between the inner and outer conductors. 
Moreover, the radiation direction of the horn depends 
on this angle and the size of the aperture. Therefore, the 
radiation pattern of the horn can be controlled by these 
parameters in a relatively wide range. It should be re- 
alized that the cone angle, which is determined by a line 
bisecting the vertex angle of the horn, not only shapes 
the curvature of the phase front emerging from the 
horn, but also has some effect on the possibility of 
higher modes being propagated. Such higher modes 
might result in cross polarized components and asym- 
metries in the far field. Therefore, the possibility of these 
modes being propagated should be eliminated by main- 
taining the cone angle at a small value. 

The near field of the coaxial horn is reflected by the 
reflector 5). This reflector can have quite general char- 
acteristics. Namely, it can be madeas a “reflecting lens,” 
which results in specified amplitude and phase char- 
acteristics in the outgoing wave. In the following, de- 
tailed experimental data will be presented only for the 
case when the reflector is a wire mesh surface. 


GEOMETRY OF THE FEED SYSTEM 


From the preceding paragraph, it follows that a point 
source is required for illuminating the reflector. This 
point source can be a coaxial horn, where a wide-fre- 
quency-band operation is assured. This horn should have 
a certain aperture size for the required pattern function, 
while the input of the antenna should be a standard size 
coaxial transmission line. Thus, a transformer section is 
required between this aperture and the input line. For 
mechanical reasons, it is convenient to feed the antenna 
at the bottom, so this transformer section should also 
invert the direction of propagation. 

The characteristic impedance at the input of the an- 
tenna is 


bo 
Zo = 60 ln — 


a 


ohm. 


(1) 


In the first section of this transition, this impedance is 
transformed to a smaller value, 


bi 
Z;1 = 60 In — 
ay 


ohm, (2) 


where bo, b; and do, a; are the diameters of the outer and 
inner conductors at the beginning and at the end of this 
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section. The attenuation per unit length within this sec- 
tion is 


ae 
F a b ee G3) 
eer b i b aes 
In— 
a 


where 7 is a constant which depends on frequency and 
material, and a and b are the diameters of the conductors 
at a particular cross section and tr =b/a. The design cri- 
terion for this section is a constant loss in any cross sec- 
tion. Thus, 

1 7+1 a 


= — = const. 
6b Int T 


(4) 


The length of this section can be about 1.5 of the maxi- 
mum wavelengths, if a conical taper is used, and some- 
what less if a Tchebycheff-type stepped transformer is 
used.? 

A toroidal transmission line can follow the input 
transformer. The geometry of this line section can be 
analyzed by a simplified theory, assuming that two loss- 
less conducting strips are arranged parallel to each 
other (see Fig. 2). Let their width be A, their separa- 
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Fig. 2—Generation of the toroidal line. (a) Strip line. 
(b) Cylindrical line. (c) Toroidal line. 


tion B, and their length extend to infinity. If BA, the 
effect of the edges is negligible and a pure TEM mode 
can be propagated between the strips. The phase con- 
stant will be equal to that for free space, 8 = 27/d. Under 
these conditions, the transverse field components £, and 
H, will exist as indicated in Fig. 2(a). This transmission 
line has an inductivity? L=y(B/A) and a capacity 
C=e(A/B) per unit length. Now let us bend this strip 
line around an axis which is parallel to the width A or 
the y-direction [(see Fig. 2(b)]. If the two strips run 
along a circular arc, the field components for these por- 
tions may be expressed in cylindrical coordinates. For 
such a system of coordinates, the E, component becomes 
an , component, with the H, designation remaining 
unchanged. However, boundary conditions necessitate 
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the existence of an Ey component. For p> A, this com- 
ponent will be very small and hence the new mode will 
_ be very similar to the previous TEM mode. For the 
_ cylindrical part of the arrangement the inductance L’ 
_ and capacitance C” per radian are, respectively, 


_ B(p2” — p21’) _ BAp 


Ih i and (5) 
er eA 
= 6 
1 p2 ( ) 
n — 
Pi 


~_ where A=p2—pi1 and p=(pitpe2)/2, and p; and p2 are the 
Z radii of curvature of the strips. Let us assume now that 
- the line will be bent around the original g-axis too. In 
_ this way, a toroidal system is obtained and Fig. 2(c) 
_ shows a small section of it. The last shaping of the lines 
- does not introduce any new field component, but the 
width A of the strips must now be measured along the 
_ new arc, which is now a full circle. If the radius of this 
second forming operation is large relative to the separa- 
- tion A between the two conducting surfaces, the char- 
acteristic impedance can be calculated from the L’ and 
C’ values defined above. At an angle of ¢ (see Fig. 2), the 
_ value of A can be computed as 


A = 2nl[rx + p(1 — cos ¢)], (7a, 


' where 


and A, is the conductor distance at the input of the to- 


- roidal line. 


Thus the value of the characteristic impedance and 
the phase constant for the toroidal section can be writ- 
ten 


ae 1/2 60A p pe 1/2 
c rx + p(1 — cos¢)LA pi 


60A | Ne | 
~ 1+ 
rz + p(1 — cos ¢) 24p? 


60A 
= (8) 
rx + p(1 — cos 4) 
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A? 
= 5o(1 a ) = Be. (9) 
p 


The characteristic impedance at the input of this sec- 

tion is 

60Ao mg 120Ao 
1K aN 


Ri = (10) 
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The design objective for this toroidal section is to 
maintain a constant characteristic impedance K = Ko 
throughout the length of this section. This means 


120Ao 120A 
re a ae (11) 
b — Ao b — Ao + 2p(1 — cos ¢) 
from which 
eA E beevelak re, | 12 
= — cos : 
0 rea ¢ (12) 
From (4), it is seen that 
r+1T 
ae Ae (13) 
bot aed 
and 
a 
eee 
phe Est b re—1T (14) 
en) 2 pint 2a 
The equation for A becomes 
r— 1 2p(1 — cos 4) 
A= 1+ BeeGlss 
a ool Le OTF | se 
l Inr a tlnr 2a | 


Eq. (15) gives the values of the separation A in terms 
of the angle ¢. In this equation there are two parame- 
ters, p and r. The radius must increase in proportion to 
the separation A, which in turn must be kept small if 
small over-all size is desired. 

If p=5A, then (15) becomes 


Eo 


oo Rig anes Soe 
in which 
Tr —1 10(is— cos 
Ai oes ee ca eee se 
ee tl 

and in this case 

p= 9 eel (18) 
1 — EoF 


Fig. 3 exhibits the typical dimensions of a feed system, 
calculated in accordance with the above equations for 
the 225—400-mc frequency band. 

The size of the mouth of the horn was determined 
from the aperture distribution required for the antenna 
system. (The remaining dimensions were calculated 
from the input VSWR and power handling capacity re- 
quirements.) 

An inspection of Fig. 1 suggests that the antenna sys- 
tem contains several reflection sources from the input 
VSWR standpoint. The possible sources are: 
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Fig. 3—Layout of small-scale experimental antenna. 


1) the central part of the reflector, 

2) the aperture of the horn, 

3) the end of the toroidal line, 

4) the input of the toroidal line, and 

5) the input of the impedance transformer. 


The larger part of these reflections can be calculated in a 
straight forward way and will not be discussed here. 
However, the results of these calculations show that the 
individual reflections are less than a few per cent and a 
very low input VSWR can be obtained in a 2:1 fre- 
quency band. Some experimental data will be presented 
in the following paragraph. 


EXPERIMENTAL RESULTS 


Though the possible use of the above discussed feed 
system and reflector structures is quite wide, experi- 
mental data will only be given for one particular applica- 
tion. In this case a medium gain (6~8 db) antenna was 
required, with slight beam shaping in the vertical plane. 

Two small-scale models were constructed for the 2.0— 
3.8 kmc frequency band. Fig. 4(a) and 4(b) show the 
two models. The model in Fig. 4(a) is a smaller gain an- 
tenna and was built mainly for pattern measurements. 
Therefore, it does not contain the toroidal line section. 
However, the input VSWR data of this model are char- 
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(b) 
Fig. 4. 


acteristic for the reflection of the horn aperture and re- 
flector and were less than 1.32 in the whole frequency 
band. At the same time, the gain of this model was in 
the order of 6 db (see curve A on Fig. 8), and the cross 
polarized (horizontal) field component was at least 14 
db less than the vertically polarized component. 

The model represented in Fig. 4(b) is closer to the 
final form of the antenna, therefore more details will be 
given about the characteristics of this antenna. 

It was experienced that the inverting transition in the 
feed horn was a rather critical component of the system. 
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— 1960 
_ This transition affects the impedance behavior of the 
_ antenna to a significant degree and, moreover, has a 
- great influence on the appearance of higher-order modes 
- and any such higher modes will disturb the gain char- 
- acteristic of the antenna. During these experiments it 
_ was found, moreover, that a rather simple transition 
_ yielded a VSWR of 1.7 or better, together with the de- 
_ sign objective for the gain of 7 db. 

If the coaxial horn has a throat impedance equal to 


+> 


wae 


. ] 


a the characteristic impedance of the input line, a match 
_ will result. This condition is more easily obtained to- 
_ wards the lower-frequency end of the band. As the fre- 
quency is increased, the energy stored in the joint be- 
comes larger and hence any junction effect will become 
more pronounced. This junction effect can be decreased 
be keeping the characteristic impedance at the junction 
_ between input line and horn at a low value. 
i Based on this consideration, an impedance of 25 ohms 
was chosen for the input of the coaxial horn. Calcula- 
~ tions have shown that a four-step Tchebycheff trans- 
_ former is capable of achieving an almost perfect match 
_ between the 25 ohms and the 50 ohms impedance of the 
_ feed line. 
The four steps are necessitated by the wide frequency 
__ band to be covered. During the experiments, however, it 
; was found that a single quarter-wave transformer al- 
_ ready gives an impedance match within the specified 
~ limits (VSWR <1.8). The actual transformer was built 
by inserting a Teflon spacer 0.9 inch long at the cross 
section of the junction, giving a section of 35 ohms, one 
- quarter-wave long in the middle of the frequency band. 
Physically, the reflector was constructed from 64 rods 
formed to the proper (nearly parabolic) shape. Only half 
_of these rods extended the full length of the cross sec- 
tion; the other half was used to fill in the spacing be- 
tween the full-length ones at the greater radii. Such a 
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reflector properly built can reasonably be assumed to 
support a surface wave besides reflecting-an incident 
wave. Thus, such a reflector can be the source of a com- 
plex radiation. The upper portion will act more as a re- 
flector, radiating a broadside pattern, with the lower 
part supporting a traveling wave, 7.e., radiating as an 
endfire antenna. 

The existence of a separate guided and reflected wave 
can be proved relatively easily by varying the distance 
between the reflector and the horn. If the separation is 
small, then practically no reflected wave can occur and 
the energy is propagated along the gradually increasing . 
center conductor of the coaxial horn. The direction of 
the radiation maximum encloses a very small angle with 
a line which is tangent to the end of the reflector. This 
radiation can be considered as an endfire radiation gen- 
erated by the guided wave along the reflector. The angle 
of maximum radiation decreases with the size of the re- 
flector.* s 

If the separation between the horn and dish is large, 
the latter can be considered as a pure reflector, and the 
radiation maximum can be determined approximately 
by means of geometrical optics. In this case the angle 
of maximum radiation is independent from the size of 
the reflector and is determined by the relative position 
of the focal point of the reflector and the radiation 
center of the horn. In the intermediate range, both radi- 
ations exist and depending on the energies content and 
phase relationship of the two modes, various pattern 
shapes can be generated. 

Fig. 5 shows some measured patterns as a function of 
the horn reflector separation. Curve A shows a typical 
guided (or endfire) beam while B represents a reflected 


4. M. T. Jones, “An annular corrugated-surface antenna,” 
Proc. IRE, vol. 40, pp. 721-725; June, 1952. 


e in the antenna system (experimental data). 
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(or broadside) pattern, with the same surface wave ef- 
fect. Further details were published by the author in 
another paper.® 

Within the dimensions of this experimental model, the 
“broadside” part of the reflector can be estimated to ex- 
tend over approximately two wavelengths with the “end- 
fire” section being of about the same dimension when 
considered for the lower end of the frequency band. 
Simultaneously, the beam emanating from the horn is 
rather wide at the lower frequencies, thus illuminating 
the outer “endfire” portion relatively well. In case of 
both parts of the reflector radiating about equal power, 
the over-all beamwidth of this system will be about 30 
per cent smaller than that of a broadside system alone. 
No particular attention was given in the design of the 
model antenna to this surface wave. However, the 3-db 
beamwidth of 22—24° obtained at the low frequencies 
clearly indicates the presence of such a surface or 
traveling wave. As the frequency is increased, the horn 
radiates a narrower beam, decreasing the illumination 
of that part of the reflector from which the end fire radi- 
ation originates. Moreover, the phase relation between 
the reflected wave and the surface wave on the re- 
flector is disturbed, resulting in a very constant gain 
over a wide frequency range (see Fig. 8). 

The separation of the rods forming the reflector is 
widest between the full-length rods just inside the half 
length rods and at the outer edge. At the highest fre- 
quency this separation is about one-third of a wave- 
length. While under these conditions some of the inci- 
dent energy most likely will not be reflected, this ar- 
rangement can be considered to be very simple struc- 
turally and so represent a smaller wind load. 

Two of the more important dimensions could be 
varied in the model. These were the spacing y between 
the horn and the reflector and the length x of the gap in 
the inverting transition. Within the time available, it 
was not possible to determine the exact optima for 
these dimensions, but the dimensions given in Fig. 3 are 
thought to be very close to these optimum values. In 
the course of the experiments, the pattern in the vertical 
plane was found to be little affected by variation of the 
horn-reflector spacing from 0.5 to 1.5 inches. 

The other dimension variable in the model, the length 
« of the gap in the transition, was found to be quite 
critical with reference to good impedance matching, 

Fig. 6(a) to 6(h) shows the measured vertical patterns 
for the model B antenna. Fig. 7 summarizes these pat- 
terns with respect to 3-db beamwidth and the direction 
of the main beam. The following information can be ex- 
tracted from these figures. 


1) The direction of the main beam varies from about 
19° above the horizon to approximately 1°, with 
an average of close to 10°. The design value was 


5 P. Foldes and S. Komlos, “Theoretical and experimental study 
of wide band paraboloid antenna with central reflector feed,” RCA 
Rev., vol. 21, pp. 94-116; March, 1960. 
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Fig. 6—Vertical patterns. 
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5° and a somewhat smaller variation had been 
anticipated. The average can be adjusted simply — 
by increasing the distance y, but a somewhat 
steeper reflector shape will be necessary to reduce 
the variation in the beam tilt. 


2) The 3-db beamwidth shows some fluctuations with 


3) 


frequency, but these are fairly small. Most likely 
these fluctuations came from spurious propagation 
(higher-order modes) and from some interference 
between reflected and surface waves. Improving 
the symmetry of the system mechanically will re- 
duce the possibility for higher mode propagation, 
and the interference problem can be minimized by 
refining the design of the horn and the reflector. 
The high angle radiation from 70° to 90° is rela- 
tively large at the low-frequency end of the band. 
This is the result of the shape of this particular 
reflector. With increasing frequency, the radiation 
from the horn becomes more concentrated towards 
the parabolic section of the reflector and the high 
angle radiation decreases. In the center of the 
band, the radiation is basically in a single beam 
and at the high-frequency end, the radiation 
through the reflector becomes greater and results 
in a relatively high field over the range —70° to 
—90°. Giving the reflector a closer meshed design 
will decrease this radiation. The sidelobe levels at 
+90° are high; however, the energy content in 
these lobes is small in relation to the main beam 
in the case of omnidirectional radiation, and are 
tolerable in many practical cases. 


The gain of this experimental model was determined 
from the pattern measurements by method of integra- 
tion and the result is shown in Fig. 8 by curve B. 
While the 3-db beamwidth varied slightly, the gain 
remains almost constant with frequency. This is con- 
sidered to be due to the high and low angle radiation 


lek ¢ 
- 90-80 70 -60-50 -40 30-20-10 0 10 20 30 40506070 80 90 e° properties discussed before. Fig. 8 shows that the gain 
of this model ranges from 6.7 to 7.1 db, with 7 db being 
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Fig. 8—Gain above isotropic radiator. 


The VSWR behavior of the antenna is depicted in 
Fig. 9. Curve A gives the performance of the system 
with the reflector removed. This curve, therefore, is a 
close indication of the VSWR characteristic of the inver- 
sion transition. a 

Curve B is a plot of the VSWR of the complete sys- 
tem. The addition of the reflector results in a slight in- 
crease only. Maxima VSWR values without and with 
the reflector are 1.45 and 1.68. From this it appears that 
the transition is the important and limiting part with 
respect to impedance characteristic. 

Comparison of VSWR and gain curve reveals an in- 
teresting correspondence. The relatively high VSWR 
at about 2.1 and 2.8 kmc is related to a relatively large 
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3-db beamwidth and corresponding low gain. This indi- 
cates the existence of higher-order spurious modes at 
these frequencies, resulting in higher reflections and im- 
proper illumination of the reflector. 


CONCLUSIONS 


The given simple methods of calculations are satis- 
factory as a design basis of the new antenna. The meas- 
ured experimental data show that the developed system 
is useful for middle gain, high-power UHF applications. 
On the other hand, it should be noted that the given 
experimental results are not the optimum and merely 
indicate the usefulness of the principle. 


Maximally Flat and Quasi-Smooth Sector Beams* 
A. KSIENSKI} 


Summary—Methods for the synthesis of sector beams are pre- 
sented and illustrated with examples. First a perfectly smooth and 
square beam, the equivalent of the maximally flat transfer function 
of network theory, is discussed and synthesized. It is then shown 
that realization of such beams by means of relatively short arrays is 
feasible when a certain amount of oscillation is permitted. The re- 
sulting beam approximates the maximally flat one in an almost equal 
ripple fashion, the amplitude of the ripple depending on the length 
of the array and the desired slope of the main beam. 


* Received by the PGAP, December 8, 1959; revised manuscript 
received, March 1, 1960. The research reported in this paper has been 
sponsored by the Electronics Research Directorate of the Air Force 
Cambridge Research Center, Air Research and Development Com- 
mand, Bedford, Mass, under Contract No. AF 19(604)-3508. 

+ Microwave Res. Lab., Hughes Aircraft Co., Culver City, Calif, 


I. INTRODUCTION 


ERTAIN antenna applications require sector 
beams of constant signal level over a certain 


range of the observation angle ©. Beyond that 
range, the signal level is supposed to decay as fast as 
possible and stay at a very low level for the rest of the 
observation range. The smoothness of the beam top is 
required for certain reconnaissance and search opera- 
tions where the errors introduced by oscillating patterns 
cannot be corrected by the system. Also, certain low- 
noise high-precision antenna applications require beams 
of this type. The permitted deviations from the desired 
signal levels are too small to permit either a Fourier or a 
Woodward type synthesis. What is presently desired is a 
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flat topped beam, similar to the “Butterworth” or “max- 
imally flat” transfer function in network theory. This 


__ paper describes methods for the synthesis of such beams. 
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First, the network and the array problems will be com- 
pared, and it will be shown why the same technique can- 
not be applied to both. 


The Equivalent Problem in Network Theory 


The expression for the maximally flat transfer func- 
tion in network theory is given by 


1 


OS aaraaes 


0O<w< x (1) 


where w is the frequency variable.! In this function, de- 
picted in Fig. 1, the first (2n—1) derivatives vanish at 
the origin which produces the flat top effect for w<1, 
while the high power of w in the denominator causes the 
fast decay of the function for w>1. 


f (w) 


wo —> 


Fig. 1—Maximally flat transfer function. 


Let us now examine the antenna pattern expression 


' representing a half-wavelength? spaced array of 2V-+1 


elements. For a pattern symmetric with respect to the 
origin, the expression is 


N 


F(u) = >> a cos neu (2) 


where u=sin 6. When (2) is compared with the general 
expression representing a network function, 
—N 
arate 
flo) =———_ 0S, (3) 
On b,w” 
0 


it is clear that the most important advantage of (3) is 


that it contains poles, which may be used to produce 


1S, Butterworth, “On the theory of filter amplifiers,” Expert- 
mental Wireless, vol. 7, pp. 536-541; October, 1930. 

2 The inter-element spacing may be varied with well known re- 
sults. See S. Silver, “Microwave antenna theory and design,” M. I. T., 
Rad. Lab. Ser., McGraw-Hill Book Co., Inc., New York, N. Y., 
vol. 12, p. 261; 1949. 
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fast changes of the function at desired points and pro- 
duce zeros for w>«. This results in the high attenua- 
tion of f(w) for w>1 and contributes to the squareness of 
the Butterworth function. 

Since in the antenna array case there are no attenua- 
tion factors corresponding to poles, one is forced to di- 
vide the available parameters between the origin and 
the end of the visible range. An examination of (2) will 
show that there are V+1 parameters, and they can be 
divided in any desired proportion between the two end 
points (4=0 and u=1). The present goal is to obtain a 
pattern that has a flat-topped beam which attenuates 
smoothly to zero at the highest possible rate for the 
given number of elements in the array. This will be 
considered as the equivalent for antenna arrays of the 
“maximally flat” function. 


Il. MAxIMALLY FLAT BEAMS 
General Discussion 


The necessary constraints may be imposed directly on 
(2), and the resulting a, computed. Such a solution 
would involve a matrix inversion, which may be quite 
cumbersome. What is more important, it is difficult to 
evaluate the behavior of a trigonometric polynomial for 
incremental values of @ or to evaluate other functional 
properties. Ordinary polynomials can be handled more 
easily since their properties are simpler in many respects 
and can be directly evaluated. Therefore the trigono- 
metric polynomial will be converted into a simple 
polynomial form and the various operations and evalua- 
tions will be carried out in that form. In the last step, 
however, the polynomial is reconverted to a trigonomet- 
ric form to avoid slow convergence. As will be shown 
later, a polynomial of a very high order is required to 
obtain a satisfactory square beam. However, when the 
polynomial is transformed into its trigonometric repre- 
sentation, it is found that the series is so highly con- 
vergent that only the first few terms are significant. 

Using this technique, rewriting (2) in powers of cos ru, 
and substituting « =cos ru gives® 


Pe), An (4) 


Since the derivative behavior of the pattern is of in- 
terest, the following expression will be discussed: 


F'(x) = > Byx. (5) 


F(x) will then be obtained by integrating the desired 
Bxye 


3 In the present discussion an array of an odd number of ele- 
ments is assumed. The even number case is discussed in Section IAW 
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Since the derivatives of the pattern are to vanish at 
the end points only, zeros of appropriate orders are 
placed there to obtain a desired beamwidth. Thus the 
general expression for the pattern derivative is given by 


Ba) = (a Ie —1)F (6) 


where a+@=N-—1, and where the ratio of a/f is 
chosen to produce a given beamwidth. 

The choice of a particular beamwidth, +30°, will 
simplify calculations and aid in the determination of 
the characteristics of such beams. Most of the details 
worked out for this example will apply to other beam- 
widths. 


Beams of +30° 


If the zeros of the derivative are equally divided be- 
tween the origin and the end of the visible range, a 30° 
beamwidth is obtained.t This can be seen from (6); 
letting a=B =m yields 


P(x) = (@ = 1)”, (7) 


which is symmetric with respect to the origin. Thus 
the angle 6 corresponding to the value x=0 repre- 
sents the beamwidth. To make sure that the desired pat- 
tern is obtained, the proper sign must be attached to 
make F’(x«) non-negative for | x| <1. Thus the correct 
expression is given by 


Fa) = 1 — 42). (8) 


F(x) is then given by 
Fic) = f (1 — x?)™dx + Ao, (9) 
0 


where A, is reserved to adjust the level of F(x) so that 


Rei = 6 (o- +7), (10) 


Example 


Let us compute the excitation distribution and the 
pattern for N=7 which corresponds to m=3. Applying 
(9) and (10) we obtain, after converting to cos nru and 
expanding, 


F(u) = 0.457 + 0.547 cos ru — 0.109 cos 3ru 


+ 0.0215 cos 5ru — 0.002 cos 7rru. (11) 
This pattern is plotted in Fig. 2 and compared to a 
Fourier approximation of the same desired sector beam 
of + 30° width. In both methods the same number of ele- 
ments was used. Note the smoothness of the maximally 
flat beam and the perfect agreement with the desired 


* Beamwidth is defined as the angular distance between the points 
of maximum beam slope which in this discussion are equivalent to 
6-db points. 
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Fig. 2—Sector beam produced by a 7-wavelength array. 


pattern for the first +15°. This, however, is obtained 
at the expense of beam slope. 

The slope at the transition from the main beam to the 
sidelobe region and its interaction with the flatness re-- 
quirements will now be discussed. The slope® at the 
transition point x=0 may be specified by choosing an 
appropriate value for ¢ in 


F(x) = c(1 — 27), (12) 
and the desired flatness may be specified in terms of 
maximum allowable value of F’(x) for | x| larger than 
some desired value. This will determine m, that is, the 
length of the array. If the specification is given in terms 
of allowable deviation of the main beam from the con- 
stant signal level, the following approximation may be 
used to simplify computation. The behavior of F’(x) in 
the vicinity of x =-+1 is affected primarily by the mth- 
order zero at x= +1 and may therefore be approximated 
by 


"(ey =a cu INGE ee ae se (13) 


where the center of coordinates was shifted to x =1, and 
where the effect of the zeros at « = —1 is represented by 
a”. For an allowable deviation A, we obtain 


1/2 
e(e).= f GO™L™ de <n, (14) 
0 
For a=2 the upper limit for m is obtained, 
g 
m <— — 1, (15) 


Now, a specification which is most likely to appear is the 
desired interval in x over which the beam has to decay, 
say to 10 per cent of its value. This is the equivalent of 
the decay time in networks. For this purpose, the follow- 


_ > The actual slope with respect to 6, the observation angle, is 
given by : 


aa = F’(«)m sin (au) cos 6. 
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ing simple approximation may be used. Since the 
change considered is for small values of ~, 
(1 — x2)™~ em” for x <1, (16) 


Hence, given the value of x» for which the beam should 


_ decay to 10 per cent of its initial value, or should reach 


90 per cent of its final value, the following relation may 
be employed: 


eo . x9 
a3 2 
i} € max — [ 6, dx 
0 0 
oO 
if emma de 
0 


Substituting (to convert to the standard form of the 


= 0.1, (17) 


- error function) 


t 
6 18 
re (18) 
and normalizing results in 
v zyvm a 
=| e *dt = — = 0.9, (19) 
Vr 9 b 


where a/b represents the fraction of the final value to be 
reached. Assuming that x)=0.1, inserting this value 
into (19), and using error function tables® results in 
m = 134.5. Thus, given a desired “rise time” to any frac- 


- tion of the final value, the required value of m can be 
~ computed. 


_ Beams of Arbitrary Width 


Maximally flat beams of varying width will now be 
derived and their characteristics evaluated. It will be 
shown that the beams closely approach a square shape 
for increasing array lengths. 

An appropriate change is made in the ratio of the 
constants a and # appearing in (6), so that the beam- 
width may be varied without introducing oscillations. 
The transition point, which is defined as the point of 
maximum slope, moves from x =0 to the point at which 
F’’(x) vanishes. Changing the exponent notation from 
(a, 8) to (n, an) results in 


F'(«) =0=[(1 — x1 + a)")’, (20) 
where 
mn+an=WN —1. 
The solution of (20) yields 
ee a-—l (21) 
atl 


6L. J. Comrie, “Chamber’s Six Figure Mathematical Tables,” 
D. Van aad Co., Inc., New York, N. Ve vol, 2, py018; 1949. 
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Thus when a is adjusted from zero to a large value, the 
beamwidth changes from a very wide and flat beam to 
a narrow main beam which is flat in the sidelobe region. 
The first beam is obtained by letting a=0; all param- 
eters are used to produce zeros of the derivative at the 
origin which yields 


F(x) = (1 — #)*-1, (22) 


This form may be useful where an approximation to 
an omnidirectional pattern is needed. The second beam 
is obtained when a is large and is given by 


F(x) = (1+ «)N-1, (23) 


This is better known as the binomial pattern, which is 
obtained when all parameters are used to produce zeros 
at the end of the visible range, and it has the property 
of being the narrowest possible beam having no side- 
lobes. 

An additional means of beamwidth variation is to 
vary interelement spacing. Thus, if no scanning is ex- 
pected, one may contract the main beam up to the 
point where secondary beams begin to appear. On the 
other hand one may widen the beam as much as de- 
sired, subject only to the minimum practical interele- 
ment spacing. 


Transition Slope and Beam Flatness as a Function of 
Array Length 


The beam height must be normalized to unity at the 
origin so that the true shape of the beam in terms of 
maximum slope will be discernible, as well as the slope 
variation throughout the pattern. 

The point of maximum slope is given by (21). The 
slope at that point is 


a—1 2 ar ie 
Mer el) ek 
a+ 1 a+1 
In order to normalize the above expression, the beam 
height is computed, 


(24) 


I 


+41 
f (1 + x)"(1 — x)o"dx 


—1 


H(a, n) 


= 2le+)nt1B(y + 1,an+ 1) (25) 


where B is the Beta function. Since am is an integer, 
factorials may be used, and 
Qn(atl+l 


n\(an)! 


Pee ntiaie (ee fal 


H(a, n) = (26) 


The resultant expression for maximum slope is given by 


a—l 
(== -) ” [(a+1)n+ 1] [a+ alt ae (27) 
H(a, n) 


‘4 ati) ni(an)! 
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As n and an increase, Stirling’s formula for the factorial 
may be used, yielding 


a—1l 
e(<—) 
a+ 1 1 


H(a, n) 2 


(a+ 1)n+1 
VJ2an 


{/* + 1 (28) 
a 
This indicates that the slope is proportional to Wn for 
fixed a. Then, as afunction of a, it again grows for both 
small a, where the denominator takes over, and for 
larger a where the numerator contributes to the growth. 
For all points other than the break point, (a—1/a-+1), 
the slope tends to zero, and the rate of decay is propor- 
tional to the power of x. 
Let A denote the increment from the point of maxi- 
mum slope. It can be shown that 


(Ss -+4) 
a+1 atl 
~n| tog (1 - ; 4) 


H (a, n) 
+ atog (1 +=" a)], (29) 


log 


plus terms not containing the factor z. Assuming that A 
is a small increment and keeping the terms up to A’ 
yields the approximation 


Be ia.) 
re lace 


J 


which is negative for all positive values of a. Hence, 


’ Gy = I 
lim F’ 4) >0, 


ye Il 
From (28) and (30), it is seen that the beams converge 
to a rectangular shape as the array length increases, 
and that this holds for any beamwidth. 
Since most of the detailed analysis presented for the 
+30° beam can easily be extended to beams having 
other widths, no further elaboration is necessary. 


(31) 


n> 0 


III. Quasi-SMootH SEctoR BEAMS 


In Section II the mathematical model was developed 
for maximally flat beams. Exact realization of the pat- 
terns considered would require very long arrays for a 
fairly square sector beam. This is the price that must be 
paid for a perfectly smooth beam. However, in practice 
a perfectly smooth beam cannot be distinguished from 
a beam whose oscillations are of an amplitude which 
cannot be detected by the measuring apparatus. Hence, 
the absolutely smooth pattern may be replaced with a 
relatively smooth pattern having oscillations sufficiently 
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small in amplitude. This is precisely where the conver- 
sion back into the trigonometric form is most im- 


portant. The rapid convergence of the trigonometric 


series permits omission of a large part of its terms, and 
the synthesis technique is consequently economical in 
number of elements. 

A more appropriate formula is first developed for ob- 
taining the excitation coefficients. Eq. (9) can be re- 
written as 


cos ¢ 
r@) =f erat 


N 
= De a, cos kp (32) 


kat 


(k odd), 


where Ao was dropped and a beam of +30° is again 
used as an example. 

Taking a derivative with respect to @ under the in- 
tegral sign, expanding in terms of sin ¢, and integrating, 
one obtains the following series :7 


(2n + 1)! 1 (2n + 1)! 
F Sr Es CL Se — 
iL. acta eee 
Bere dee Dy 2 : 
+ + Salas pa! n +1). (34) 


Since the absolute values of the coefficients are of little 
importance, only the relative magnitude is required, 
and this relationship is given by 


2n—r)+1 r 
dy Tar) ieee 


ay—1 


(34) 


Thus the first coefficient, that of cos ¢, is assumed to 
be unity, and the others are obtained by means of (34). 
Also, from the same equation, an upper bound can be 
obtained for the error involved in truncating the series 
after a given term. For example, when r=27/3, 


; ( *) ec 
Qy-1 3 


a, o(*) +3 4 
3 ee 


et 
3 


cide 3 
a tees 


Thus all terms following @2n/3 will be multiplied by 
factors smaller than one half. This implies that a ge- 
ometric series whose common ratio is g=4 has terms 
larger than the series under consideration, and thus can 
be used as an upper bound for convergence and error 
estimates. 


’ For details, see A, Ksienski “Maximally-Flat and Quasi-Smooth 
Beams,” Hughes Aircraft Co., Cul er City, Cali ienti er 
No: 3084s atoee Vv ity, Calif., Scientific Rept. 
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Exam ple 


Consider = 50 and compute the first ten coefficients.® 
The resulting expression is given by 


F(u) = Ao+ cos ru — 0.321 cos 3ru + 0.178 cos Sru 
— 0.113 cos 7ru + 0.075 cos 9ru 
— 0.050 cos 117u + 0.034 cos 137u 
— 0.022 cos 15ru + 0.014 cos 177ru 


— 0.009 cos 197u, (36) 


When (36) is plotted as shown in Fig. 3, it is seen that 
the actual error or oscillation amplitude is equal to the 
omitted term 0.006 rather than 0.016, as the estimate 
would indicate. The reason for that is the alternating 
nature of the series, which decreases the significance of 
the omitted portion. Also note that, because of dom- 
inance of the first omitted term over the following ones, 
the pattern has an almost equal ripple characteristic. 
The above discussion was concerned with +30° 
beams, but the results indicate the behavior for other 
beamwidths. The main difference is in the computation 
of the excitation coefficients. Using the expression cor- 
responding to (32), and converting to cos ¢, we obtain 


cos } 
FQ) =f +o — pera, (37) 
x /2 
which, upon substituting 20=¢, yields 
F' (6) = 2(at2)n41 sin 241 9 cos?”+1 9 
2N 
= >> by sin k6 (k even). (38) 
k=2 


The coefficients b, may be cbtained either by expan- 
sion or by Fourier sine inversion 
Tv 
hb =A f sin22"+1 9 cos?"+! @ sin k0d8, (39) 
where A is a constant common to all coefficients. The 


resulting series is integrated as before, after @ is sub- 
stituted for 2 in (38). 


IV. SynTHEsIs BY MEANS OF ARRAYS OF AN EvEN 
NUMBER OF ELEMENTS 


A symmetric pattern produced by an array of 2N 
elements spaced a half wavelength apart 1s given by 


F(u) = 3 dn COS | oh = ae H (40) 


This may be transformed into a polynomial given by 


F(x) = 5 Agxer}, (41) 


8 Ay is chosen as to make F (u) vanish at u=1. 
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Fig. 3—Quasi-smooth sector beam produced by a 
19-wavelength array (21 elements). 


where x=cos(ru/2). Taking the derivative of the 
above expression, we obtain 
N 
F’(x) = >> Bnx?”2, (42) 
1 
To obtain a maximally flat beam the zeros have to be 
placed at x=1, and x=0 corresponding to u=0 and 
u= +1. The expression, therefore, for the derivative is 
given by 


F(a) = c#*e(1 —ix7)8, (43) 


where a+6=N-—1, as in the previous case of an odd 
number of elements. The reason that (43) is given in 
even powers is, of course, to conform to (42), which con- 
tains only even powers. This formulation does not: dis- 
tort the pattern in any way, since all the zeros of the 
derivative appear at the correct points, namely at 
x=0, 1. The zero at x= —1 corresponds to u=2, which 
expresses the behavior of the pattern in the invisible 
range where similar beams are formed above and below 
the u axis. The details of analysis are parallel to the 
ones concerning arrays of odd number of elements and, 
therefore, will not be discussed further. 


V. Sector BEAMS DERIVED FROM 
OptiMuUM PATTERNS 


The synthesis method presented in this section has 
the advantage of extreme simplicity and practicality. 
It is obtained almost directly from a Tchebycheff pat- 
tern, or similar pencil beam patterns, and can draw on 
a large backlog of information and tabulated data 
available for these patterns. It can thus be very useful 
for the antenna engineer as a source of readily accessible 
design information or as a means of estimating the per- 
formance of an array without elaborate computation. 

The synthesis procedure is based on the reduction of 
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a desired pattern into a combination of simpler pattern 
functions of forms that are either easier to synthesize 
or that may have been already synthesized and opti- 
mized.? This procedure may involve differential and 
integral operators, operating on elementary pattern 
functions to produce the desired shaped beam. Thus, 
all operations are carried out in the u=sin @ domain 
until the final results are obtained. Then the cor- 
responding operations are carried out on the excitation 
distributions; this procedure is permissible since the 
distributions are themselves connected by a linear op- 
erator to the patterns, namely a Fourier transform. In 
many cases the required operations may be rather 
simple, as for example the Woodward synthesis, which 
involved only the addition of elementary patterns. An- 
other simple operation is the integration or differentia- 
tion of an elementary pattern of any Tchebycheff or 
sin x/x form that may yield a desired pattern. 

In the present case of a sector beam, the derivative 
with respect to u=sin 6 yields two delta functions dis- 
placed from each other and of opposite sign (see Fig. 4). 
Now the required operation is quite clear; two pencil 
beams of Tchebycheff, Taylor, or other appropriate 
nature are directed to +4, respectively, and sub- 
tracted from each other. The resultant pattern is then 
integrated, yielding the desired approximation to a 
sector beam. 

The excitation distribution is derived by means of 
the corresponding operations carried out on the original 
pencil beam distributions, as will be demonstrated in 
the following example. The patterns obtained by such 
operations will, of course, not be of optimum nature 
themselves, even if derived from optimized beams. How- 
ever, as the results indicate they are of very satisfactory 
form, the errors can be easily estimated, and, what is 
rather important, very few calculations are necessary. 

To provide a fair comparison, the sector beam syn- 
thesized here will be, by means of an array of the same 
number of elements, as in Section III. Let T(u) repre- 
sent the appropriate Tchebycheff pattern given by!° 


N 
T(u) = >> A, cos nru. 


n=0 


(44) 


If uo represents the desired beamwidth of the sector 
beam, the pencil beam combination will be 


N N 
F’(u) = >> A, cos na(u — uo) — >. A, cos nm(u + up) 
n=0 


n=0 


N 
= 0 A, sin (nmruo) sin NTU, 


n=0 


(45) 


POR Price and R. F. Hyneman have independently conceived 

a similar idea which they applied to the synthesis of a monopulse 

pattern. See “Distribution Functions for Monopulse Antenna Dif- 

ae Patterns,” Hughes Res. Lab., Malibu, Calif., Res. Rept. 
0. : 

© This expression assumes an odd number of elements; however 

for an even number of elements, the procedure is identical, 
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Fig. 4—A sector beam and its derivative. 


from which we obtain the desired pattern 


uf JAG, 
Fa) = SO Sein (nruo) Cos nu, (46) 
n 


n=0 


where the A, are the original excitation coefficients of 
the Tchebycheff pattern.!! For the case of up = 4 we ob- 
tain 


n—1 


Btu) = Ay ye oe cn) 


111 MW 2 


cosuru (nodd), 


(47) 


in which A is chosen to make F(u) vanish at u=1. In- 
troducing the proper coefficients obtained from a 38-ele- 
ment, 40-db Tchebycheff, we obtain 


F(u) = Ay + cos ru — 0.319 cos 3ru + 0.176 cos 5ru 
— 0.116 cos 7ru + 0.071 cos 9ru 
— 0.046 cos 117% + 0.029 cos 133% 
— 0.017 cos 15ru + 0.009 cos 177u 
— 0.008 cos 197. (48) 


"L. B. Brown and G. A. Sharp, “Tchebycheff Antenna Distribu- 
tion, Beamwidth, and Gain Tables,” Naval Ordan Kaba 
Calif., NAVORD Rept. No, 4629; 1958, a i eee 
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A comparison of (36) shows the closeness of the two re- 
sults. F(u) is plotted in Fig. 5, which shows that we 
have an equal ripple for almost the whole range of u, 


indicating that the pattern is close to the optimum “rise 


_ time” for a given oscillation amplitude. 


a PE PALS 


> 
—< 


* ue * 2 


ATS 


es. 3 


5 


We shall now investigate the relation between the 


_ main beam to sidelobe ratio of the derived sector beam 
_ and the corresponding ratio of the original pencil beam. 


Since the sector beam represents the area under the 
pencil beam, the sidelobe level is not directly related to 


_ the sidelobe peak values of the original pattern, but 


rather to the area under those lobes. Thus in the region 
where the oscillation of the Tchebycheff pattern is rapid 
the resulting sidelobes are smaller, and when the oscilla- 


_ tions are relatively slow, the sidelobe level will increase. 
_ When we have the summationof two (or more) beams, the 
oscillation rates will be different in the two beams for 
the same value of uw, and the resulting sidelobes will 
oscillate at a higher rate and with higher amplitudes. 
- This increased rate will tend to reduce the integrated 


sidelobe level from the value which would have resulted 


if the sum of the peaks carried over directly. The side- 


lobe level can be computed in approximate fashion by 
the use of the ideal space factor” cos (4+/u?— A?) as the 
pencil beam, and the resulting ratio is multiplied by a 


_ factor of 2 to compensate for the increase in sidelobe 
level due to sidelobe combinations and of course for 
_ the fact that the distribution is not ideal. Since the area 


of the slowest sidelobe is of interest, one may assume 


that u>>A and the sidelobe area is given by 


ar/2 
A A f cos m+/u2 — A? du 
—1r/2 


w/2 » 
f cos rudu = — - 


—r/2 Tv 


(49) 


The area under the main beam represents the height of 


the sector beam and is given by 
A ee 
Ams. f cos mv/u2 — A* du. (50) 
—A 
- Substituting 
u = Asin, 
~ we obtain 
F aw /2 
Aes Af cos (ir A cos ¢) cos ¢do 
—1r/2 
= ArQi(ri A), (51) 


2 T. T. Taylor, “Design of line source antennas for narrow beam- 


3 width and low sidelobes,” IRE TRANS. ON ANTENNAS AND PROPAGA- 


Be tae TS 


TION, vol. AP-3, pp. 16-28; January, 1955. Note that in this paper 


= is the aperture width. Although the space 
Be donot ‘ i be used for a discrete 


- factor refers to a continuous aperture, 1t may 


array with negligible error. 
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Fig. 5—Tcheybcheff-derived sector beam produced by a 
19-wavelength array (21 elements). 


where 2),(Z) is the Lommel-Weber function of mth order. 
If a table of these functions is not readily available, the 
following series expansion may be used:¥ 


Ax? A*ir* 


oe it @— 1)(42 — 1) 


| (52) 


This series converges quite rapidly. Applying these re- 
sults to our example where a 40-db sidelobe ratio was 
considered, the value of A is given by 5.3, and the cor- 
responding value for (52) is 50.37. The main beam to 
sidelobe ratio is 


G5 
ai) = E a 


A476 
e (= hie = ee) or 


ASE. AQ, (ir A) 
At 2/1 


= 266. (53) 


The actual sidelobe level is at a ratio of 1:158 to the 
main beam, which indicates that the sidelobes have not 
added directly; direct addition would have resulted in 
a ratio of 1:133. Thus, the sidelobe level has been re- 
duced somewhat. Actually the sidelobe level may be 
further reduced by a slight alteration of the shift angle 
uy, which determines the way in which the sidelobes of 
the pencil beams combine. This is particularly effective 
in arrays of a small number of elements where most of 
the sidelobes are spaced fairly uniformly and, by a 
proper shift, the sidelobes of the two patterns may be 
made to cancel each other almost completely. 

The following relations between the Tchebycheff and 
sector beams may be used. Given a desired “rise time” 
for the sector beam, a Tchebycheff beam of the cor 


18 E, Jahnke and F. Ernde, “Tables of Functions,” Dover Pub- 
lications, New York, N. Y., p. 211; 1945. 
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responding beamwidth is chosen. Thus, for’ example, if 
the sector beam is to reach to within 10 per cent of its 
final value within an interval 0, this interval is chosen as 
the beamwidth for the “parent” Tchebycheff and where 
the beamwidth is presently defined as between 10-db 
points. One can easily estimate the area under the pen- 
cil beam to determine what portion of it constitutes 90 
per cent, or any other percentage specified, and thus 
obtain the appropriate interval 6 defining the beam- 
width. Once the beamwidth is determined, the array 
length is obtained from the permitted oscillation ampli- 
tude or sidelobe ratio. 


VI. SYNTHESIS BY MEANS OF MODIFIFIED 
FOURIER TECHNIQUES 


In order to provide a more complete exposition of 
synthesis methods applicable to sector beam synthesis, 
alternative methods will be briefly reviewed. If one de- 
sires only to eliminate or reduce the ripple contained in 
a Fourier synthesized sector beam, there are several 
methods which will accomplish that task. These meth- 
ods apply a taper to the Fourier coefficients, and thus 
insure a smooth convergence to the desired pattern. The 
Fejer summation technique" applies a linear taper given 


by 
N-—-n+1 


A’, = 
N 


A, (n21), (54) 


where A, is the mth Fourier coefficient, and the A’, 
is the modified one, and where the series has N+1 
terms. The results of such a taper are shown in Fig. 6 
and compared to the sector beam synthesized earlier. 
It is seen that the ripples completely disappeared. 
However, the beam shape suffered substantially. One 
may apply a more sophisticated taper or summation 
scheme, called the “Cesaro m” summation, in which the 
weighting factor is given by 


(N — 1)(N —n +m)! 
(N — n)\(N +m — 1)! 


(55) 


By varying m between zero and unity, various tapers 
between the original Fourier (m=0) to the Fejer 
(m= 1) can be obtained, thus trading between the ripple 
content and beam squareness. The disadvantage in the 
above techniques is that one has no way of predicting 
the oscillation amplitude for a given taper or the result- 
ing beam steepness; thus, the method involves a trial 
and error procedure to obtain the desired beam shape. 


4 E. A. Guillemin, “Mathematics of Circuit Analysis,” John 
Wiley and Sons, Inc., New York, N. Y., p. 496; 1949, . 
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Fig. 6—Comparison between a Fejer and Tchebycheff 
derived sector beam. 


VII. CoNcLUSIONS 


We have shown how a radiation pattern equivalent 
in nature to the maximally flat transfer function of net- 
work theory can be synthesized. This implies that the 
sector beam synthesized approaches a square beam 
without any oscillation (or sidelobes) occurring through- 
out the whole range of the observation variable. The 
“squareness” of the beam, 7.e., its flatness and rate of 
cutoff, is related to the length of the array, and will ap- 
proach the perfect shape as closely as desired for in- 
creasing array lengths. Since the array lengths required 
for perfectly smooth beams of reasonable flatness and 
cutoff rate are rather large, and since perfect smooth- 
ness is undetectable, a more practical approach is con- 
sidered. It is shown that the series representing the ex- 
citation coefficients converges so rapidly that only the 
first few terms are needed to obtain the desired beam 
with almost negligible oscillations, and only a relatively 
short array is then required. Because of the nature of 
the series, the first neglected term approximates the 
oscillation level. 


The last part of the paper presents an alternative 


synthesis technique for obtaining, in a rather simple 
fashion, a sector beam of desired oscillation level and 
cutoff rate. This technique uses pencil beams as “par- 
ent” patterns, from which the sector beam is derived by 
integration. The sidelobe level of the resulting beam is 
indirectly related to the original pattern and is, in fact, 
somewhat reduced. The main advantage of the method 
is that most of the tabulated data relating to Tcheby- 
cheff and Taylor patterns can be almost directly applied, 
and the results approach an equal ripple approximation. 
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A Geometrical Optics Method of Pattern 
Synthesis for Linear Arrays* 


HOWARD E. SHANKST 


Summary—A unique phase expression is given which allows 
synthesis of a wide variety of shaped radiation patterns for a linear 


_ array. The accompanying amplitude distribution is found by a sta- 


tionary phase evaluation of the radiation integral and is shown to be 
functionally related to the desired radiation pattern. Because of the 
optical type of approximation (stationary phase) used to evaluate the 


radiation integral, the method is most applicable to antennas with 


large size-to-wavelength ratios. 


INTRODUCTION 


HE problem of pattern synthesis from array type 
Blegic. has received considerable attention for a 
~~ number of years. Dolph! and Taylor,? among 
others, have contributed greatly to the understanding of 


- directive pattern synthesis from linear arrays. In addi- 


tion, the work of Woodward’ and Dunbar‘ has estab- 
lished methods by which patterns with shaped char- 
acteristics can be synthesized. Although this body of in- 
formation has many applications, situations arise for 
which these methods are either not physically realiza- 
ble or not of optimum utility. 

This paper describes a technique for shaped beam 
synthesis which, although related to that outlined by 
Dunbar, is applicable to a variety of practical antenna 


- configurations. Both schemes utilize a stationary phase 
~ evaluation of the radiation integral and are therefore 


equivalent from the standpoint of ultimate accuracy. 
The present theory differs from Dunbar’s in that a 
judicious choice of the phase distribution is made and 


the amplitude excitation is adjusted to produce a de- 


sired radiation pattern. Moreover, once the range of 
synthesis is fixed, the phase expression remains un- 
changed for any type of pattern; the pattern under this 
condition is controlled only by the amplitude distribu- 
tion. The required phase distribution is a simply vary- 


_ ing function, and hence the present method may be of 


considerable value as an alternate technique or as the 
preferred technique in special situations. 


* Received by the PGAP, October 16, 1959; revised manuscript 
received, February 22, 1960. This work was supported by the USAF 
under Contract No. AF19(604)-3508 with the Air Force Cambridge 
Reserve Center. Similar material has been published as Scientific 
Rept. 3508/1 on this contract. : ‘ 

+ American Systems, Inc., Inglewood, Calif. This work was per- 
formed while the author was with the Hughes Aircraft Co., Culver 
City, Calif. . ; 

1C. L. Dolph, “A current distribution for broadside arrays which 
optimizes the relationship between beamwidth and sidelobe level, 
Proc. IRE, vol. 34, pp. 335-348; June, 1946. ; 

2'T. T. Taylor, “Design of Line Sources for Narrow Beamwidth 
and Low Sidelobes,” Hughes Aircraft Co., Culver City, Calif., Tech. 
Memo. No. 316; July 31, 1953. ge? 

3P. M. Woodward, “A method of calculating. the field over a 
plane aperture required to produce a given polar diagram,” J. IEE, 
vol. 93, 1946. Pk: 

‘4A S, Dunbar, “On the theory of beam shaping,” J. Appl. Phys., 


vol. 23, pp. 847-853; August, 1952. 


THEORY 


It is common in the analysis of many types of arrays 
(particularly arrays of closely spaced elements) to ap- 
proximate the array excitation by a continuous distribu- 
tion. This is assumed in the present theory and results 
in little loss of generality since most practical arrays 
with discrete distributions approximate the continuous 
case rather well. The array orientation and geometry 
to be used in the analysis are illustrated in Fig. 1. 


x24 


x20 


Fig. 1—Linear array geometry. 


The integral representation for the radiation pattern 
from such a line source is given by 
I 
g(0) = A@ W(x) ef? @ eike sin Ody (1) 
0 
where A(6) is a factor related to the polarization of the 
array, and y(x) and (x) are the amplitude and phase 
distributions, respectively. The element factor is not 
essential to the synthesis procedure and may be elim- 
inated by redefining f(@) =g(@)/A(@). Under this condi- 
tion (1) becomes 
U 
f(8) = p(x) et? @ este sin Ody, (2) 
0 
This expression can be considered the basic equation 
for synthesis from line sources of either the discrete or 
continuous variety. The synthesis problem is then one 
of determining W(x) and ®(x) so as to produce a pre- 
scribed pattern f(6). 

Because of the problems inherent in synthesizing a 
given pattern over the entire real spatial domain, it is 
common to specify a restricted angular domain as the 
region of interest. For definiteness, it will be assumed 
here that the pattern is specified in the range 4o <6 <0 nn 
As a result, the procedure will allow synthesis of a given 
pattern only within this range and will provide no con- 
trol outside of this region. 
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A unique phase expression utilizing this angular 
range, which is the key to the entire synthesis proce- 
dure, may be written: 


Rl Am — 9o 
cos i x + Oo |. 


m 0 


(3) 


&(x%) = : 


It should be noted that this expression is a function only 
of the array length in wavelengths and the angular 
range. For the present, this relation is presented without 
justification. However, a plausible derivation based on 
intuitive considerations is given in the Appendix. Its 
utility will become apparent in the stationary phase 
evaluation of the radiation integral. 

When (3) is substituted into (2), the radiation integral 


becomes 
l Om — Bo 
cos | x + A 
6 i 


tn = 0 


10) = f Va) exp jk 


+ x sin ab dx. (4) 


Considerable simplification is realized when the array 
length is normalized to unity by letting n=x/l. Eq. (4) 
then becomes 
cos [(@m — 00) + | 

(Am ee 90) 


f@) = if TONS jn} 


+ 7 sin sh dn. (5) 


The general synthesis result is obtained when the 
integral of (5) is evaluated by stationary phase meth- 
ods. The stationary phase method of integral evaluation 
is outlined by 


ate 
i h(x) 78) da 


(er={5 


[pees] Ho enfmo st] +0(2) 0 


where u’(a) =0; that is, a is the value of « for which the 
phase is stationary. When the first term in this expan- 
sion is used to approximate the integral, an error of the 
order of 1/8 is made. This approximation, then, be- 
comes better as the value of 6 increases. In the present 
formulation [see (5)]| 6B corresponds to ki, the array 
length in wavelengths. The consequences of this error 
in terms of the array length will be mentioned later. 

Applying (6) to (5) results in the integrated radiation 
pattern, 


f(0) =1 (sect 
ieee = bs 
wi jut] igs oe sino | (7) 
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where the stationary phase point ns is given by 
0 — 6 


ES 


Since, in most practical antennas, only the magnitude 
of the field is of interest, the phase term is unimportant 
and may be neglected. In addition, the proportionality 
constants are not essential and can be neglected without 
redefining f(@). Under these conditions, 


| £0) | = ¥(—* ) sects 
lite Alain 


Close examination of (8) reveals that if the argument 
of w is chosen as a special function of 7, namely 


v(m) = v(m cs 90)” rT, ol, 
then (8) reduces to 
| (0)| = ¥(@) sect/2 9, 


This relation may be solved to obtain the required 
amplitude distribution y in terms of the pattern f(@). 
When this is done, 


¥(0) = |f(@)| cost’? @. 


With the result from (10), the complete aperture ex- 
citation necessary to synthesize a desired pattern [(6) 
is given by 


Fm — 40)n + 80] cos*!? [(8m — 0) + 80] 


(8) 
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It should be noted that once the length of the array 
(kl) and the range 0)<0@<6,, are fixed, the pattern is 
controlled only through the amplitude distribution. As 
mentioned earlier, the stationary phase evaluation of 
the integral leading to (11) involves an approximation 
based on neglecting a term of the order of 1/k/. The ac- 
curacy of the synthesis procedure can thus be expected 
to increase as the array length (k/) increases. This fact 
is verified in the following section where several pat- 
terns are computed exactly, using the synthesized dis- 
tributions. 


NUMERICAL COMPUTATIONS AND PATTERNS 


The accuracy and utility of the present technique 
are best illustrated by graphical comparisons of the de- 
sired and synthesized radiation patterns. This is ac- 
complished by an exact numerical evaluation of the 
radiation integral using the distribution obtained by 
the synthesis procedure. This calculation has been car- 
ried out for a cosecant-squared pattern and a sector 
type pattern using various values of the synthesis pa- 
rameters k/, 0m, and 4). As an example of the calculation, 
consider a sector pattern which is to be uniform over the 


range —45°<@<45° for an array with kl=100. The 


mi 
-. 
. 


a 
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radiation integral under this condition is 


0 A erdetal 


-exp [7100p sin 6]dn, 


200 
€xp7—— cos 
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__ which is readily integrated by numerical methods. Table 
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J is a compilation of the patterns calculated with the 
values of the parameters used in each case. 


TABLE I 
Desired Pattern kl Oo Om, 
Sector 100 —45° 45° 
Sector 1000 —45° 45° 
csc? 100 10° 90° 
csc? 1000 10° 90° 
csc? 100 10° 80° 
csc? 1000 10° 80° 
csc? 100 107 102 
csc? 1000 103 70° 


These computed patterns have been graphed and are 
displayed along with the desired patterns in Figs. 2-5. 
The results for the sector beam pattern shown in Fig. 2 
deserve special mention. For k/=1000 an extremely 


- close approximation to the desired form is realized, and 


the marked increase in accuracy for k/ increasing from 
100 to 1000 is evident. 

The cosecant-squared results of Figs. 3-5 demon- 
strate the increase in accuracy with increasing array 
length and also provide an indication of the effect of 

angular range on the procedure. The closest approxima- 
tion to the desired csc? pattern is obtained for the larg- 
est range (10°<0<90°). This might have been expected 
from consideration of the synthesized distribution as the 
range approaches zero (8m —8—0). In this limit a uni- 


form distribution is obtained, giving a pattern with 


many peaks and nulls. 
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COMPARISON WITH OTHER SYNTHESIS PROCEDURES 


It is of interest to compare the present technique 
with conventional synthesis methods, in particular 
those of Dunbart and Woodward.? One way of ac- 
complishing this comparison is by computation of the 
exact radiation pattern from a distribution obtained 
from these synthesis methods; this is done for Dunbar’s 
method. In addition, it will be shown mathematically 
that the inherent accuracies of all methods are similar. 
In the first case it will suffice to show that the present 
distribution is in fact mathematically compatible with 
Dunbar’s. For comparison with Woodward’s technique, 
it will be demonstrated that both methods are exact in 
the limit of array length increasing to infinity. 

The basic equations of Dunbar’s synthesis in the 
present notation are 


d®(m) _ 
dn 


kl sin [6(n) | (12) 
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Fig. 3—Desired and calculated csc? patterns; 0.=10°, @m=90°. 
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Fig. 4—Desired and calculated csc? patterns; 09 =10°, 4, =80°. 
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Fig. 5—Desired and calculated csc? patterns; 6)=10°, 0m =70°. 


and 


cos 6[f(6) 740 = y(n) dn (13) 


where (12) is the differential equation for the phase dis- 
tribution and (13) is an equation which must be solved 
to find @(n). If the amplitude distribution obtained in 
this paper, (11), is substituted into (13), there results 


cos 6f7() d0 
= 608 [(Om — 90) + 90] f2[(6mn — 80) + Oo] dn. (14) 


Noting that the functional form of both sides of this 
equation are identical, we can immediately conclude 
that 


(7) = (8m — 60) + 8. 


Substituting this result into (12) and integrating gives 


i ; cos [(8m — 00)n + Oo]. (15) 


rr 


Rv 


This result is identical with the phase in (11) and veri- 
fies that the present results are completely compatible 
with the procedure of Dunbar. 

In Dunbar’s paper an example of his synthesis tech- 
nique is presented. For this example the desired pattern 
was | (0) *=csc? cos 6 for the range 10°<@<50°, The 
array length chosen in this case corresponded to k/ = 130. 
This same pattern with identical parameters was syn- 
thesized using the present theory and the exact radia- 
tion pattern was calculated by numerical methods. Both 
Dunbar’s result and the present data are displayed in 
Fig. 6. 
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Fig. 6—Comparison of computed patterns by two synthesis 
methods; kl=130, 00=10°, 0, =50°. 


From these data both methods can be seen to provide 


an adequate approximation to the ideal pattern. De- 


_ pending on the error criterion which is chosen as a meas- 
ure of the accuracy of the approximation, one result may 


fare better than the other. As a result a definitive state- 
ment cannot be madeas to the relative utility of the pro- 
cedures. However, in a particular application one 


~ method may provide a better approximation than the 
- other and may also be more readily realizable with a 


Amal bby 


- practical antenna. 


The method of Woodward is based on the expansion 
of a desired shaped beam pattern f() in a series of the 
following form, 


n=N/2 sin kilsin § — sin On| 
) 


(16) 


@) = 
g(6) TE sin 6 — sin 6, 


where the 0, are taken in increments of pattern beam- 
width between the first nulls. By this procedure the AG 
are directly related to the desired pattern by 


Az = 7 (Gn). 


This procedure obviously involves an approximation 
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which can be shown to be more accurate as the array 
length increases. In the limit of infinite length, 0, be- 
comes a continuous variable w and the summation be- 
comes an integral. That is, 


: ; om sin kl|sin @ — sin 
lim g(@) = lim f(a) ul ; a1 da. (17) 
kl kine J a, sin @ — sina 
But, 
sin A/[sin @ — sin a| 
lim : : = (9 — a) 
k1 sin @ — sina 
so that 
tim g@ =f fle)@—a)da=f0). (18) 
Sa) ag 


Thus, Woodward’s method is exact in the limit of in- 
finite array length. 

Although the above work demonstrates the equiva- 
lence of Woodward’s method and the present technique 
in the limit of infinite length, this in itself is not enough. 
To show complete correlation between the inherent ac- 
curacies of the two methods, it would be necessary to 
demonstrate that both approximations approach this 
limit in the same way, namely 0(1/k/). It has not been 
possible to show this, and it is doubtful that this is true. 
It is entirely possible that Woodward’s method con- 
verges to the limit faster, implying a more accurate ap- 
proximation for finite length arrays. Nevertheless, for 
finite but large arrays (k/~~1000), the difficulties in 
working with Woodward’s technique may not be 
worth the additional accuracy obtained. 


CONCLUSIONS 


A relatively simple and straightforward synthesis 
technique has been demonstrated which may find appli- 
cation in certain antenna problems. The technique 
should find particular application for extremely large 
antennas requiring shaped beam characteristics. In ad- 
dition, the fixed phase nature of the synthesized phase 
distribution suggests applications to variable pattern 
antennas; the possibility of amplitude “scanning” 1s 
an important special case of this class of applications. 

Because of the approximate nature of this method, 
in acritical application an exact calculation of the radia- 
tion pattern should be made to determine whether the 
result is compatible with design specifications. How- 
ever, for arrays with (kl=1000) the results should be 
well within the specifications for any practical shaped 
beam design. 
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APPENDIX 


In the body of this paper the phase distribution which 
is the key to the entire synthesis procedure is written 
without justification. Although the fact that this ex- 
pression is usable is sufficient justification, a heuristic 
derivation is of interest. 

This derivation is based on the following reasoning. 
Let the phase distribution along the array be such that 
the energy from two points at « and «+dx on the array 
add in phase in the far-field direction (@) associated with 
the position x. Mathematically this condition is ex- 
pressed by 


d(x) = — kdx sin [6(x)| (19) 
where ®(«) is the phase distribution and 6(x) explicitly 
denotes that the far field point must be a function of the 
position along the array. This relation can be taken as 
the simplest functional form, namely linear. We then 
obtain 


a(x) = af pees (20) 
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Substituting (20) into (19) results in 


soda = — ksin | es 


dx 
which may be readily integrated to give 


cos | a + 0 | 
Om — Oo L 


= by 


+ 00] 
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This result is identical to that used inm-the radiation | 


integral earlier. 


Implicit in the above formulation is the assumption © 


that the energy radiated from the region dx at x will 


! 
| 
| 
| 
| 


contribute to the far field only in a small vicinity around — 
the direction 6(x). Obviously this is not true since the — 


element pattern of an element dx is essentially omni- 


directional. However, if the array becomes very long, — 


the elements dx can be considered to have finite length. 


The element pattern can then be reasonably directive — 
and will confine the energy mainly in the desired far-_ 
field direction. In the mathematical derivation this ef- _ 
fect occurs by virtue of the increasing accuracy of the | 


stationary phase method with larger kl. 


Gain Limitations of Large Antennas* 


ROBERT C. HANSENT 


Summary—Arrays of many elements and continuous aperture 
antennas are very often designed using Taylor’s “approximate” space 
factor which is optimum in resolution, i.e., it offers minimum beam- 
width for a given sidelobe ratio. This space factor has a number of 
equal level sidelobes followed by tapered lobes, the latter being 
necessary to make the aperture function physically realizable. For 
conventional antennas, these designs offer gain which, while not 
that of a uniform aperture, is quite satisfactory. Very long (or large) 
antennas with their narrow beamwidths, however, may exhibit a 
Severe gain degradation due to large energy contained in the side- 
lobes. Thus the sidelobe behavior must be carefully controlled to 
realize the full gain available. Data are presented for the gain varia- 
tion of various Taylor equal sidelobe designs, and of the modified 
sin x/x one-parameter family as a function of aperture length. 


I. INTRODUCTION 


N recent years, considerable effort has been ex- 
pended towards the design and realization of an- 
tennas optimum in the sense of having the narrow- 
est beamwidth, hence best resolution, for a given side- 
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lobe level. The Dolph-Tchebycheff array distribution 
has been shown to be optimum in this sense [1 ]-[3]. 
To recapitulate, the Dolph design provides that all side- 
lobes be of equal magnitude; this has the effect of nar- 
rowing the main beam compared to that of a pattern 
with tapered sidelobes where the largest is of the same 
amplitude as the Dolph case. The Dolph design is ac- 
complished through forcing a correspondence between a 
Tchebycheff polynomial and the array polynomial. 
These aperture distributions are widely used and con- 
siderable literature exists to aid designers. These in- 
clude the extensive tables [4] of Brown and Scharp 
and approximate formulas [5]. In the limit as the num- 
ber of elements becomes very large, and the element 
spacing very small, the Dolph-Tchebycheff distribution 
becomes a Taylor “ideal” continuous distribution [6]. 
These distributions, which also have the equal sidelobe 
and optimum resolution characteristics, are widely used 
for the design of continuous line sources and apertures, 
and for arrays where the number of elements is large and 
the consequent Dolph-Tchebycheff computations tedi- 
ous. Excellent tables of these aperture functions by 
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Spellmire are also available [7]. The aforementioned 

references are primarily concerned with line source or 
rectangular aperture-type antennas; similar develop- 
ments have been pursued by Taylor for circular arrays 
and circular apertures [8], [9]. For most antennas, the 
; optimum beamwidth design offers a satisfactory gain 
value. However, in recent times, advanced system re- 
‘quirements have necessitated the realization of larger 
_and larger antenna gains with construction of corre- 
_ spondingly larger and larger antennas. Thus, the ques- 
_ tion of optimum gain again needs to be investigated. 
i It has been suggested by H. Wheeler that the opti- 
mum beamwidth antennas do not provide optimum 
% gain, especially where the antenna dimensions are large 
_ in wavelengths.! It is easy to see physically where a 
_ gain limit might occur. Take, for example, a long Dolph 
array or Taylor line source with fixed sidelobe ratio. 
The ‘gain (strictly speaking, the directivity gain) is de- 
_ fined as the quotient of peak intensity (which we will 
keep fixed) to the total radiated power. The latter is the 
5 sum of the main beam power and the sidelobe power. 
_ For modest length antennas, the main beam power pre- 
_ dominates, and as the beamwidth is narrowed (by in- 
creasing the length, of course), the radiated power de- 
.. creases and the gain increases. However, the sidelobe 
power is essentially constant, since the sidelobe level is 
> fixed, and all sidelobes have nearly the same shape. 
When the beamwidth becomes very narrow, the main 
beam power becomes small compared with the sidelobe 
power, and as the latter is nearly constant, so, therefore, 
" is the gain. To remove this gain limitation, one must ob- 
_ viously taper the far-out sidelobe envelope. It is the 
_ purpose of this paper to investigate the effect upon gain 
_ of various tapers for some of the widely used aperture 
3 illumination functions. Similar reasoning shows that 
_ the gain cannot exceed the average sidelobe level. If 60- 
_ db gain is required, the average sidelobe level must be 
- below —60 db. Ruze and others have pointed out that 
% aperture errors affect the sidelobe structure strongly, 


- and this represents another gain limitation. 
a 


a al 


a ee 5. > 


, a Il. OptrimuM GAIN 


% The term “optimum gain” is used here in a very spe- 
~ cial sense. If no restrictions are placed on the aperture 
illumination function, infinite gain is, of course, possible 
% [10], although this fact has not always been appreci- 
~ ated [11]. Bouwkamp and de Bruijn demonstrated that 
arbitrarily large directivity can be obtained for a line 
~ source, and Riblet [12] extended this demonstration to 
a two-dimensional distribution. The arbitrarily large 
,. gain is produced by an interference process in which 
large phase changes and large currents are used to pro- 
~ duce a low value of effective radiating current.” This 
~ well-known phenomenon is called supergain and is gen- 


1 Private correspondence. : ; 

2 For an array with a fixed number of elements, a maximum gain 
does exist. This problem has been solved for a small number of ele- 
ments by Pritchard; see reference [13]. 
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erally avoided because it involves large energy storage 
and thereby high losses, and because of high sensitivity 
to excitation errors [14]. Tucker has pointed out that 
supergain designs are also undesirable because of their 
adverse signal-to-noise performance [15]. Thus, in 
practice, it is necessary to place some limit on the 
amount of supergain to be allowed. This limit can be an 
upper bound on the supergain ratio, which is the ratio 
of total power (reactive and radiated) to radiated power 
[6], or it can be a restriction on the allowable total 
phase change over the aperture, or a restriction upon 
the maximum slope of the aperture function. Determin- 
ing the illumination function which yields maximum 
gain subject to one of these constraints represents an 
exceedingly formidable Euler-Lagrange problem, even 
for the simplest case of a line source. This problem has 
been solved by Proctor.? There are synthesis techniques 
where the gain is maximized subject to a different con- 
straint, such as constant ohmic loss [16]. However, 
these do not appear to relate directly to limiting the 
amount of supergain. A different approach has been 
taken by Solymar [17] and Harrington [18]. Solymar 
determines the maximum gain for a specified supergain 
ratio and number of harmonics of the distribution func- 
tion; Harrington relates the gain to the number of 
spherical harmonics used in expanding the pattern, and 
heuristically relates this to physical size. 

Because of the extreme difficulty of these synthesis 
problems and because of the practical limitations of 
utilizing supergain, almost all aperture antennas are 
designed with a constant phase over the aperture. (For 
beam directions other than broadside, the constant 
phase is of course replaced by the linear phase cor- 
responding to the desired beam direction.) Thus, in this 
paper, only constant phase distributions will be con- 
sidered and “optimum gain” will be used in this context. 
It is readily demonstrated [19] by use of the Schwartz 
inequality that the constant amplitude (usually termed 
uniform amplitude) distribution has the largest gain of 
all constant phase designs.t However, the 13-db side- 
lobe ratio obtainable with the uniform illumination is 
often unsatisfactory for both radar and communica- 
tions purposes. Thus, this paper will consider the varia- 
tion of gain with length for two types of often used dis- 
tributions: 


1) tapered envelope distributions which afford greater 
sidelobe ratios, and 

2) “approximate” Taylor distributions in which the 
far-out sidelobe envelope is tapered, but close-in 
sidelobes are of nearly equal amplitude. 


Il]. UnirormM [LLUMINATION 


Gain factors will be given for the array or line source 
pattern only; the element pattern will be assumed to 


3. K. Proctor, manuscript to be published. 

4G. C. McCormack [20] shows that for an end-fed line source 
with ohmic losses, the uniform distribution must be modified to ob- 
tain optimum gain. 
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be isotropic. Because the gain of a long (of length L) uni- 
form line source is 2L/X, it is convenient to normalize the 
gain, and to consider \G/2L. For a square aperture an- 
tenna (again not including the element factor), the gain 
is found from 


(1) 


2 
Geanire = TG lines 


For a uniform line source, the normalized gain is 
readily written as 


AG 1 
ois — 
U 
1b, )N TU 
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where the normalized variable u=L/\ sin 0, K=27/), 
and Si(x) is the conventional sine integral. 

The gain function is plotted in Fig. 1. For small KL, 
the normalized gain is greater than unity because the 
absolute gain G must equal unity for KZ=0 in ac- 
cordance with the assumption of an isotropic element 
factor. This is just observable in Fig. 1 for KL =10. The 
curve also contains imperceptible oscillations. 
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Fig. 1—Normalized gain vs length for uniform 
and ideal space factors. 


IV. Tayior “IDEAL” SPACE FACTOR 


A continuous source equivalent to the Dolph-Tche- 
bycheff array design is the Taylor “Ideal” aperture dis- 
tribution and space factor [6]. The space factor has 
equal level sidelobes and is given by 

P(u) = cos*(rv/u? — A?), (3) 
where the voltage sidelobe ratio 7 is related to the pa- 
rameter A by 7 =cosh 7A. The aperture distribution cor- 
responding to (3) is given as an infinite sum by Taylor, 
and has delta function singularities at each end. Thus 
the “ideal” distribution can only be approximated in 
any physical situation. 
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The “ideal” space factor gain is given by 


a /2 
G= Puss i f P(u) cos 6d6. 
0 


In order to evaluate the integral, it is separated into two ) 
at the change from cos to cosh, at w#.=A. Call 0; the 


angle related to m, and let 6;=7/2—A, where for long 


antennas A is small so that AX) A/L. So the integral is | 


mw/2—A Ss 
yi+1,= f cos? +/72u2 — m?A? cos 6d 
0 


w/2 
+ cosh? +/72A?2 — 22u? cos 6d0. 


mw /2—A 


The second integral is evaluated by numerical fitting, 
using a sixth-degree polynomial with economized 


(Tchebycheff) coefficients. Error should be <2 per 


cent for sidelobe ratios between 20 and 40 db. 
Do lute MBO Tas 2 An?.482 
n= f cosh? +/72A2 — 4? du = —————_ - 
15 0 Kes 


For J;, approximate the integrand by cos? ru and inte- 


grate from 0 to 7/2. This is equivalent to adding an ex- ~ 


tra fraction of sidelobe in the main beam position. For 


long antennas with many sidelobes and narrow main 


beams, this approximation is good. 
i) sin KL 


X pLir 
fae ahs cos? rudu er cera, (7) 
Thus 
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Fig. 1 depicts these gains for design sidelobe ratios of 


20, 25, 30, and 35 db. The abscissa is L/\, the line | 


source length in wavelengths. As expected, it can be ob- 
served that the normalized gain falls off with increasing 
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(4) 


(5) 


(6) 


length; thus the gain does not increase linearly with | 


length for long lengths and. narrow beamwidths. For 
example, with 25-db sidelobes and L=100 \, gain is 
down to 138 from 178. This would correspond to a line 
source gain of 21.4 db, or a square source gain of 47.8 db. 
Of more interest, however, are the “approximate” Tay- 


lor designs, which do not have singularities in the illu- ~ 


mination function. 


V. Taytor “APPROXIMATE” SPACE FACTOR 


Taylor showed that there is a relation between the far- 
out sidelobe envelope and the behavior of the aperture 
distribution at the line source ends. To eliminate the 
infinite values of end excitation, the space factor must 
have a tapered far-out sidelobe envelope. The “approxi- 
mate” design contains (nearly) equal sidelobes out to a 
transition point NV, beyond which the envelope tapers 
down [6]. Thus the “ideal” pattern is approximated 
arbitrarily closely. Fig. 2 gives three examples of illumi- 
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nation function; as V becomes large, the end peaks tend 
toward a Delta function, giving the “ideal” distribu- 
tion. The space factor is given as a continued product 
over the zeros; only the sidelobe integral J, of (5) is 
markedly different. This is now 


a/2—A N—-1 U2 2 
0 n=1 Ona 


=F (1 — “) cos 6d6, (9) 


n=N nN 


where a,?=A?+(n—}4)? and o is a main beam dilation 
factor close to unity produced by the shifting of zeros 
for M>N to produce the taper (see Fig. 3). For M<N, 
_ the envelope is nearly constant; for 2>>N, it is sin X/X. 
In the region ” near N, both envelopes are modified to 
join together. The integral over the sidelobes is evalu- 
ated by approximating each lobe by a sinusoid. Then 
the integral from 2 to wy is closely 
N-1 
area —— >) (uny1 — Un) 


i 


po eae 
dhe) rien ae 


APERTURE FUNCTION 


DISTANCE FROM CENTER 


Fig. 2—Aperture distributions for 25-db 
sidelobe level for several NV. 


SIN 
x 
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Fig. 3—Sketch of “approximate ” space factor. 
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since 


A? 
Unt = Un = 1 ce ee 
2n? 


(11) 


The error in this approximation has not been evaluated, 
but is thought to be small. For the tapered sidelobes, 
the peaks are at w=n—4 with amplitude 2/m(2n—1). 
In this case, un41—Un=1, and the integral from uy to 
Um 1S 


XN BM 4(uay1 — Un) (M — N)d 


abled, ~~ 5 
yy peer 0) eh Ly cy Gi es bY 8 


(12) 


Now, neglecting the last partial lobe, w»=L/X, so the 


value of J; becomes 
AN 
At i1-— 
IL, 


pee E 1 =| at (13) 
PE oye 12 2n?L(N + 1) 
The expression for normalized gain then becomes 
AG 271? 
phe) 


2L Xr 
2\ 1——— 
mw? A? JE, 
1.939An?-+22 |. N—1——— 4 
12-J- 2 a(N 41) 


Because 2/rN 1.93 7An?, and because 1*/6<1.93n’, 
result (14) can be simplified with small loss in accuracy 
to 


AG 2n? 
TE 93 Aetee BON 21) 


(15) 


If the sidelobe “break point” N remains fixed as L in- 
creases, the normalized gain is constant; that is, the 
gain increases directly with L as desired. A more practi- 
cal design, however, would allow WN to increase with L 
to efficiently use the aperture in producing a narrow 
main beam. Two cases have been computed, for NV 
=M/2 and N=M/4. Each of these cases will exhibit 
a decrease of normalized gain for long lengths of line 
source. Figs. 4 and 5 display this for the two cases and 
for the 20-, 25-, 30-, and 35-db sidelobe levels. 

As an example of the use of these data, suppose a line 
source of length L/N=100 is desired, with 25-db side- 
lobes. For a choice of N= M/2, or half of the sidelobes 
equal level and half tapered, the normalized gain is 0.8 
compared to 0.91 for short line sources. The 0.8 value 
gives a line source gain of 22 db or a square source gain 
of 49 db. For N=M/4, or three-fourths of the lobes 
tapered, the normalized gain is 0.85. 

Eq. (15) indicates that for a fixed length antenna, 
gain is maximized by using a minimum value of N. 
This is incorrect, as (15) does not include the beam- 
width dilation factor [6] ¢ which is 


N 


= ——= * 16 
7 ATH (WN 3) os 
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As WN is decreased, the sidelobe integral decreases, 
thereby increasing gain. However for small N, @ is no 
longer approximately unity, and as WV is decreased to 
this range, the main beam integral increases, thereby 
decreasing gain. Taylor [6] has given the values of NV 
for maximum gain for the limiting case of infinite 
length. Examination of (14) shows that only the last 
term of the denominator contains the length, and if 
N>3,L/\>4, this term will be less than one per cent of 
the preceding term. Thus for any practical length, Tay- 
lor’s results are applicable. Similar results can be ob- 
tained by including o? in the first term of the denom- 
inator of (14), and locating the maximum value of NV. 
Because both main and sidelobe integrals change only 
slightly for small values of NV, the approximation that 
sidelobes are of unit height out to N, whereas these 
actually decrease, produces a shift in maximum calcu- 
lated from (14). Table I gives the values of NV which give 
maximum gain for a fixed length of source. Values of 
normalized gain may be compared with the uniform 
source (L=10 X), where sidelobe ratio is 13.2 db, gain 
1.01, and beamwidth 5.04°. 

For short line sources with a few sidelobes, the ap- 
proximations made in obtaining the gain are sufficient 


10 20 30 40 60 


80 100 200 300 400 600 800 1000 


Fig. 4—Normalized gain vs length for Taylor 
approximate space factor. 
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Fig. 5—Normalized gain vs length for Taylor 
approximate space factor. 
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to indicate a value of gain for the 20-db Taylor case 
higher than uniform; this is, of course, not valid. Good | 


accuracy is obtained for L>20 . 


VI. TAYLOR ONE-PARAMETER SPACE FACTORS 


A family of tapered sidelobe distributions has been ‘| 


derived by Taylor [21]. These are of the modified 
sin X/X type. The space factor (pattern) is given by 


(17) 


The sidelobe ratio 7 is 4.6 sinh 7B/7B; the aperture 
function is 


F(p) = Jo jBV/n? — 9°), (18) 


where Jo(jx) is the Bessel function for imaginary argu- 
ment. 

These distributions offer the designer patterns with 
tapered sidelobe envelopes and with sidelobe ratios 
greater than the 13.2 db obtained for the uniform case. 
Because of the taper, the normalized gain its nearly con- 
stant, and the asymptotic value for large L may be 
used. This is [21] 


AG =s- 2 sinh? 7B 

2L ° wBIo(2eB) 
where I(x) = /oIo(€)dE; see Rothman.’ Values of nor- 
malized gain and beamwidth per wavelength of source 
length are in Table II. These sources have a highly 
tapered illumination function, and may be used where 
a sacrifice of gain and beamwidth can be tolerated in 


order to obtain rapidly decreasing sidelobes, e.g., low- 
noise antenna designs. 


(19) 


TABLE I 
TRANSITION POINTS FOR MAXIMUM GAIN 


Beamwidth 
Sidelobe nee Nex | Maximum | Ber Inverse 
Rati a mr cleng 
atio (Taylor) | Eq. (14) AG/2L @ pertear 
Dilation) 
20 db 6 9 0.97 Soe2 
25 db 12* 18* 0.90 S827 
30 db — 36* ORS O32 14 
35 db —_ 64* 0.67 67 35° 
* Very broad maxima. 
TABLE II 


TAYLOR ONE-PARAMETER SOURCES 


Sidelobe Rati Beamwidth per Inverse 
idelobe Ratio G/2L Wavelength 
20 db 0.933 58 .6° 
25 db 0.863 63 ..9° 
30 db 0.801 68 .8° 
35 db On7st (omer 


5M. Rothman, “Table of Zo (x) for 0(.1) 20(1) 25,” oh 
and Appl. Math., vol. 2, pt. 2, pp. M2 ott: yes ee oa 
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The “ideal” space factor, where all sidelobes are equal 
and the one-parameter family can be considered as spe- 
cial cases of a two-parameter family, developed by Bick- 
more and Chu [22]. This two-parameter space factor is 


e J (aV/u2 — C?) 


F(u) = ey oD ae (20) 


where y is real. C is related to the sidelobe ratio and 
u=L/\ sin 6. The second parameter is y. When y= —4 


2) 
(20) reduces to Taylor’s “ideal” space factor cos 


 (xVw—C), whereas for y=-+4 we get the modified 


one-parameter family sin (4-V/w2— C?)/1V/w—C?. In- 
termediate values of y yield different sidelobe envelope 
tapers. This versatile family of space factors allows in- 
sight into the interrelationship between sidelobe taper, 
beamwidth, and zeros of the aperture function. 


VII. ConcLUSIONS 


The widely used Taylor “approximate” line source 
designs must be carefully used for dimensions very large 
in wavelengths. At the same time, this approximate 
space factor yields the best compromise between gain 
and resolution. It is possible, using the curves given, to 
choose the sidelobe break point so as to realize nearly 
all of the gain available from a stipulated sidelobe ratio 
and yet maintain close to the optimum beamwidth. 
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Some Equivalences Between Equally and 
Unequally Spaced Arrays* 


SHELDON S. SANDLERT 


Summary—General analytical expressions are presented for 
unequally spaced arrays. These relations allow for the analysis of the 
nonuniformly spaced array in terms of its equivalent uniformly 
spaced array. The inherent broad-band qualities of the nonuniformly 
spaced array are discussed. Some equivalence is made between the 
amplitude and spatial variation in the uniformly and nonuniformly 
spaced array. The general synthesis problem is discussed. An array 
with monotonically increasing interelement spacings is presented 
as an example of the theory. 


INTRODUCTION 


N GENERAL, a linear array with uniformly spaced 
] elements produces a radiation pattern which has 

several major maxima of equal height (grating 
lobes). This reproduction of the principal maxima pre- 
vents operation over a broad frequency range. A solu- 
tion to the problem of reducing the minor principal 
maxima has been proposed involving a linear array 
with nonuniformly spaced elements.! This type of radi- 
ator will be referred to as a nonuniformly spaced array. 

A general discussion of nonuniformly spaced arrays 
has been given by Unz.” The classical expansions were 
performed on the far zone pattern. The complexity of 
this formulation limits its practical application. It is 
believed that the present analysis offers some advan- 
tage in understanding problems associated with non- 
uniformly spaced arrays. Furthermore, some important 
properties of nonuniformly spaced arrays may be 
readily demonstrated and should find usefulness in nor- 
mal array applications. For example, there is an equiv- 
alence between the variation of the individual element 
distances from the array center and amplitude tapering 
of a uniformly spaced array. 

As an introduction to the theory of nonuniformly 
spaced arrays consider the uniformly spaced array 
shown in Fig. 1. The far zone field for this array may. be 
found by summing up the contributions due to the indi- 
vidual radiators. With no loss of generality, the analysis 
may be restricted to the case of N isotropic elements 
driven in identical time phase. For this case the radia- 


3 Received by the PGAP, February 5, 1960; revised manuscript 
received, March 25, 1960. This work was supported in part by Rome 
Air Development Center, Rome, N. Y., under Contract No. AF 
30(602)-1776. 

t Electronic Communications, Inc., Timonium, Md. 

1D. D. King, R. F. Packard, and R. K, Thomas, “Broadband 
Steerable Linear Arrays,” Electronic Communications, Inc., Timo- 
nium, Md., prepared for AFCRC, Contract No. AF 19(604)-5234: 
October 7, 1959. 

? H. Unz, “Linear Arrays with Arbitrarily Distributed Elements,” 
Electronics Res. Lab., University of California, Berkeley, Ser. No. 
60, Issue No. 168; November 2, 1956. 


tion field is directly proportional to the uniformly 
spaced array factor given by 


(N—1) /2 


A(®@)=1+2 >> cos 2rmd cos ®, (1) 
m=1 


where 


d=element spacing in fractions of a wavelength, 
N=total number of elements. 


It is easily seen from (1) that the array factor is 
periodic since it repeats when the argument changes by 
an angle of m radians. The usual method of reducing 
(t.e., to zero) the level of the secondary maxima is to 
restrict the spacing. However, the secondary maxima 
may not be reduced when the element spacing exceeds a 
half wavelength. For example, an array with two-wave- 
length spacing will have its second principal maxima at 


&=7/3 repeated at 6=0. 


1 
| 

| 

| 

I 7 

$14 
Fig. 1—Uniform array. 


NONUNIFORMLY SPACED ARRAYS 


A useful physical picture of the nonuniformly spaced 
array is concerned with the concept of spatial frequency. 
Spatial frequency is a measure of the variation of each 
element contribution with far field angle. With this 
definition, the lowest array frequency may be associated 
with the center element, the next highest with the m=1 
term in (1), etc. The concept of spatial frequency may 
be visualized by considering two sources separated by 
a distance d. The radiation pattern for the two sources 
will have an increasing number of visible lobes as d is 
increased. The linear array consists of a large number of 
such couplets since the center element serves as a ref- 
erence. With respect to the center element, the uniformly 
spaced array is seen to contain a spectrum of frequencies 
which are integral harmonics of the fundamental fre- 
quency determined by the first element position. In this 
case, each spectral component is located at the physical 
element position, and the locations are uniformly dis- 
tributed over the length of the array. 


4 
q 


a 
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Consider the nonuniformly spaced array factor given 


by 


A(%) = 1+ 2 >> cos (21dm cos ®), (2) 


m=1 
where 


dm= distance from the mth element to the center of 
the array in fractions of a wavelength. 


The nonuniformly spaced array is characterized by 
spatial frequencies which are not necessarily related by 
integers. The spatial frequency spectrum is now un- 
evenly distributed over a range which extends from zero 
to wy, where 


wu = 2mdy cos ®. (3) 


The lowest nonzero array frequency plays a major 
part in determining the location of the principal and 
minor maxima. Generally speaking, when a function is 
composed of slowly as well as rapidly varying functions, 
the slower variations will predominate. This is true 
since it is always possible to find some high-frequency 
maxima in the vicinity of a low-frequency maximum. 
The total number of high-frequency maxima in the vi- 
cinity of the low-frequency maximum then determines its 
magnitude. The major maxima differ from the minor 
maxima in that they represent simultaneous maxima of 
both the low- and high-frequency components of (2). 

Unfortunately, compact expressions are not available 
when the argument of (2) isa nonlinear function of m. A 
perturbation analysis could be performed to give an in- 
dication of small nonlinear changes in element spacing. 
However, some of the useful properties on nonuni- 
formly spaced arrays depend on large nonlinearities in 
the element spacings. A method is needed which will 
lend insight into the gross behavior of the nonuni- 
formly spaced array and make it susceptible to a linear 
analysis. Fortunately it is possible to make a cor- 
respondence between the nonuniformly spaced and 
equivalent uniformly spaced array (EUA). The EUA 
is the best mean square representation for the original 
array. Some superficiality is introduced since the EUA 
does not physically exist. Furthermore, it will be shown 
that the spacing in the EUA is somewhat arbitrary. 
The EUA involves the transformation of higher fre- 
quency terms in (2) toa set of integral harmonic terms. 
The lowest nonzero superficial array frequency in the 
EUA is a, where 


w, = rR cos ®. (4) 


The scale factor R determines the element spacing in 
the EUA. For example, the value K=1 corresponds to 
an EUA with half-wavelength spacing and R= cor- 
responds to quarter wave spacing. The transformation 
to the EUA is made by finding the Fourier representa- 
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tion for each higher harmonic.’ Each term in (2) is re- 
duced to a series of integral harmonic terms with vari- 
able amplitude coefficients. For example, consider the 
general term in (2) given by 


COs wu, (5) 
where 


Wm = 2mrdm cos ®&. 


Since m>1, w@m>w, (5) may be expanded in integral 
harmonics of the scaled variable w; where 


w, = tR cos ® 
O <a RoS (6) 


Once the lowest EUA spacing is chosen, each term in 
(2) must be expanded into an infinite number of uni- 
formly spaced equivalent elements. The amplitude of 
the equivalent elements varies considerably across the 
array. Hence, for practical computations only a few 
terms of the expansion need be considered. Each term in 
(2) will have an integral harmonic cosine representation, 
and the individual expansions may be added term by 
term. When the element spacing of the EUA is de- 
creased toward zero, only those EUA elements which are 
very close to the original will have an appreciable mag- 
nitude. In the limit as R approaches zero, only one 
term will appear in the EUA representation of the 
original. This will also be true of the other EUA repre- 
sentations for the original elements. The sum of all the 
uniformly spaced representations of each term will then 
be identical to (2). 

Note that the EUA has all the general characteristics 
of the normal uniformly spaced array (e.g., grating 
lobes). Different values of the scale factor RK may be 
chosen for specific applications. From the standpoint of 
reproducing the actual radiation pattern, an EUA with 
a small value of R would give the best fit. Large values 
of R would yield an EUA which represents the gross pat- 
tern characteristics. The fine detail of the radiation pat- 
tern is filled in by a large number of EUA elements 
(i.e., small values of R). 

The first step in the EUA representation is to expand 
the mth term in (2) in an integral harmonic cosine ser- 
ies. When (6) is substituted in (5), the result is given by 


Lins 
COS Wm = COS 75 W1 = COS MmW1, (7) 


where 

Ben, = 20m IR 
Now the mth term may be expanded into v integral har- 
monic terms by the following expansion: 


3R, Courant, “Differential and Integral Calculus,” Interscience 
Publishers, Inc., New York, N. Y., vol. 1, p. 444- 1955. 
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(oe) 


COS Minit = >, Any COS Yuri, (8) 


v=0 


” TR 
Any = =f COS Mmw1 COS vw1du1, 
T 0 


2 E SIN MmW1 COS PW1 + VY COS MmW1 SIN =)" 


a — i 
Con = ) IR = 1 
sd D) 2 
T Mm” — V 


Paniee Op se dae 2 tt stare 


The Fourier expansion given by (8) must be applied to 
each component in (2). The amplitude coefficients are 
added term by term to give the EUA, or 


iP c-} 
A=14+'2 >} Do GrniCOs rast 


m=1 v=0 
iP 

A=1+2)>0 A, cos roy (9) 
r=0 


where 


Ao = @19 + Goo + G30 + °°: 
Ay = G11 =F Gor + Gi + *- ~, etc. 


Note that when wm is an integer, (8) may be replaced by 
a trivial identity. The representation given by (9) is 
not unique. Many choices are possible as a function of 
the scale factor R. However, from the least mean square 
error property of the Fourier expansion, each repre- 
sentation will be the best possible for a given a. 

The final form of the array factor in (9) is readily re- 
ducible to another analytical form. The reduction of 
(9) is performed by representing the values of A, at P 
points by a Pth degree polynomial in m.4 The array 
factor given in (14) is now given by a finite series of 
form 


P 
A=B+2)>> Am cos ma, 


m=1 


P P 
=B+ ys eS a;m" COS Mw, 


m—=1 r=0 


P ie 
A =B+24 Y ancosmos + > am cos mw; + toe 


m=1 m=1 


P 
+ >) amm? cos mosh se (Cis), 


m=1 


The individual terms in (10) may be summed exactly 
as sums of successive derivatives of the functions S(w) 


4F. B. Hildebrand, “Introduction to Numerical Analysis,” 
McGraw-Hill Book Co., Inc., New York, N. Y., pp. 261-269; 1956. 
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and C(w;),> where 


A ee 
sin — w; sin ——— 
M. 2 
S(w1) = Do sinma: = pont ds) 
m=1 ered 
sin — 
ee 
MM M+1 
sin —w; Cos — oy 
M 2 
C1) = SS COS Mw, = (12) 
m=1 ne) AOD 
sin — 
For example, the third term in (10) is given by 
M 
>> am cos ma, = a1 — S(w1). (13) 
m=1 dw 


A simple approximate method for making an analysis 
of unequally spaced arrays may be based on the EUA 
and the work of Spencer. The EUA may be approxi- 
mated by an equivalent line source distribution by 
interpolating with continuous functions. These func- 
tions may be obtained empirically or by least squares 
fit. Note that the original form of the array factor given 
by (2) does not lend itself to this type of representation. 


EXAMPLES 


Amplitude Tapering With Equally Fed Elements 


The sidelobe level and beamwidth of a uniformly 
spaced linear array may be controlled by varying the 
element driving currents. This amplitude variation cor- 
responds to the variation of the 4,’s in (9). A general 
result of uniformly spaced array theory is that the side- 
lobe level may be reduced by decreasing the current 
amplitudes toward the extremes of the array. This re- 
duction is accomplished at the expense of the array 
beamwidth. Now it will be demonstrated that there is 
an analogy between amplitude tapering of a uniformly 
spaced array and the space variation in a nonuniformly 
spaced array.’ 

As a specific example, consider the array of Fig. 2. 


a =/,0000A 


d, =2,0053A o=> 
ds =3./660A a 
dy 74, 5388A > 
=o O-—6—9-8-6-—6—o oe 
4’ } 2 te } 3 3 4 


Fig. 2—Logarithmically spaced array. 


5 F. B. Hildebrand, “Methods of Applied Mathematics,” Prentice- 
Hall, Inc., New York, N. Y., p. 260; 1952. 

®R. C. Spencer, “Fourier Integral Methods of Pattern Analysis,” 
Mass. Inst. Tech., Cambridge, R.L. Rept. No. 762-1; January 21, 
1946. 

7 This result was experimentally demonstrated” during” World 
War II at the Navy Underwater Sound Laboratory, Harvard Uni- 
versity, Cambridge, Mass. (Prof. F. V. Hunt, Harvard? University 
private communication). a i 
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The elements in this array are arranged so that the 
d;’s increase logarithmically from the center of the ar- 
ray. The value of w; for this discussion is chosen to be 
proportional to half the electrical distance to the first 
element, or 


w1 = 7 Cos ®, R= 1. 


The EUA is readily constructed by applying (8) to 
each high-frequency component. The result is shown 
graphically in Fig. 3. The EUA may be decomposed into 
two arrays with one-wavelength spacing. The presence 
of amplitude tapering may be seen from the amplitude 
variation of the subarray with an odd number of ele- 
ments. Since this array yields the major contribution to 
the radiation pattern, its effect will be more pronounced 
than the “even” array. The radiation pattern shown in 
Fig. 4 demonstrates this reduction. Note that the radi- 
ation pattern is plotted as a function of the broadband- 
ing factor U=R cos ®, where R=2. The utility of U 
may be seen from the following example: consider a uni- 
formly spaced array operating at frequency f; with an 
element spacing of one wavelength. For this case R=1, 
and the visible range of U is —1<U<1. Now if the fre- 
quency is doubled to 2f,, then R=2, and the visible 
range of Uis —2<U<2. Returning to Fig. 4, note that 
the sidelobe level is less than 16.8 db over a +45 degree 
range, compared with 12.8 db for a uniformly spaced 
array of nine elements with one-wavelength spacing. 

The radiation pattern shown in Fig. 4 was calculated 
on the basis of the original element positions. The 
truthfulness of the EUA representation of Fig. 3 is 
shown in Table I. 

The error in the EUA representation is greatest at 
the rapidly varying portions of the radiation pattern. 
An increase in the number of EUA elements would 
improve the approximation to the actual radiation pat- 
~ tern shown in Fig. 4. 


BROADBAND PROPERTY OF NONUNIFORMLY 
SPACED ARRAYS 


A useful property of the nonuniformly spaced array 
is the reduction in the magnitude of the major prin- 
cipal maxima. The broad-band property of nonuni- 
formly spaced array may be crudely explained by the 
division of the EUA into an even and odd array. It is 
well known that an array with an odd number of ele- 
ments has principal maxima of the same sign, while an 
array with an even number of elements has alternate 
signs. Returning to the example of Fig. 3, the array has 
a comparatively large “odd” array and a smaller “even” 
array concentrated near the extremes. Due to the alter- 
nation in sign of the second principal maxima for the 
two types of subarrays, the even array increases the 
major maximum of the odd subarray while decreasing 
the height of the second principal maximum. Since the 
“even” array is located near the end elements, there will 
be an increase in the over-all sidelobe level. 
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The EUA for a Small Value of Element Spacing 


In the limit as R becomes vanishingly small, the EUA 
and the original nonuniformly spaced array become 
equivalent. This result is demonstrated in Fig. 5. A non- 
uniformly spaced array of length 38.6 \ has been ap- 
proximated by an EUA with 0.222 Xd spacing. The geo- 
metrical equivalence between the dominant EUA ele- 
ments and the original nonuniformly spaced array is 
clearly demonstrated. 


ENVELOPE FOR SUB~ ARRAY WITH 
ODD NUMBER OF ELEMENTS 


SUB—ARRAY 
WITH EVEN 


NUMBER OF 
7 ELEMENTS 


HALF WAVELENGTHS 


RELATIVE MAGNITUDE OF DRIVING CURRENTS 
S 
ES 


X- ACTUAL ELEMENT LOCATIONS 
la 


d, «1, 00004" 
IN HALF ARRA 


dy =2,0053 A 
ds =3.1660 A 
dy =4.5388 A 


Fig. 3—Equivalent uniform array for array shown in Fig. 2. 
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Fig. 4—Radiation pattern for array shown in Fig. 2. 
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0.100 —25.1 db —23.5 db 
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DISTANCES (ORIGINAL) 
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dg =/5.8 A 
Gq 2/9-3 A 
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d,= LOA 
de = 214 
ds = 354A 
dg =4.92 
d5=6.64 


Fig. 5—Equivalent uniform array with 0.222 \ spacing. 
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SynTHESIS oF NONUNIFORMLY SPACED ARRAYS 


The complicated nature of the A,’s in (9) does not 
allow a one to one correspondence between the actual 
element positions and the EUA driving currents. For 
example, if R=1, the following set of equations must be 
solved simultaneously for the element positions, given 
the value of A,: 


sin wir = SiN Wow = SIN gr 
= eer aera ’ 
Mit _ bet Mest 
y 24 SIN wir 2u2 SiN wor 2u3 SID ws 
‘ltr ay nm ae Le 
rut — 1) eu? — 1) et — 19 
y 2u1 SIN wir 2us SIN Mor 2u3 SiN 37 
2 RP SSS a: je ie 
(are 27) uy =" 27) rug” 2?) 
ple eWwere lo 1 2S) cue: 


Note that the number of A,’s that may be independ- 
ently chosen is less than the number of original ele- 
ments. Hence, if the EUA has more elements than the 
original array, only N uniform driving currents may be 
specified, where 2N-+1 is the total number of elements. 
If the mode of subdivision is changed to an EUA with 
a larger number of elements the same limitation holds 
true. A reasonable approach to the synthesis problem 
would involve an EUA having a total number of ele- 
ments commensurate with the original nonuniformly 
spaced array. 
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In practice, it is rather easy to see the effect of moving 
an array element graphically on a representation of the 
type shown in Fig. 3. When the element frequency is 
close to a harmonic of the fundamental frequency, that 
harmonic will be large. When the element frequency is 
between two array harmonics, each adjacent harmonic 
will be large, though not as large as in the preceding 
case. 

If the array elements are widely spaced, it is difficult 
to find the optimum distribution for supressed minor 
principal maxima. For example, consider the operation 
of a nonuniformly spaced array over a two to one fre- 
quency band. If the elements in the EUA were spaced 
one wavelength apart at the low-frequency range, they 
will be spaced two wavelengths when the frequency is 
doubled. It follows that it is difficult to make any rigor- 
ous mathematical statements about optimum non- 
uniformly spaced arrays. However, the EUA analysis 
may serve as a guide in specific applications. For ex- 
ample, consider the design of a nonuniformly spaced 
array with a minimum of one-wavelength spacing. The 
array must operate over a two to one band with no 
visible grating lobes. The optimum nonuniformly 
spaced array would be given by an element distribution 
which had an EUA identical to a normal uniformly 
spaced array with less than half-wavelength spacing 
Ce. a he 
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The Efficiency of Launching Surface Waves on a 
Reactive Half Plane by an Arbitrary Antenna* 


JULIUS KANE} 


,Summary—The following two-dimensional problem is considered: 
An arbitrary antenna is located at the origin which launches surface 
waves on a reactive half plane y=0, x<0. It is assumed that the 
other half plane, y=0, x>0, is occupied by a perfect conductor. 
The antenna is approximated by a distribution of multipoles easily 
found from the antenna’s radiation pattern. With this approximation, 
simple algebraic expressions are obtained for the ratio of the power 
delivered to the surface wave as compared to the power in the radia- 
ated cylindrical field. 
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List oF SYMBOLS 


6(x) =dirac delta function, 
Ey(v) = pw(cosv/k) =an entire function, a polynomial, 
HH) (kp) = Hankel function of the first kind 
of order zero (a time factor 
e~*** is suppressed), 
k= | | e*®=propagation constant with a 
small imaginary part propor- 
tional to 6, 
K=a constant proportional to the 
source strength, 
/k? —y* =the radical is chosen by the pre- 
scription »/2= +| | BED. 
A=)/k, 
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m= JR, 
M=Vi+hi, 
v=k cos ¢ Fourier transform vari- 
ables, 
p.o.(3) =radiation pattern of an arbitrary 
antenna, 
bw) =truncated Fourier expansion of 
p.o(@), 


P(!) = power radiation pattern, 
II= Poynting flux, 
x = 
‘ a(v) = ——— = af ) 
1v/k2+ py? o(v) 


= Wiener-Hopf decomposition 


U(#) = Heaviside unit function. 


I. INTRODUCTION 


4 T IS well known that a reactive surface which sup- 
ports surface waves can be represented mathe- 
matically as 


—-+ hu = 0, 
on 


where 0u/dn is the normal derivative and ) is propor- 
tional to the surface reactance [1]. In this paper we con- 
sider rigorously the excitation of surface waves on such 
a half plane by both a line source and an approximation 
to an arbitrary antenna. In Sections II-IV, we solve 
the problem of excitation of a surface wave by a line 
source. Using Kay’s results [1], it is shown that the 
ratio of power delivered to the surface wave as com- 
pared to the radiation field is an elementary algebraic 
function of X. 

In Section V, we study the surface wave excitation 
produced by an arbitrary antenna. The antenna is ap- 
proximated by a distribution of multipoles obtained 
from a knowledge of the antenna’s radiation pattern 
p.(0) when placed on an infinite metallic plane. The exci- 
tation produced by these multipoles is then calculated. 
It is shown that the ratio of energy delivered to the sur- 
face wave as compared to the cylindrical field is an ele- 
mentary algebraic function of p..(8) and X. 

We omit several details of analysis which are covered 
thoroughly by Kay [1] and Bazer and Karp [2]. 


Il. FORMULATION AND SOLUTION 


The problem (cf. Fig. 1) we shall solve is posed by the 
three equations: 


(V2 + B)u(x, y) = — 28(x)6(y), alla,y20 (1) 
PF 
See eet pe Oy =10 (2) 
oy 
F) 
“= 0, n>0,y=0 (3) 
oy 


where u(x, y) represents the transverse magnetic com- 
ponent H., and H,=H,=0. For convenience, K has a 
small positiveimaginary part k= | k| e®, and a time factor 


Kane: Launching Surface Waves on a Reactive Half Plane 


501 


Veu+k*u = -2 8m Sy) 


Fig. 1. 


e-*t is surpressed. Owing to the normalization! chosen 
for the source term, the free space Green’s function of 
(1) is 


eietiVe—y*ly| 


i 1 
Halk -—f aT 4 
Re Ss) eres . 


where the contour C is taken along the real axis in the 
v-plane as shown in Fig. 2. Since (4) already has the cor- 
rect singularity at the origin, it follows that if we seek 
a solution of the form 


1 eietiVe yy 
ula) => [ aX) se) 
where 
Ax) = 0) “as | >| > 0, Imv20, (6) 


then (1) will continue to be satisfied. At the boundary 
we assume that we can differentiate freely under 
the integral sign, and then if At+(v) is an analytic func- 
tion in the upper v-plane, the boundary condition (3) 
will be satisfied, since we can close the contour C for 
x>0. Again, it is sufficient that if G-(v) defined by 


co) = Are) 14 (7) 


wESa 


is an analytic function in the lower v-plane and O(1) as 
|y| >, Im »y<0 the boundary condition (2) will be 
satisfied. Kay [1], and Bazer and Karp [2] show that 
the function a(v) can be expressed as 


ot (y) 


o-(v) 


(8) 


o) =[14+ ae - 


1 Recall that the usual source strength of —475(x)5(y) leads to the 
Green’s function iH“ (Rp). 
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Imv 


@ —-argvk* v7 =5 + 3) 
O-—-arg/k*v? =-248 


Fig. 2—Complex v plane. The magnitude of 6 is 
greatly exaggerated for purposes of clarity. 


where ot+(v) is analytic in the upper v-plane Im »y2=Im 
— kand o~(v) isanalytic for Inv $Im—R. It follows that 
(7) can be written as 


G-(v)o-(v) = AtO)ot(r). (9) 


It can be shown? that |o+(v)| and |o-()| are bounded 
above and below by positive constants in their respec- 
tive domains of analyticity. Owing to the postulated be- 
havior of G-(v) and At(v), the left side of (9) is bounded 
and analytic for Im vy S0 and the right side is bounded 
and analytic for Im y=0, implying that each side is the 
analytic continuation of the other. If we invoke the 
Liouville theorem, we discover that each side of (9) is 
equal to some constant K proportional to the strength 
of the source, whence 


G(v)o-(v) = At()ot) = K 


or 
ACA) = 
(v) eh eR 
ore 
a (v) 


Observe that (10) yields kernels A+(v) and G-(y) having 
the required asymptotic behavior for large |v]. 
The integral representation of the solution (5) is then 


—f 1 
2rid ¢ atv) 


or alternately by (8) 


eirtiv ky? y 


a/ k? = p2 


u(x, y) = dy (11) 


eivebiv k?—y* y 


K 1 
re ie a (v)a(v) Wk? — yp? 


u(x, y) = ie GES 


* For example, see Bazer and Karp [2]. 


IRE TRANSACTIONS ON ANTENNAS AND PROPAGATION 


September 


From (11) we observe that the poles of the integrand 


are the zeros of ot(v). Since |ot(v)| is bounded above © 


and below in the upper half plane, these can only arise 
for Imyv<Im &. For reference we locate the poles of the 
integrand of (11). We calculate that 


1 f 1 ik? — v? 
oO) oH OheGY cue eee 
1 in/ kh? — v?(X — in/k? =v?) 
on) (VFB VFB + 9) 


The only pole that can arise in the lower v-plane is at v9 
given by 


vo = — Vk? +)? (13) 
The corresponding residue of (11) is found to be 
= e iV ey (14) 
o(—VEP EM) DOL e 
which represents the surface wave contribution. 
Il]. Tae Far FIELD 
Introduce the transformations 
y= kcos@ (15) 
and 
x = pcos, y = psin’d (16) 
so that (11) becomes 
(6,8) = -—f em enay (10) 
u(p, alg Se etke cos (o— 
: 2mid p ot(k cos ¢) - 


in which the contour C has been transformed into the 
path P (Fig. 3). If we transform the alternate represen- 
tation (12) to the g-plane, we obtain 


1k sin vg ete cos (oS) dp 


(0, 9) ~ f 
U\p, = 
: Oni Jp Ooh CoseiGR oe Cent 


(18) 


from which it is seen that the pole at vy) which generates 
the surface wave now arises for go given by 
rN 


sing = —— or 
tk 


Cosgo = — VW1+ (A/k)2. (19) 
An asymptotic development of (17) or its equivalent 
(18) is readily obtained by deforming the path P into 
the steepest descents contour S(?). For a given & this 
contour is defined as the locus of points in the g-plane 
which satisfy the equation 


cos (g — #) = 1 + is?, (20) 


with s as a real parameter. It is easily established that 
S@) originates at Re g=d3—7/2, Im g=~, passes 
through the saddle point g=0 on the real axis, and 
terminates at infinity Re g=3+7/2, Im p=—~, In- 
spection of Fig. 3 reveals that for O<$0<7/2 we can 
deform C freely to S(@) without picking up the residue 


1960 
iS Steepest descents 
“ contour S(@) for 
= $ 5 BUN Sear 
Bla 
ss ¢ = cos) 
~ 
re 
sue 
ce) = 2 Oia. Re ¢ 
lou 
Z 7 
S(@) for a 
o<@< - ~~ 
e-m Ra 
+P 


Fig. 3—Complex ¢ plane. 


of go, and for some # in the interval w/2<0S7 we pick 


up the residue. With the aid of (19) and (20) it is estab- 


~~ & 


lished by comparing the real parts of (20) that this crit- 
ical angle 3 is given by 


= —1 Cie ) (AD) 
‘pee Cesar: 


where 1/250, 7, so that we only obtain the surface 
wave on the left for x<0. 
For kp>>1, we obtain the far-field behavior 


en ima—hy 
ute). = eee U(o — 9.) 
een a (—m)m 


ei(ket+n]/4) 


—_ Baye (22 
4 at(k | uy ol( P) | 2) 


where 
m= S/R + (23) 


and U(#) is the Heaviside unit function. The result (22) 
should be modified somewhat if \~0 to allow for a pole 
in the vicinity of the saddle point for 0 =m. This detail 
is fully discussed in [2] and [4] where the procedure for 
isolating the pole is explained. ; 
Apart from terms o[(kp)*], we see that the radia- 
tion field contains a cylindrical wave p(d)e*/+/kp with 
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a pattern function p(@) proportional’ to 
1 cae 

ee COS) 24 

at(k cos?) 3 ) me 


and in addition a surface wave field whose amplitude is 
proportional to 


1 


em) 


(25) 


For each field the power carried off to infinity varies as 
the absolute square of these coefficients. As Kay and 
Lurye have observed, it is an extraordinary occurrence 
in that the absolute squares of (24) and (25) are simple 
elementary functions. From Kay’s paper (66), we obtain 
in our notation: 


we 
| o-(—k cos#) |? = se en PO (26) 
cost 4/1 + (\/%)? 
and 
| o-(—m) |? = | o(-V# #») |? 


VPP MER mk 
Can 


Note that our \ corresponds to Kay’s a. 


(27) 


IV. THE ENERGY DISTRIBUTION 


In this section we shall obtain an explicit expression 
for the ratio of the power delivered to the surface wave 
field as compared to the cylindrical radiation field. From 
(22) we see that the surface wave is 


\K Vat 
Ua aa eH iVRM cr (F — 3,). (28) 
surf mo (—m) 
The Poynting flux II transported by this field is readily 
calculated if we recall that u(x, y) represents the trans- 
verse magnetic component H,. 


2 1 Tou 

Wert ff |B | mice dy ae —f = u* dy 
0 twed 9 LOX Sure 
nN 


K 2 


 Umwe a (—m) 


If we use Kay’s results, our (26) and (27), then 
A| K|? 
we(s/1 + A? +1) 


Isurt = =tanh(4sinh A) (29) 


where A=A/k. 


In the same’ fashion we can compute the radiated 
power of the cylindrical field by using (12), 


K 
ot(k cost) \/2rkp 


Weyl = eithotal + + Ol (kp)—*/?]; (30) 


8 Since o(y) is even in », we have o*(v) =1/o-(—»), [2]. 
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whence 


1D basen =), ise dd.* oy pad: 
0 


1 T 
az —f [upue*]ey1 pad 
UWE 0 
ele dd 
Qawe J | ot(k cosy) |? 
|K|? ¢* cosd + 1 
-——f See Bea) 
Qrwe Jy cosd + V1 4+ (A/k)? 


The integral that appears in this expression is eval- 
uated in Appendix I; we find 


| K|? oe 
cyl — 2, pene a KS 31 
“i sa 2+ aa ES 6) 
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the excitation produced by such a radiator? If we re- | 


place the surface wave structure by a metallic plane, 
then the antenna @ will produce some radiation pattern 
p.(8) which will characterize the antenna. The pattern 
p.(8) is known for 0 $0 S7, so that it is a simple matter 
(Appendix II) to approximate it arbitrarily by a cosine 
polynomial 


N 


Pol?) ~ pr) = DI an(k cos 3)", 


n=0 


(36) 


since the cosine functions cos nd or cos" J are complete 


on the interval 0 $0 <7. In particular, the sine functions _ 


sin nv are expressible as infinite series of cosine functions 
for 0S0 S77. We observe then that 


eivatiV ke, —y* 2y 


The ratio Hgsur:/Tey: represents the efficiency of launching He > Opal wae ea (37) 
the surface wave. With the aid of (4) and (31) we obtain ni @ nae Lp k? — y? 
I sure A | A(A = /1 a= ik) | (32) 
Me (+ VIF A)L2A(A— Vi+ A) +A-VIF A? +1 


For A>>1, the bracketed expression behaves like 2A so 
that 


eee 


(33) 


lim SOE 


Nagiee Wey 


Consequently, even without using a directive source, the 
major portion of the radiated energy is apportioned to the 
surface wave as we increase the surface reactance A. 

It is perhaps of greater use to compare the surface 
wave energy to the total radiated energy 


Tey ae ll erste 


Tf we use (4) and (31) a little manipulation reveals 
ay A , AM yt 
we LM +1 


2A(A — M) 
where M=+/1+ A? so that 


TWietal — 


Tltotal = 


+1], 69 


is a function which possesses the approximate far-field 
behavior (36). It is also clear that (37) satisfies the wave 
equation (cf. Fig. 5) 


N Gis, 
(V2 + u(x, y) = —2 Dy =, FP(@)80), 20, (38) 


where 6) (x) =0"/0x" 6(x), and the boundary condition 


Seles (39) 
oy ¥ 


at y=0. 

Eq. (38), reveals that what we have done is to replace the 
antenna by an approximate equivalent—a distribution 
of multipoles. We shall assume that the a, are real, 
which corresponds to feeding the multipoles in phase in- 


Toure A | 
Tltotat 1 + V1 + A? 


which is plotted in Fig. 4 (opposite) as a function A. 


V. ARBITRARY EXCITATION 


In actual practice one would not excite a surface wave 
structure by a line source. We would use a directive an- 
tenna @ such as an exponential horn to focus the max- 
imum energy into the surface wave. How can we obtain 


A aan tt] 
M+1 


2A(A — M) k32) 


tervals of 7/2. One need only solve (38) subject to the 
two-part boundary condition 


= y=0,x%>0 (40a) 


yf = 10 nae 0; (40b) 
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Fig. 5—Approximation to an arbitrary surface wave launcher. 
A Tsurt/Ttot Tey1/Ttot Tey1/Usurt 4 
0.00 0.0000 anno as to obtain the effect of introducing the surface wave 
0.02 0.0198 0.9802 49 5050 structure in the left half plane. 
ae ape oo te Our procedure is to choose 
0.08 0.0768 0.9232 12.0200 1 pivetiVr—vty 
0.10 0.0950 0.9050 9.5249 wl Y) =< | A*(y) En() yee dv, (41) 
0.12 0.1128 0.8872 7.8632 De bis “ 
0.14 0.1302 0.8698 .677 : ; : , 
0.16 0.1473 0.8527 aa where Ey(v) is a known entire function of v, 4.e., a 
0.18 0.1639 0.8361 5.1002 Tchebycheff polynomial approximation defined by 
0.20 0.1802 0.8198 4.5495 N 
O22 0.1961 0.8039 4.0998 = 
0.24 0.2116 0.7884 3.7258 En(7) Xu ro (42) 
0.26 0.2268 ORTIS2 3.4101 ee: 
0.28 0.2415 0.7585 3.1401 En(k cos 8) = py(8) 
nee aeehe oe sas If A+(v) is an analytic function of » for Im v0, the 
es ee pee ee boundary condition (40a) will be satisfied. In addition, if 
0.38 0.3102 0.6898 2.2234 Ato) = 00) inl, (43) 
ave yore ey: ee: then (41) will have the singularity at the origin de- 
0.44 0.3475 0.6525 1.8779 manded by (19). For «<0, y=0 we observe that if 
0.46 0.3593 0.6407 * 1.7834 
0.48 0.3708 0.6292 1.6971 x 
G-(v) = AtW)E i ———— 44 
0.50 0.3820 0.6180 1.6180 @) @) | Ae 5| Ce 
0.60 0.4338 0.5662 1.3052 
nate ate Lote Aree is an analytic function of algebraic order in the lower 
0.90 0.5546 0.4454 0.8030 y-plane, (40b) will be satisfied. If we assume (43), then 
1.00 0.5858 0.4142 0.7071 we see by (44) that 
1..25 0.6492 0.3508 0.5403 vs es e 
1.50 0.6972 0.3028 0.4343 G-(v) = 00"), |r| a. (45) 
175 0.7344 0.2656 0.3616 : 
Use of (8) yields by (44) 
2.00 0.7639 0.2361 0.3090 
2.25 0.7878 0.2122 0.2694 G-(v)o-(v) = Atv) E()o* 0). (46) 
2.50 0.8074 0.1926 0.2385 
2.75 075238 0.1762 aah) Once again we repeat the argument that each side of 
3.00 0.8377 0.1623 0.1937 (46) is analytic in a half plane with a common region 
on 4 eur hee sees (at least the real axis) of analyticity and hence is the 
6.00 0.9172 0.0828 0.0902 analytic continuation of the other. Owing to the as- 
8.00 ee aes see sumed growths, we obtain by an appeal to the Liouville 
Se : theorem that each side is equal toa polynomial ®%(y) of 
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order NV, hence 

At(v) EQ@)ot(v) = OF (rv) 
or 
PN(v) 
EQ) 


Since Ey (v) is a polynomial with real coefficients, its zeros 
are either on the real axis or conjugate pairs. Unless 


ANY) a*(0) = (47) 


ON ANTENNAS AND PROPAGATION 


From (21) we have 
E(—m) = py (cos *) = py (cos! V/1 + A?) 


so that the amplitude of the surface wave is propor- 


tional to the reactive power delivered by the antenna to — 
the complex angle } =cos! +/1+ A. wes the cyl- — 


indrical radiated power becomes 


@(v) =K En(v), where K is some constant, the zeros of ‘ ; 
Ey(v) in the right side of (47) will contradict the ana- en er yes | K| — f" pv) 
lyticity of Atot for Im v2=0, whence a aol 2rwe Jo |ot(k cos’) 
Ap) AQ OT 
Oa (48) . 
Kote at cos? + 1 
Iey1 = i ite di, 55) 
3 2 i 2mwe J 9 | pul) igs osd + fi + A? ( 
K E(v) eivativir—o™ yy a 49) 
UA) ie = ot(v) Vk — v® a ( and then 
Isur Qa a cos? + 1 
“= = | E(-m) bf | (8) |? dg (56) 
IIcy1 (Pa ay igi) ae 0 cosd + +/1+ A? 


and again the constant K is proportional to the strength 
of the source. A straightforward application of the 
method of saddle points yields 


lim u(p, 8) = wart U8 — 8.) + weyi + O[(Rp)-8/2], (50) 
p72 
where the surface wave field wsurt is 
Usurt(X, y) are e ics ) e-ivk + any (51) 
a (—m) 
and the cylindrical field uy) is 
K py(8)  eithotn/4) 
Ueyl = igs ———— a) (52) 


at(kcosd) +/2rkp 


The power pattern P(?) of the cylindrical field is dis- 
tributed as a function of 3 by 
cos? + 1 


PO) = ; 
osd + V1 + (A/k)? 


| pw() |? : (53) 


It is interesting to observe that this expression 1s 
| pw) |? modulated by the previous power pattern func- 
tion of the line source. 
The Poynting flux II transported e each field can be 
readily calculated; we find first 
ne 
surf 


1 “Tou 
Tsurt = ae SS? 
WWE 9 Ox 


| E(—m) |?. 


(54) 


we lee 


In Appendix I we evaluate integrals of the form 


1B =f cos” 2 
0 


which materially eases the labor in computing (56) for 
any specified radiation pattern. 


cosd + 1 d 
cos? + /1+ A? 


bf 


APPENDIX | 
Evaluation of the integral Ip where M=1/1+(\/k)?: 


= cosd + 1 1 2x cosd + 1 
t= f ———— dj = f dd, 
0 cose? + M 2/79 cosd0+M 


| Af | Sie 


This integral can be evaluated very simply by residues; 
define the complex variable 


z= e”, 
whence 


dz = iz dd 


1 1 
cos dv = 5 («+ -), 
a z 


‘so that Ip becomes 


I 7 ec 2a 1a ide 
==> "i 
lel 27 + 2Mz+1 2g 
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‘The denominator can be readily factored 


fee Zits Vi (aa) (a9), 


“where 
eng a . g = us ie 1. 


A 


and since | M| >1, clearly |@| <1, || >1. If we add 


the residues at z=0 and z=a, we obtain 


- 
3 


a ie lee ate fo 1 
a an 
2 a(a — B) a8 


or in terms of A=X/k 


A-— 
? | Mesh |. 
: ea A) 


It is a simple matter to note that [)>—>0 as A>~, and 


Ij as A->0. 
If the integral is of the form 


a cost + 1 
ji af ces? t= — 0, 
0 cost + M 


the same procedure yields 


3 ra 
7 ee 
a : lela a 22+ 22+ 1 ( e) 


. ee a cee 
< Qri 2° f2|[=1 2? ae 2Mz Se 1 z 


Observe that ifn=1 


1 1 Ne il 
Zz Zz 2? 


1 LNe 2 
itna2—(s+—) eg? to = I, 
Z 2 z 


F 


efficient of 1/z is 


—) AOD ee 8 


and in general for ” odd the coefficient of 1/z vanishes 
for the factor (z+1/z)"/z, whereas for m even, the co- 


Hence for m odd, the only pole of the integrand is at 
2=a; we find 


AP A/T on? Ne 
I, = Ip Soe) n odd, 
Aiea / PE? 
ee (See Y aT'(n + 1) 
Sd (Crs CERT ey 
T2 
2 
n even. 


APPENDIX II 


First, observe that the cos m3 functions form a com- 
plete set on the interval 0S0S7. Expand p.(#) as a 
Fourier cosine series 


Pol®) = >> an cos nd (57) 
n=0 
in the interval 0<0Sz7 either from an analytic knowl- 
edge of p.(#) or numerically [3] from a measured radia- 
tion pattern. Truncate the series at some n= N and ex- 
press cos uv as a polynomial of order 1 in cos #, e.g., 


cos 43 = 8 cos‘? — 8 cos? 0 + 1. 


A collection of terms in (57) will yield 
a / 
pr) = D7 an cos" 8, (58) 
n=0 


or as a function of v=k cos 3 


N p\ 

Ey) = >o@ (=) (59) 
n=0 Rk 

The expression (59) corresponds to a Tchebycheff ap- 

proximation [3] of £..(v) in the interval |v/k| Si bya 

polynomial of order NV. 


ACKNOWLEDGMENT 


The author gratefully acknowledges the many helpful 
remarks and suggestions of Prof. S. N. Karp of New 
York University. 


REFERENCES 


[1] A. F. Kay, “Scattering of a surface wave by a discontinuity in 
reactance,” IRE TrANs. ON ANTENNAS AND PROPAGATION, 
vol. AP-7, pp. 22-33; January, 1959. 

[2] J. Bazer and S. N. Karp, “Propagation of Plane Electromagnetic 
Waves Past a Shoreline,” Inst. Math. Sci., New York University, 
New York, N. Y., Res. Rept. No. EM-46; July, 1952. 

[3] C. Lanczos, “Applied Analysis,” Prentice Hall, Inc., New York, 
N. Y.; 1956. 

[4] B. L. Van der Vaerden, “On the method of saddle points,” 4 ppl. 
Sci. Res. vol. B-2, pp. 33-45; 1951. 


508 


IRE TRANSACTIONS ON ANTENNAS AND PROPAGATION 


. 


September 


Surface-Wave Luneberg Lens Antennas* 


C. H. WALTER}, SENIOR MEMBER, IRE 


Summary—It is demonstrated that a surface-wave structure 
having isotropy in the plane of the structure can be operated as a 
lens. The index of refraction is given for a dielectric slab on a ground 
plane, a bed of metal posts on a ground plane and a parallel-plate 
structure with one plate filled with holes. The necessary radial varia- 
tions in index for several Luneberg-type lenses are given. These in- 
clude the conventional Luneberg lens, Gutman’s modification, the 
rim-fed Luneberg adapted to a spherical surface and a modifica- 
tion adapted to a spherical surface and collimation of the rays at an 
arbitrary angle with respect to the plane on which the lens is mounted. 
An approximate analysis of surface-wave lenses based on optics is 
described. 


in a medium of variable index of refraction re- 

sulted in a type of lens that has many applications 
in microwave antennas. In the spherical Luneberg lens, 
an image ona sphere in a region with spherically symmet- 
ric index is focussed diametrically opposite onto a con- 
centric inner sphere. Usually the outer sphere is taken 
to be at infinity. In the cylindrical lens the spheres are 
replaced by cylinders in a region having cylindrically 
symmetric index. The electromagnetic theory of the 


Ble reian s work by Luneberg! on the optics - 


Luneberg lens has been worked out by Jasik? for the 
cylindrical lens and by Wilcox*® and Tai‘ for the spher- 
ical lens. Eaton’ has considered some extensions of 
Luneberg-type lenses, and recently Morgan® presented 
a general analysis of the Luneberg lens. 

Another basic study that has been applied to micro- 
wave antennas is that of surface-wave propagation. An 
electromagnetic surface wave can be defined as an elec- 
tromagnetic wave that propagates along an interface 
between two media. The earliest work on this subject 
appears to be that of Sommerfeld,’ who discussed the 
propagation of a transverse magnetic surface wave 
along an infinitely long cylindrical wire of finite conduc- 


* Received by the PGAP, April 1, 1960. This work was performed 
under Contracts AF 33(616)-3353 and AF 33(616)-6211 between 
Wright Air Dev. Center and The Ohio State University Research 
Foundation. 

} Antenna Lab., The Ohio State University, Columbus, Ohio. 

*R. K. Luneberg, “The Mathematical Theory of Optics,” Brown 
University Press, Providence, R. I., pp. 208-213; 1944, 

*-H. Jasik, “The Electromagnetic Theory of the Luneberg Lens,” 
AF Cambridge Res. Center, Bedford, Mass., Rept. No. TR 54-121; 
November, 1954, 

3 C. H. Wilcox, “The Refraction of Plane Electromagnetic Waves 
by a Luneberg Lens,” Lockheed Aircraft Corp., Van Nuys, Calif., 
Rept. No. MSD 1802; June, 1956. 

“C. T. Tai, “The Electromagnetic Theory of the Spherical Lune- 
berg Lens,” Antenna Lab., The Ohio State University Research 
Foundation, Rept. No. 667-17; August, 1956. 

5 J. E. Eaton, “An Extension of the Luneberg-Type Lenses,” 
ae Res. Lab., Washington, D. C., Rept. No. 4110; February, 

6S. P. Morgan, “General solution of the Luneberg lens problem,” 
J. Appl. Phys., vol. 29, pp. 1358-1368; September, 1958. 

_ A. Sommerfeld, “Fortpflanzung Electrodynamischer Wellen an 
Tae Zykindrischen Leiter,” Ann, Phys, Chem., vol. 67, p. 233; 


tivity. Important contributions have been made by 


Cutler® with his work on electromagnetic waves guided | 
by corrugated conducting surfaces, Goubau*® with his” 
work on electromagnetic waves guided by a dielectric. 
coated wire, and Attwood!® with his work on surface- 
wave propagation over a dielectric coated plane con-_ 
ductor. A good summary and extensive bibliography — 
on surface waves have been presented by Zucker.™ | 

It is the purpose of this paper to show that a surface-_ 
wave structure can be made to perform as a two-dimen- 
sional Luneberg-type lens. By definition, the index of re- | 
fraction y of the lens is given by 


(1) 


where c is the velocity of light in free space, and v is the | 
phase velocity of a wave in the medium under consider- | 
ation. In this case v is the phase velocity of the surface — 
wave. Surface-wave structures consisting of both nat- 
ural and artificial dielectrics will be considered. 


n = c/, 


SURFACE-WAVE STRUCTURES 


A surface-wave Luneberg lens is a two-dimensional 
lens. Lens action is achieved by varying the phase veloc- 
ity of the surface wave in the proper manner. The sur- 
face-wave structure must be isotropic in the plane of the 
structure. 

Examples of surface wave structures having the nec- 
essary isotropy are dielectric sheet on a ground plane, 
bed of metal posts on a ground plane, and parallel plate | 
structure with holey top plate. 


A. Dielectric Sheet on Ground Plane 


This problem has been treated by many au- 
thors.'°-® A sketch of the dielectric sheet structure is 
shown in Fig. 1. Both transverse magnetic (TM) and 
transverse electric (TE) surface waves are possible. A 
TM mode will produce polarization perpendicular to the 
ground plane, whereas a TE mode will produce polari- 


*C. C. Cutler, “Electromagnetic Waves Guided by Corrugated 
Conducting Surfaces,” Bell Telephone Labs., New Yorks Nias 
Rept. No. MM-44-160-218; 1944. 

*G. Goubau, “Surface waves and their applications to trans-_ 
nig os lines,” J. Appl. Phys., vol. 21, pp. 1119-1128; November, 

0S. S. Attwood, “Surface wave propagation over a coated plane 
conductor,” J. Appl. Phys., vol. 22, pp. 504-509; April, 1951. 

uF. J. Zucker, “The guiding and radiation of surface waves,” 
Proc. Symp. on Modern Advances in Microwave Techniques, Poly- 
pee Inst. of Brooklyn, Brooklyn, N. Y., pp. 403-436; November, 

CC. H. Walter, “Surface-Wave Luneberg Lens Antenna,” An- 
tenna Lab., The Ohio State University Research Foundation, Rept. - 
No. 667-32; May, 1957. 

8 J. H. Richmond, “Surface Waves on Plane Dielectric Sheets 
and Sandwiches, ” Antenna Lab., The Ohio State University Research 
Foundation, Rept. No. 786-10; September, 1958. 
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Ground Plone 


Fig. 1—Infinite dielectric sheet of thickness ¢ and dielectric constant 


= « onan infinite perfectly conducting ground plane. The sheet is 
assumed to be lossless, isotropic and homogeneous. 


he 


zation parallel to the ground plane. The index depends 


‘on the mode, the thickness, ¢, of the sheet, and the rela- 


“tive dielectric constant. For TM modes, the index for a 


ie 


sheet of dielectric of infinite extent is obtained from the 
equations 


Bo? = ko? + kz’, (2) 
€Bo® = ki? + ka’, (3) 
Band 
e-ko = ki tan yi, (4) 
where 


ate toe hk od Ae 


—Bo=2m/do = free space propagation constant, 

2 do =free space wavelength, 
ky = propagation constant in z direction in €, 
ki = propagation constant in z direction in 4, 
kz = propagation constant in x direction, 


€. = €,/€0. 


Eq. (4) may be obtained by the usual boundary value 


solution or by application of a transverse resonance 
- method." 


“2 


ee 


RTE SA 


hn i ihn Sas 


Eqs. (2) and (3) reduce to 


Ro\2 
eae Fan (3) (5) 
Bo 
and 
By 
=e | (6) 
: & 
respectively. 
Eliminating ko gives 
——4 7 
n= y/1+(— tan bi) ; (7) 
Bo€r 


which is a convenient form for comparison with the 
~ metal post structure to be described later. Elimination 


LS 


of both Ro and hk; gives 


Pata Veseialves el © 
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If the required index is known, the thickness ¢ is 
given by 


t/ro = tan! [a = | (9) | 


Vea Voliveoum 
The use of (9) assumes that the variation in ¢ will be suf- 
ficiently gradual so that at any point on the lens the in- 
dex is the same as it would be for a sheet of that thick- 
ness and dielectric constant and of infinite extent. 

For TE modes the index is obtained from (2) and (3) 
and 


Ro = ky cot kit. (10) 


Again, if the index is known, the thickness ¢ is given by 


fected 
COG aS et 
Vv €; oul 7? 

Higher-order modes are possible for both TM and 
TE propagation. The lowest-order TM mode degener- 
ates into a TEM wave as ¢ goes to zero. The lowest-order 
TE mode, however, has a cutoff. For a given dielectric 
there is a critical value of t/Ao, below which TE surface- 
wave propagation is not possible for the structure in 
Fig. 1. This does not mean that the structure of Fig. 1 
is not suitable for a TE surface wave Luneberg lens; it 
merely means that such a lens would have finite thick- 
ness at the rim. 

In practice it may be most convenient to keep ¢, con- 
stant and obtain the necessary variation in 7 by varying 
t. However, this makes 7 a function of Xo. Satisfactory 
operation of a 45 \» diameter lens has been obtained 
over a 10 per cent band in this manner. Operation over 
greater bandwidths may be achieved by increasing ¢ to 
reduce the effect of No and varying €, to obtain the nec- 
essary index. 


t/\o = (11) 


2av/er — 7" 


B. Bed of Metal Posts on a Ground Plane 


The metal post structures illustrated in Fig. 2 may be 
considered as sheets of artificial dielectric. Such struc- 
tures will support the dominant TM _ surface-wave 
mode.45 Furthermore, these structures are very nearly 
isotropic. The index depends on the post size and spac- 
ing and the height 4. An exact analysis is not readily 
obtained. An approximate analysis based on a trans- 
verse resonance solution gives’® 


n= V1+ (W tan kh)’, (12) 


4 A. F, Kay, “Excitation Efficiency of Surface Waves Over Cor- 
rugated Metal and Doubly Corrugated Metal and in Dielectric 
Slabs on a Ground Plane,” Technical Research Group, New York, 
N. Y., Scientific Rept. No. 5; December, 1956. ; 

15H. B. Querido, “Surface Wave Fields and Phase Velocity 
Variations of Grounded Dielectric Sheets and of Periodic Structures 
of Metal Posts on a Ground Plane,” Antenna Lab., The Ohio State 
University Research Foundation, Rept. No. 667-46; November, 1958. 
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Fig. 3—Holey-plate surface-wave structure. 


SQUARE POSTS field would be everywhere normal to the common walls, 
these walls may be removed to form the structure in 

(a) Fig. 3. An expression for the index may be obtained by 

using a transverse resonance procedure.'® The resulting 


Y : : 
equations are given by 


12bc ; 
bl? — (e, — 1)Bo?] — 7 Ses — j2« cot xa = 0. (13) 
and 
kz = Ve6o? — x = B — ja, (14) 
where 


k=propagation constant in z direction in parallel 
plate region 


Tv 
= = e) 
a 
CYLINDRICAL POSTS 
; k. B =Bon, 
(b) a = attentuation constant of wave traveling in x direc- 
Fig. 2—Metal post surface-wave structures. tion. 


where k&1.058.=propagation constant in the direc- By a perturbation method the index is given aDPrOx 


tion parallel to the posts in the post region, mately Dy ‘ ere 
: 2 
W = weighting factor ie {1 +— |, (15) 
. 47a? 
aren for square posts where 
gtw No 
PD=—D/4 San Macys 
Serta a LO circulamposts, /« Sa (5°) 
PD 2a 
Eq. (12) is found to be a good approximation for struc- 3 BB 
tures with five or more posts per wavelength. AS G? + B? 
C. Holey-Plate Structure G’ = normalized aperture conductance 


A holey-plate surface-wave structure capable of pro- 
viding an index from nearly zero to a value approaching = be [x2 — (e, — 1)Bo?] 
the square root of the dielectric between the plates is il- ‘ 
lustrated in Fig. 3. This structure gives horizontally normalized aperture susceptance 
polarized radiation (electric field parallel to the sur- 6bc 
face). An approximate expression for the index may be Ber 
obtained by assuming that the holey-plate structure is i 


composed of stacked 1 ee waveguides with holes in a % L. O. Goldstone, and A. A. Oliner, “Leaky Wave Antennas I: 


: ‘ Rectangular Waveguides,” Polytechnic  Instit t f Brookl 
narrow wall of each of the guides. Since the electric Brooklyn, N. Y., hoot No. R-606-87, PIB 534; AoE 1957. a 


w 
Tl 


7 
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The structure is assumed to be lossless. However, for 
indexes less than unity a wave supported by such a 
structure is attenuated due to radiation. The attenua- 
tion constant in nepers per wavelength is given by 


ee Rt 


a 


where 
eee ee 
G+ B? 


(16) 


_ Eqs. (15) and (16) become very poor approximations 
_ for operation in the region near cutoff. For such opera- 
_ tion it is necessary to find an exact solution for (13) and 


(14). In all cases the use of the above equations is re- 


“stricted to c<d/2 and d<b. For isotropy, 6 is made 


equal to c. 


LENS DESIGN 


The surface-wave lenses that have been considered to 
_ date are of the Luneberg typehavingaradiallysymmetric 
_ index. A typical rimfed lens is shown in Fig. 4. The di- 
ameter of the lens will be determined by the desired 
bs azimuth beamwidth (lens assumed to be horizontal). 
_ The sidelobe level in the horizontal plane may be con- 
trolled by proper choice of feed pattern.’’ The vertical 
‘ beam is not easily controlled. Its shape and position are 
influenced by the shape and extent of the lens and 
- ground plane as well as the location and extent of the 
~ lens feed. For surface-wave lenses with indexes equal to 
- or greater than unity, the vertical pattern is essentially: 
that of an end-fire line source with a nonuniform phase 
velocity. The effect of ground plane shape and extent is 
_ yery pronounced for the case of vertical polarization.'® 
For lenses 15 to 60 wavelengths in diameter and ap- 
_ proximately equal extent of ground plane, the vertical 
_ beam maximum ranges from 12° to 6° with respect to 
the ground plane, respectively. Higher elevation angles 
for the vertical beam maximum may be obtained by 
_ using lenses with indexes less than unity. One may also 
control the vertical beam position by partially covering 
_ the aperture to form a rim radiating lens as illustrated 
- in Fig. 5. In principle, the vertical beam position for a 
ring radiator may be controlled by means of the phase 
distribution around the ring. This would be achieved by 
the proper selection of index. The radial variations in 
index for various lenses are given below. 


us fess D. M. Peeler and D. H. Archer, “A Two-Dimensional 


Microwave Luneberg Lens,” Naval Res. Lab. 


Rept. No. 4115; March, 1953. 


, Washington, D. C., 


18 B, T. Stephenson and C. H. Walter, “Endfire slot antennas,” 


IRE TRANS. ON ANTENNAS AND PROPAGATION, 
April, 1955. 


vol. AP-3, pp. 81-86; 
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Fig. 4—Sketch showing a TM surface-wave lens on a ground plane. 


Shield 


Feed Contoured Lens 


Ground Plane 
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Fig. 5—Sketch of a rim radiating lens. 


A. Conventional Luneberg Lens 

In the conventional Luneberg lens an incident plane 
wave is focussed at a point on the rim diametrically op- 
posite from the incident wave. The index is given by 
the relation 


AG NN NOR te (17) 


where r is the normalized radius. Eq. (17) is not the only 
index variation producing focus at the rim,® but it is 
perhaps the most practical. 
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B. Luneberg Lens with Internal Focus 


There is no unique index variation for focussing a 
plane wave at a point within the lens. The problem has 
been considered by many authors.®§!9,20 One of the 
simplest relations is that obtained by Gutman?® in 
which the index for a lens of unit radius with focus at 
r1 is given by 


n(r) = : (18) 


C. Spherical Cap Lens 


A solution to the problem of a lens capping a spher- 
ical dome on a ground plane has been obtained.2! Such 
a structure is shown in Fig. 6. However, the analysis 
will apply also for a spherical depression, such as the 
one shown in Fig. 5. A plane wave traveling parallel to 
the ground plane and striking the spherical surface will Focus 
focus diametrically opposite at the intersection of the 
sphere and ground plane. The index variation is given 


by 


Surface- Wave Structure 
Capping Spherical Dome 


Incident 


Wave 


AV2R = 27 Parallel Plat 

Bs arallel Plate 
n(r) R? , (19) Transmission Line 
Capping Spherical 


Dome 
Foc 
wheré g ae Incident 
Wave 


1 
A=1i+ We qs Fig. 6—Spherical cap lenses. 
TE. 


ET EN ee ’ where 
R=1+4/1-(4). 
a p= nr 


In the limit as the radius of the sphere becomes infinite 1 i : -1 COS) C Fs es ies a 
‘ = tan eeneeeeeeet (| 3 
(19) reduces to (17) for the conventional Luneberg lens. Td rs p?(r) — cos? rf/a—r ? 
and P(r) is arbitrary except for the restrictions 


Q(p) 


I 


If the lens is covered to restrict radiation to the rim 
as in Fig. 5, an index variation for a feed at m<1and 
rays collimated at an angle £ with respect to the ground 


plane may be obtained from the analyses in footnotes 6 P(r) = p(r) > cosé Tri 
and 21. A resulting expression that reduces to Lune- 
berg’s relation [see (17)] for a= ©, 4,=1 and £=0 is 7 1 


IV 


1 
; : See cost f d 
given by 2 rn 1Va* — r/p2(r) — cos? £ : 


_ 60s E(cos £ + +/cos? — — p?) (a + Va — 712)? + pr ,2¢22() 
2arr/cos t(cos + Veos*§ = pi) a+ Vato riea@) a 


19 J. E. Brown, “Microwave wide-angle scanner,” Wireless E 
vol. 30, pp. 250-255; October, 1953. ; as 


MEASUREMENTS 
20 A. S. Gutman, “Modified luneberg lens,” J. Appl. Ph Surface- 
25, pp. 855-859; July, 1954, Ppl. Phys., vol. e-wave lenses have been constructed for both 


zs Adachi and C, H. Walter, “A Spherical Cap Lens,” Antenna aM on ie ope Alan X-band, Ky-band and ie 
re The Ohio State University Research Foundation, Rept. No. band with diameters from 15 to 60 wavelengths in order 
3-3; April, 1959, to show that lens action is indeed achieved, When the 
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lens is fed by a point source at the focus, the beam- 

width and sidelobe level are very close to the expected 

values. A point source feed, however, may produce an Deleciri 
excessive interference pattern in the vertical plane. For a 
TM dielectric-sheet lenses a tapered-depth source!® as 

illustrated in Fig. 7 has been found to work quite well. 

The aperture of the antenna in Fig. 7 has the necessary 

field components for TM surface-wave excitation and 
the velocity ratio c/v (the index) along the aperture 


ene Tapered Depth 


Ground Plane 


Rectangular 


“may be controlled by varying the width and depth. The wen ees 
t ol 
antenna may be placed under the lens and its phase pears 
_ velocity adjusted to match that of the lens at the point Fig. 7—A tapered-depth antenna suitable as a feed 
of focus. For the TE case, however, a comparable ta- element for a TM surface-wave lens. 


_pered-depth feed is equivalent to a magnetic current 
_ source and produces a strong cross-polarized component 
_ in the far-field pattern of the lens antenna. 


Lens. 
Patterns of a 60 \,) diameter TM dielectric sheet lens BES ae 


of the conventional Luneberg type with tapered-depth Topered on 


Plane 


feed are shown in Fig. 8. Polystyrene dielectric was used age 


and the necessary thickness vs radius was obtained 
from (9) and (17). Measured patterns for a 15 Xo di- 
~ ameter lens with feed at 7:=0.65 are shown in Fig. 9 
(next page). This lens was designed from (9) and (18). 
_ Measured patterns of a metal post lens designed from 
~ (12) and (17) are shown in Fig. 10. An open ended wave- 
_ guide positioned at the rim of the lens provided the 
_ feed in this case. 
These results and additional measurements which in- 
cluded TE lenses demonstrate that surface-wave struc- 
- tures can be made to operate as lenses. The vertical pat- 
~ tern (plane normal to ground plane) is that of a traveling 
wave source with a nonuniform phase velocity and 
_ with significant influence from the ground plane for TM 
_ modes. Cross sections of the beam show the lens action. 
- Although there is some defocussing due to the use of an 
extended source the lens operation compares quite 
favorably with that of conventional two-dimensional 
_ lenses.!” In some cases there may be significant deteriora- 
tion in lens action due to direct radiation from the feed. 


Voltage in db 
‘ 


ANALYSIS OF SURFACE-WAVE LENSES 
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- The radiation patterns of a surface-wave lens may be 


Horizontol Or Lens Pattern Vertical Pattern 


_ obtained by application of Huygens’ principle. By Through 0-0 Of Vertical Pottern 
Huygens’ principle the electromagnetic waves generated Fig. 8—Measured patterns at 35 kmc of a 60 \ diameter TM surface- 
in a source free region by any set of sources can be gen- wave Luneberg lens with 2 o tapered-depth feed. The ground 


i a x i bey he lens. 
erated by a sheet of electric currents of surface density plane extends approximately 60 do beyond the lens 


J="nxXH and a sheet of magnetic currents of surface 

‘density M@ =EX2 spread over the surface S which sur- i 5 ,-ikR gS 

rounds the sources, where E and H are the electric and 7 =f f Sa evel (21) 
magnetic fields, respectively, and is a unit vector nor- 4nd s ae 

mal to S and pointing into the source free region. Two 

vector potentials may be defined in the usual manner,” and 


] 


; 


2 A 8B. Bronwell and R. E. Beam, “Theory and Application of dey Me-**dS (22) 
Microwaves,” McGraw-Hill Book Co., Inc., New York, N. Y., ist R= ne aaa he 
ed.; 1947. TJs 
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HORIZONTAL PATTERN 
(cross SECTION b) 
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100° = 90° Cet Cs 
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Fig. 9—Measured patterns at 10 kmc of a 15 \) diameter TM surface- 
wave Luneberg lens with Gutman modification. A 2 Xo tapered- 
depth feed was matched to the lens at 7, =0.65. 


Cross Section a—a 


es a ee te at 10 kmc of a 15 Xo diameter cylin- 
rical metal post lens. e ground plane extend i 
15 \o beyond the lens. ee pee 
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where 


k=propagation constant, 

s= permeability, 
e=dielectric constant, 

R=distance from a point on S to point of observation. | 


The electric field on the source free side of Sis given by” 
; 

at 1 ad ke new , 

Ea NK Fey eee 


(23) 
€ We 
where w is the angular frequency of the source. © | 

It is usually most convenient to take S as the surface | 
of the structure and in most cases a satisfactory approxi-_ 
mation is obtained by assuming that the currents on the 
rear surface are negligible compared to those on the 
front surface where the lens is located. 

Tangential E and H over the lens and ground plane 
may be found approximately by optical methods. Re- 
ferring to Fig. 11, one can assume that all of the energy — 


Fig, 11—Coordinate system used in approximate analysis 
of surface-wave Luneberg lens antenna, 


in an incremental angle Ao at the source, which is as- 


sumed to be located at a point P;, travels as a surface 


wave between two surfaces designated as rays.in the 
figure. If the lens is focused at infinity, the angles r and 
o are equal. From Fermat's principle and the applica- 
tion of the calculus of variations to a radially symmetric 
medium, the ray paths are found to satisfy the relation 


nr sin y = sin r. (24) 


This is a generalization of Snell’s law. An equivalent 
equation in differential form is 


1 /dr\2 
nr = 4/3 + = (+) sin 7. 
r\do 


(25) 


. 
, 


a 
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The optical approximation enables one to find the 
relative surface-wave amplitude, from which the rela- 
tive tangential electric fields, E;, and magnetic fields, 
_ Hf, can be found for a particular surface wave structure. 
The relative phase, y, of the fields at any point P on the 


‘surface with respect to the feed point P, may be ob- 


tained from the relation 


P 
y = 360r f nds, 


Py 


(26) 


where ds is an element of path length along the ray 


through P; and P and n is the radius of the lens in 
_ wavelengths. 


If only the lens pattern of a slow wave (7>1) lens is 


2 desired, a good approximation is to replace the lens by 
~ an equivalent line source such as line 0 in Fig. 11. As- 
_ suming no g variation (direction out of page) conserva- 


VON" Le Re 


tion of energy gives 


U(a)Ao = Us(y)Ay, (27) 
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where 


U(c) =energy distribution per unit ul in the pat- 
tern of the feed, 
U>(y) =energy distribution per unit length along 0. 


If the lens is focused at infinity the phase is constant 
along b and 
E*(c) E*(c) 


EQ) = = ee 


28 
mri cosa) Vn? — y oy 


where 


E(c) =electric field pattern of the feed, 
E,(y) =electric field amplitude along 8, 
i =index at point of feed. 


If the lens is used as a ring radiator, the amplitude 
and phase at the rim of the lens may be found by the 
above methods. 
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Correction 


E. J. Martin, Jr., author of “Radiation Fields of Cir- 
cular Loop Antennas by a Direct Integration Process,” 
which appeared on pp. 105-107 of the January, 1960, 
issue of these TRANSACTIONS, has called the following 
to the attention of the Editor. 

In (1), the subscripts P should be lower case, in order 
to be consistent with the text; thus, (1) should read 


r = fp — asin 6, cos (d — $7). 


In (2), (4), and (6), the y components at the vectors 
should be proportional to sin 6), rather than i/sinds. 

In the second sentence of the last paragraph, the sym- 

bol V appears inverted in three places. These relations 


should read 


E = (A/joue) V(V- A) — jwA 


and 


H = 1/u(V X A). 
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Phenomenological Concepts in Beyond-the-Horizon 
Tropospheric Propagation* 


B. J. STARKEY} 


HE PAPER by Bullington! gives an excellent 
“| Cape survey, but one rather puzzling 
aspect of the author’s approach is a somewhat 
confusing vagueness with which he treats the phenom- 
enological issues involved. 
A reader anxious to interpret the relevant part of the 
paper in a more definite way would have to draw the fol- 
lowing conclusions: 


1) The theory generally used by radio engineers in 
their beyond-the-horizon tropospheric propaga- 
tion work is the one based on the assumption that 
scatter of EM waves caused by atmospheric 
“blobs” is the principal mechanism involved. 

2) It appears that a so-called “reflection” approach 
might provide a rather easier explanation of ex- 
perimental results. 

3) Nevertheless, let us forget about the physical 


* Received by the PGAP, April 5, 1960. 

} Dept. of National Defence, Ottawa, Ont., Can. Formerly 
with RAE, Farnborough, England. 

‘K. Bullington, “Status of tropospheric extended range trans- 


mission,” IRE Trans. oN ANTENNAS AND PROPAGATION, vol. AP-7, 
pp. 439-440; October, 1959. 


mechanism of the phenomena involved, keep on 
using equations “that fit the experimental results” 
(which in this context seems to mean equations of 
the “scatter” theory conveniently manipulated to 
“come into reasonable agreement”- with experi- 
ment) and “minimize the fine shades of meaning 
in the interpretation of the equations.” The his- 
torical fact that Isaac. Newton discovered a for- 
mula that accounted for the motion of the planets, 
but did not explain the mechanism of gravitation, 
is quoted as a good precedent for such attitude. 


Of course, Newton explained a large number of com- 
plex phenomena by means of a single one, gravity, 
which was left unexplained. It is rather difficult to place 
tropospheric propagation in the same category of basic 
concepts as gravitation. An unreserved application of 
Dr. Bullington’s recommendation in this respect to 
other fields of scientific endeavor might lead to rather 
unexpected results. 

The paper also tends to create an impression that the 
difference between the concepts of “reflection” and 
“scattering” can be regarded simply as a difference be- 
tween the “layman” and the “physicist” interpretation 


September 


‘~ 


1960 


of the same phenomenon. Certainly, there is more to it 
than that! Physical factors behind this difference in the 


case of tropospheric propagation have been discussed 
_ elsewhere.” 


Reply to Starkey’s Discussion? 


I agree with Mr. Starkey that an understanding of 
the physical mechanism is a worthwhile objective. How- 
ever, until that goal is reached, it seems important to 
describe the experimental results with a minimum of 
- theoretical interpretation. 


F * B. J. Starkey, W. R. Turner, S. R. Badcoe, and G. F. Kitchen, 
“The effects of atmospheric discontinuity layers up to and including 

the tropopause on beyond-the-horizon propagation phenomena,” 

_ Proc. IEE, vol. 105, pt. B, suppl. 8, pp. 97-105; January, 1958. 

3 Received by the PGAP, April 20, 1960. 
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Since my paper! was the opening address of a sym- 
posium, its principal purpose was to describe the general 
characteristics and usefulness of tropospheric beyond- 
horizon transmission without taking sides on the con- 
troversial question of the theoretical explanation. The 
vagueness to which Starkey refers was caused in part 
by the neutral purpose of the paper and in part by my 
own reluctance to accept the scattering theories because 
of the arbitrary parameters involved. The generalized 
equation used in the paper was intended to represent 
experimental data. The use of the angle @ does not nec- 
essarily imply scattering, since it is possible and in 
some cases convenient to describe diffraction and reflec- 
tion, as well as scattering, in this way. 

K. BULLINGTON 
Bell Telephone Labs. 
Murray Hill, N. J. 


The Magneto-Ionic Theory for Drifting Plasma* 
H. UNZt 


later date. However, it is felt that one of the main 

results should be mentioned as early as possible. 

It has been found that the refractive index of the mag- 

neto-ionic medium, for the case of plasma drifting witha 

constant velocity %o, is given by the following algebraic 
- equation: 


A PAPER titled as above will be published on a 


C[A? — Yin? — 1)?| 
= AV en? = 1)? + Br? X*nr(1 — nB1)?| 
— 2XV1(Br-¥r)n(1 = Bx) (0? — 1)? = 0; (1a) 


A= (n? — 1)(1 — iZ) + X(1 — 081); and (1) 


GC = (n? — 1)(1 — X = iZ) + w XB’, (1c) 


* Received by the PGAP, April 25, 1960. This work has been 
supported in part by contract No. DA-44-009-ENG-3769 at the 
Midwest Research Institute, Kansas City, Mo. 

+ Elec. Engrg. Dept., University of Kansas, Lawrence, Kans. 


where 


) voT and VOL 


are the transverse and the longitudinal components of 
the constant drift velocity vector of the plasma. All the 
rest of the symbols are defined as in the classical mono- 
graph by Ratcliffe." 

In the particular case of stationary plasma, 4.€., 
Br =B.=0, the solution of (1) will give the classical 
Appleton-Hartree equation’ of the magneto-ionic theory. 
Several other more complicated cases have been solved 
explicitly from (1). The above theory has been applied 
to emissions generated in the earth exosphere,” whistler 
propagation,’ etc. The results will be published soon. 


1]. A. Ratcliffe, “The Magneto-Ionic Theory and Its Applications 
to the Ionosphere,” Cambridge University Press, Cambridge, 
England; 1959. 

2R, M. Gallet, “The very low frequency emissions generated 
in the Earth Exosphere,” Proc. IRE, vol. 47, pp. 211-231; February, 
1959. 

3R. A. Helliwell and M. G. Mogran, “Atmospheric whistlers,” 
Proc. IRE, vol. 47, pp. 200-208; February, 1959. ; 

4, E. Barrington, “Whistler propagation in the presence of ion 
Streams,” (Abstract) IRE TRANS. ON ANTENNAS AND PROPAGATION, 
vol. AP-7, pp. 287; July, 1959. 
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An Endfire Array Continuously Proximity-Coupled 
to a Two-Wire Line* 


GEORGE R. FORBEST{ 


N THE new type of VHF-UHF communication an- 

tenna designed by Sletten and Forbes in 1957,1 

radiating dipoles are capacitively coupled to a two- 
wire transmission line. The degree of coupling is a func- 
tion of the angle that the dipole makes with the axis of 
the transmission line, and is maximum at an angle of ap- 
proximately 20°. The broadside antenna developed with 
dipoles spaced \/2 apart along the axis of the two-wire 
line and alternately making an angle of +20° and —20° 
with the axis is excellent for VHF communication ap- 
plication, since it requires no individual dipole feed 
lines. 

The principle of this broadside antenna has been ap- 
plied to develop an endfire array of proximity-coupled 
dipoles. Fig. 1 shows \/2 radiating dipoles, spaced \/4 
apart, capacitively coupled to a two-wire line. Each di- 
pole is transversely perpendicular to the axis of the line. 
When the coupling between the dipoles and two-wire 
line is small (dipoles spaced \/4 or less apart), the phase 
velocity along the line approximates that of free space 
and the currents induced in the dipoles produce a single 
beam in the direction of propagation of the traveling 
wave. (Spacing of \/2 produces a bidirectional pattern.) 

Although the coupling of a dipole perpendicular to 
the axis is considerably less than that of a dipole inclined 
to the axis, it is sufficient for use with an array 7.5 \ long 


Fig. 1—)/2 radiating dipoles spaced \/4 apart capacitively 
coupled to a two-wire line. 


* Received by the PGAP, April 12, 1960. 

} Electromagnetic Rad. Lab., AF Cambridge Res. Ctr., Bed- 
ford, Mass. a 

1C, J. Sletten and G. R. F orbes, “A New Antenna Radiator for 
VHF-UHF Communications,” AF Cambridge R. Cc 
No. AFRC-TR-57-114; June, 1957, Tear 


(30 dipoles spaced \/4 apart).? The coupling depends on © 


the spacing between the two wires of the line, the length 
of the dipole, and especially on the height of-the dipoles 
above the line. 

In testing the admittance of a single dipole shunted 
across the transmission line, it was not possible to ob- 
tain a purely resonant load. This condition was also 
found by Seshardri and lizuka.’ 

Fig. 2 shows the measured conductance of a \/2 di- 
pole plotted against the height of the dipole above the 
transmission line. To determine impedance characteris- 
tics, an equivalent circuit was derived for the dipole 
[see Fig. 3(a) ].4 Each wire couples to the"dipole_by a 


DIPOLE LENGTH 9 IN. 
FREQUENCY GOO MCPS 


o 
0 125 250 .378 500 625 .750 .875 1.0 50 150 178 2.25 2.50 278 30 


Le 20 
= [HEIGHT OF DIPOLE aBove TRANSMISSION LINE IN INCHES ] 


Fig. 2—Reciprocal of VSWR vs height of dipole 
above transmission line. 


a 
Oe 


i? de 


Fig. 3—Equivalent circuit of a dipole coupled to a two-wire line. 


*G. R. Forbes, “An Endfire Array 
Coupled,” AF Cambridge Res. Ctr, 
59-368; December, 1959. 

SSE IRS Seshardri and K. lizuka, “A Dipole Antenna Coupled 
Electromagnetically to a Two-Wire Transmission Line,” Cruft 
Nac ae oo Cambridge, Mass., Tech. Note No. 

~tN-90-994, Contract AF19(604)-786, Scienti 
21; October 3 ios, (604) cientific Rept. No. 
ig. 3(a) and 3(b) and the mathematical analysis are due t 
J. Ruze, Consultant to the Electromagnetic Radiation Paberaion a 


Continuously Proximity- 
Tech. Rept. No. AFCRC-TR- 
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capacitive reactance X,/2. The two capacitors are 

shorted by an inductive strip created by the radiating 

dipole, the impedance Xz of the dipole coupling electro- 
magnetically with the strip. With this coupled imped- 
ance denoted as Rp+jX p, the equivalent circuit is as 
shown in Fig. 3(b). The input impedance Z is 


iX7(R iX 
= — 5x, 4 7 Beale! D) 
3 Z RpX 1? 

. R*p + (Xp + Xx)? 

Se 


[R2p + Xv(Xv + Xz)] 


aye 
wan Fen (xpath) 


When X,>X, (small coupling and short strip), 


RpX?1 


Ea Xs, 

R?’p + (Xp = Xy)2 
j The series resistance is then small and the reactance ca- 
' pacitive; reactance cannot be tuned out by changing 
the dipole length. 
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Fig. 5 shows typical radiation patterns of the array, 
both terminated and unterminated. Patterns measured 
over a range of frequencies from 580 to 670 mc showed 
very little change in their basic characteristics. From 
the results of these measurements, it is apparent that 
the effective phase velocity along the array is slower 
than the velocity of the light since the Hansen-Wood- 
yard condition is approximately fulfilled. 

Radiation from the line, checked by removing all the 
dipole elements, is about 17 db below the beam peak 
when the line is terminated. The level can be lowered 
by decreasing the spacing between the wires. 
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When the impedance of the entire array was meas- 
ured, it was found that the array was not coupling as 
much energy as expected, probably because mutual 
coupling between parallel dipoles is high in an endfire 


marray. — 


At an operating frequency of 600 mc the most prac- 
tical length for an antenna of this type is the 7.5 A (13- 
- foot) array constructed of 30 dipoles (see Fig. 4). The 

amplitude distribution is effectively tapered by varying 


the height of the dipoles above the two-wire line. 


——— 
if ON eves 
Tet 


Fig. 5—Thirty-element endfire array with 25-db taper. 


Arrays of the type in Fig. 1 were used to study 
“modulated” surface wave antennas,® in which two sur- 
face waves—one along the two-wire line and one along 
the array of dipoles—are coupled in an attempt to con- 
trol the aperture illumination. This approach differs 
from the one described in both design philosophy and 
results. 


5D. K. Reynolds and R. A. Gigelmann, “Research on Traveling 
Wave Antennas,” AF Cambridge Res. Ctr., Tech. Rept. No. AFCRC- 
TR-59-160 (ASTIA Document No. AD217821), Final Rept. Contract 
AF19(604)-4052; June, 1959. 
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Direction Finding in a Two-Component Field* 


D. N. TRAVERS}, MEMBER, IRE, AND W. M. SHERRILL}, MEMBER, IRE 


finding is sometimes severely complicated by re- 

radiated secondary signals of frequency equal to 
that of the primary, but of different azimuth, ampli- 
tude, and phase. An example of considerable practical 
importance is naval shipboard direction finding. A 
method for analyzing a two-component field, using a 
twin-channel receiver and crossed loops, has been pro- 
posed by Gabler.! The method makes use of the fact 
that the phase difference between the primary and sec- 
ondary components varies throughout the resultant 
field, and provides a means of displaying both azimuth 
and relative amplitudes of the components incident on 
the antenna in the form of an envelope of ellipses. The 
purpose of this note is to show that the use of a twin- 
channel receiver and indicator is not essential to achieve 
a display which provides the same information for a two- 
component field. A spinning loop antenna, or crossed 
loops and goniometer with a single-channel receiver, 
can be used to analyze a two-component field and dis- 
play both the azimuths and the relative amplitudes of 
the two incident waves. 

Consider two spinning simple loop antennas, A and 
B, of radius r, small with respect to the wavelength of 
two incident plane waves arriving from infinitely dis- 
tant sources. Let the antennas be located about a co- 
ordinate center, point O, from which phase in the re- 
sultant field of the incident waves is measured. The 
radial distances, p4 and pz, from the origin to the loops 
are arbitrary. 

Let two vertically polarized plane waves of equal fre- 
quency, but of arbitrary amplitudes EZ; and E2, phases 
g: and ¢», and azimuths 6; and 6, arrive at horizontal 
incidence upon the spinning loops as shown in Fig. 1. 
Since the sources of the incident waves are remote, their 
azimuths are the same for each of the two spinning loop 
antennas. Thus the only difference between the signals 
induced in the two loops is a phase difference. 

The total signal (Vr) induced in either of the loops 
may be derived in the following manner: 


[: is well known that high-frequency radio direction 


Ver = jhe, Ey sin [Bor sin (04.7 — 51) 


- ei let—Bopa,B cos (61—-7A,B8)+41] 


-++ EK, sin [Bor sin (4,2 — 6.)| 


» ei l@t—Bopa,B cos (2-742) (1) 


* Received by the PGAP, April 12, 1960. 

} Dept. of Electronics and Elec. Engrg., Southwest Res. Inst. 
San Antonio, Tex. 

1H. Gabler, and M. Wachtler, “A new method of determining 
the components of radio bearings from coherent waves,” Elektrotech. 
Z. (ETZ), vol. 79A, pp. 385-388; June, 1958. National Research 
Council of Canada Technical Translation 819. 


Fig. 1—Geometry for a single channel direction 
finder in a two-component field. 


where h, is effective height, £ is field strength, 6+ wt is 
the field phase at point O, and By)=w/c. 

Since the wavelength is much greater than the loop 
dimensions, it can be assumed that 


sin [Bor sin (04,3 — 51,2) | = [Bor sin (64.2 — 51,0) | 
and 


Vor = = jheBor\ Ex sin (04,3 mee 5;)e—7 [Bor .B cos (61—TA,B)—¢1] 


+ E.sin (04.2 as 52) e—7 [Bora ,B cos (32-74,B)—o2] } eft, (2) | 


Let the relative phase between the two component volt- 
ages induced in either loop be dac4,), so that 


Vr= — jh.Bor\ Ey sin (O4.n — 63) 

+ Ey sin (64,2 — 52)e%al4.» } 

-eiletBopa,b cos (bi-r4,B)+6il (3) 
where 
a(4,B) = Bopa,a|cos (6: — t4,8) — cos (52 — Ta,n)| 
+ (¢2 — $1). 

The two patterns displayed by the two loops are 
| V7’ la = [E," sin? (04 — 81) + Ey? sin? (04 — 52) 

+ 2£,F» sin (64 — 61) sin (04 — 83) cos daa] *!2, 
| Vr'|s = [E:? sin? (63 — 8;) + Ey? sin? (0p — 8s) 

+ 2EF2 sin (82 — 61) sin (8g — 82) cos daz|'/?, (4) 


a / 4 . 
where | Vr'| includes the constant —jh.Bor and the 
transmission characteristics of the receiver. 


The polar display of either antenna pattern is of the 
form 


= | Vr’ | cos 6 
= | Vr’ | sin 0. 


& 
| 


(S) 


S eptember 


1960 


In order to use a single-channel receiver it is necessary 
either to switch between antennas A and B so as to al- 
low #4 and 6g each to move through 27 radians sepa- 


rately, or to use the same antenna at positions A and B- 


at different times. The two patterns will intersect at a 
common value of 0, say 0; At this intersection the 
amplitudes | Vr'|4 and | Vr’| x are equal so that 


sin (0; — 51) sin (0; — 52) (cos ddA =COS dap) = (0) (6) 
‘Therefore, 
sin (0; = 51) = 0, that is 0, => 61, 01 + TT 


= 


e sin (0; == do) = 0, that is i, = Oe, 6+ a7 


+A 


cos daa = COS dap, that is the patterns are coincident. 


eb Vo cds 2 
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In the case of coincident patterns, there may be more 
than one condition which allows ¢as=¢ap to be ob- 
tained; however, a necessary condition for determining 
the bearings is ¢das#daz. Thus, the intersection azi- 
muths 6; are the directions of the field components with 
a 180° ambiguity. Furthermore, when 6;=46,, 


| Vr’ |ons, = | Ee sin (61 — 85) | 
and when 0;=6:, 


| Vr’ | ocatts = | iy sin (62 = 61) | A 
Since 
| sin «| = | sin (—x)| , 


the display amplitudes at these points are proportional 
to FE; and Ee, the field amplitudes. 


Experimental Comparison of Image Line Radiators 
and Polyrod Antennas* 


g S. P. SCHLESINGER}, MEMBER, IRE, AND A. VIGANTST 


Summary—Experimental comparison is made for some polyrod 
antennas and the corresponding HE, excited dielectric image line 


‘radiators on a finite image surface, dipole and reflector launched. ~ 


; There is a sufficient amount of correspondence in the radiation in- 
_tensities to justify the use of polyrod design data in design of image 
-Jine antennas. The amount of back radiation is smaller for the image 


‘line antennas. 


section of an HEy excited dielectric image line.® 

It is instructive to compare the practical cases of 

-a polyrod antenna launched from the end of a wave- 

guide and an image line antenna ona finite image surface 
launched from a dipole with a reflector. 

- In this investigation the polyrod antennas were 

launched from windows at the end of RG-96/U wave- 


(te polyrod antenna’ has an analog in a finite 


A * Received by the PGAP, April 20, 1960. This work was sup- 
ported by the AF Cambridge Res. Center, Cambridge, Mass., under 
Contract No. AF 19(604)3879. i 
| + Dept. of Elec. Engrg., Columbia University, New York, N. Y. 
1G. E. Mueller, and W. A. Tyrell, “Polyrod antennas,” Bell Sys. 
Tech. J., vol. 26, pp. 837-851; October, 1947. 
2D. G. Kiely, “Dielectric Aerials,” John Wiley and Sons, Inc., 
’ New York, N. Y.; 1953. > . 
8 J, Brown, and J. O. Spector, “The radiating properties of end- 
_ fire aerials,” Proc. IEE, pt. B, vol. 104, pp. 27-34; January, L957 
4L. W. Mickey, and M. G. Chadwick, “Closely spaced high di- 
electric constant polyrod arrays,” 1958 IRE NATIONAL CONVENTION 
- REcoRD, pt. 1, pp. 213-224. ; Le fe 
5S. P. Schlesinger, and D. D. King, “Dielectric image lines, 
IRE Trans. ON MICROWAVE THEORY AND TECHNIQUES, vol. 


-MTT-6, pp. 291-299; July, 1958. 


guide into which the dielectric rods were fitted except 
in the case of a large diameter rod of relative dielectric 
constant 1.67 which was launched from a window at 
the end of a small horn. A cross section of the image 
line radiator is shown in Fig. 1. The reflector is located 
do /4 behind the dipole where \» denotes the free space 
wavelength. The reflector diameter is about 0.9 Xo and 
the image surface size is about 32 X» by 18 Xo with the 
long dimension parallel to the rod and the dipole at 
the center of the image surface. The frequency in all 
cases was 34.86 kmc. The length of the rod L is the 
distance from the window to the end of the rod for the 
polyrod antennas and it was arbitrarily taken to be the 
distance from the center of the dipole to the end of rod 


IMAGE SURFACE 


WAVE GUIDE 


Fig. 1—Cross section of HEy excited dielectric 
image line radiator. 
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for the image line radiators. The radiation patterns were 
taken with a horn polarized to detect the electric field 
vector in the E plane of the exciting waveguide for the 
polyrod antennas and in the plane formed by the rod 
axis and the dipole axis for the image line radiators. 
Three sets of radiation patterns obtained in the man- 
ner described in the previous paragraph are shown in 
Fig. 2. The angle 6 is measured from the rod axis and e€ 
and 2a denote the relative dielectric constant and diam- 
eter of the rods. Since there are diffraction effects in the 
main lobes of the image line radiator patterns, an arbi- 
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trary basis for comparison is arrived at by matching 
first sidelobes of corresponding image line and polyrod 
patterns. Masking and diffraction due to launching 
mechanisms can be seen to begin at about 90° with 
masking more pronounced for the image line radiators. 
There is a sufficient amount of correspondence between 
the radiation patterns to justify the use of polyrod de- 
sign procedures for image line antennas. An attractive 
feature of the dipole and reflector launched image line 
section is the reduced level of back radiation as com- 
pared to the polyrod antenna. 
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Fig. 2—Comparison of radiation intensities for HE excited 
dielectric image line and polyrod antennas. 
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Horn Antennas for HF Long-Range Communication* 


H. BRUECKMANNY, FELLow, IRE, AND B. G. HAGAMANT, sENIOR MEMBER, IRE 


; HE well-known pyramidal horn antennas with a 
| rectangular aperture whose sides are parallel to 
those of the rectangular waveguide feed section 
produce radiation patterns with relatively strong side- 
lobes in the E-plane. This undesirable characteristic is 
explained by the fact that the amplitude distribution 
of the electric field in the aperture parallel to the E- 
plane is essentially uniform. In the H-plane the ampli- 
tude distribution is sinusoidal and the sidelobes are 
much weaker. B. G. Hagaman has proposed a simple 
modification! which reduces greatly the sidelobes of 
horn antennas. It has been investigated under a con- 
tract sponsored by the U. S. Army Signal Research 
and Development Laboratories, Fort Monmouth, N. J., 
which was guided in the direction of technical effort by 
the co-author. Independently, G. F. Koch has also 
‘made a proposal along similar lines for reflector type 
antennas.2 Hagaman proposed to align the diagonals, 
rather than the sides, of the aperture parallel to the 
‘sides of the waveguide feed section, and to shape prop- 
erly the aperture. 
_ The accompanying change of the amplitude distribu- 
‘tion in the aperture was measured using the unsym- 
‘metrical configuration shown in Fig. 1. For the sake of 
brevity, this configuration was named TAHA, short 
for Tapered Aperture Horn Antenna. The lower corner 
‘of the aperture was cut off in order to keep its vertical 


dimension to a minimum. This is a modification of prac- - 


Fig. 1—Dimensions of model TAHA at 2000 mc (160:1 model). 


* Received by the PGAP, May 6, 1960. This paper is the sub- 
stance ofa presentation to the NTG, Prof. Group “Antennas (a 
German professional society of electrical engineers), Munich, Ger- 

many, on April 29, 1960. The complete paper will be published in the 

TG Journal. 

2 ip u S. Army Signal Res. and Dev. Lab., Ft. Monmouth, N. J. 
¢ Developmental Engrg. Corp., Ft. Evans Lab., Leesburg, Va. 

1 B. G. Hagaman, U. S. Patent No. 2,851,686, filed June 28, 1956. 

2G. F. Koch, “Richtantennen mit besonderen Aperturformen, 

Nachrichtentech. Z., vol. 70, pp- 175-186; April, 1957. Also, Fliichen- 
-strahler mit kleinen Nebenmaxima,” Fernmelde tech. Z., vol. 7, pp. 


498-509; October, 1954. 


tical importance, as will be discussed later, but has no 
bearing on the TAHA principle. Disregarding this de- 
tail the aperture is a rhombic with an acute angle of 
about 60 degrees. The function of the peculiar shape of 
the aperture and the particular alignment of its excita- 
tion is explained in Fig. 2, if only qualitatively. It shows 
the lines of the electric field across the aperture as one 
would expect them to flow, from a theoretical consider- 
ation of the conducting boundaries and the chosen ex- 
citation parallel to the long diagonal. Obviously, the 
vertical components of the electric field vectors above 
and below the long diagonal symmetrical to it have 
opposite directions and essentially cancel with respect 
to radiation. The horizontal components, however, 
have the same directions and add with respect to radia- 
tion broadside to the aperture. Moreover, Fig. 2 allows 
realizing qualitatively that the integral of the excita- 
tion over a vertical narrow strip tapers from a maximum 
at the center gradually to zero at the ends of the long 
diagonal. The results of measurements at three differ- 
ent frequencies shown in Fig. 3 confirm this explana- 
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Fig. 3—Measured E-field distribution of TAHA aperture. 
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tion and also prove that the effect is little dependent 
on frequency. 

It might be argued that the explanations and evi- 
dence given for the #- and H-planes do not preclude 
strong new sidelobes arising in directions between the 
two principal planes. This argument is disproved by 
the experimental evidence presented in Fig. 4, which 
shows contour map plots of the radiation in the entire 
forward hemisphere over a wide range of frequency. The 
radiation plotted is the resultant of the field strength 
components measured separately in orthogonal planes 


of polarization. In calculating the resultant it was arbi- 
trarily assumed that these components are in phase. 
This is the most unfavorable phase relation possible 
from the viewpoint of sidelobe suppression. The actual 
resultant may be up to 3 db weaker, which occurs when 
the phase difference is 90°. Note that all sidelobes are 
more than 20 db down from the peak of the main lobe. 


The sidelobes of an optimum-design TAHA made of 
sheet metal are actually in the order of 30 db down over 
a frequency range of more than 1:3. : 

An immediate practical application of TAHA has | 
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Fig. 4—Resultant radiation pattern of experimental symmetrical TAHA 
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Shown are all contour lines for —20 db field strength level and above, 
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been found as a receiving antenna for HF long-range 
communication systems, making use of the reflection of 
electromagnetic waves at the ionosphere. A simple and 
effective solution for the design of the feed section was 
found in the so-called ridge type waveguide depicted in 
Fig. 5. The two ridges taper gradually into the throat 
of the horn. The cone between them is isolated at its 
tip and represents the output terminal. Another diffi- 
cult design problem, if not the most difficult, is the de- 
: ign of the wire grids which replace the sheet metal walls 
of conventional horns for obvious practical reasons. 

- The electrical performance of the completed antenna, 
A photograph of which is shown in Fig. 6, conformed 
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Fig. 5—Ridge waveguide feed system. (Perspective view.) 


Fig. 6—Full-scale TAHA. 
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essentially to the expectations derived from model 
measurements. The major dimensions are given in the 
sketch of Fig. 7. A SWR of 4:1 can be tolerated in re- 
ceiving applications because it corresponds to a mis- 
match loss of only 2 db. Using this criterion, the fre- 
quency range covered is 6.5 to 23 mc. Fig. 8 shows a 
sample of radiation patterns measured by means of a 
helicopter. Fig. 9 summarizes the pattern measure- 
ments as to directivity gain and half-power beamwidth 
in azimuth and elevation vs frequency. 


3H. Brueckmann, “Helicopter measures antenna patterns,” 
Electronics, vol. 28, pp. 134-136; November, 1955. 
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Fig. 8—Comparison of measured patterns of the TAHA and of the 


rhombic at 12 mc. ——— Azimuthal pattern of rhombic (E¢); 
envelope of patterns for y=6°, 9.2°, 14. a Seer: 
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Fig. 9—Characteristics of TAHA vs frequency. 
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Fig. 10 illustrates a certain pattern characteristic of 
TAHA which distinguishes it favorably from most other 
type antennas, including rhombics. The characteristic 
referred to is the change of the elevation angle of the 
main lobe maximum vs frequency. Usually, this angle 
is about inversely proportional to frequency. Such a 
strong dependance on frequency is not desirable from 
the viewpoint of adaptation to normal ionospheric 
propagation. The “angle of fire” of TAHA depends less 
strongly on frequency and comes closer to an optimum 
dependence. The reason is that the center of gravity of 
the excitation in the aperture decreases in height above 
ground with increasing frequency. The main advantage 
of TAHA over a rhombic of equal gain lies, of course, 
in the much lower sidelobes which differ by a factor of 
10 to 15 db.* The disadvantage is the higher cost of 
construction. The experience with the first full-scale 
antenna has uncovered a number of possibilities to re- 
duce the cost; for example, by arranging the supporting 
towers differently. Others are being investigated, such 


4H. Brueckmann,” Suppression of undesired radiation from 
directional HF antennas and associated feed lines,” Proc. IRE, vol. 
46, pp. 1510-1516; August, 1958. 
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Fig. 10—Comparison of TAHA vs rhombic. Measured gain: 1) 
TAHA directivity. 2) Rhombic directivity. 3) Rhombic power. 
Elevation angle: 4) TAHA major lobe. 5) Rhombic major lobe. 
Rhombic referred to has leg length of 400 feet, tilt angle of 71° 
and height of 130 feet. | 


| 
as shortening of the horn. A current model investiga- 
tion of horns for dual-polarization employing the 
TAHA principle shows promise of success. This would | 
permit substituting one horn antenna operated in po-. 
larization diversity reception for two large rhombic an- : 


tennas operated in space diversity reception. 


The Measurement of Self and Mutual Impedances* 
E. ALTSHULER} 


Summary—There are two well-known methods which can be 
used to measure the self and mutual impedances of identical, parallel 
antennas, the symmetric-antisymmetric method, and the open cir- 
cuit-short circuit method.’ In principle, these methods are quite 
simple and straightforward. In practice, however, major difficulties 
arise in both methods. In this paper, it is shown that by using a part 
of each method, the self and mutual impedances can be measured 
quite easily and to a good degree of accuracy. 


ONSIDER the set of linear equations which re- 
late the driving voltages, input currents, and self 


and mutual impedances of identical antennas. 


Vio = (Zs1 += Zri)I10 = Z 01 09 
V20 = ZL 16 SP (Zs2 = Z12)1 9. 


' eres the PGAP, May 19, 1960, 

ordon McKay Lab., Harvard University, Cambridge, M 
TR. W. P. King, “The Theory of Linear tenting ete 
University Press, Cambridge, Mass., pp. 346-347; 1956. 


Zrii and Zz2 are impedances which may be placed in_ 
series with Zs1, or Zs2. Since the antennas are identical, 
Zs1=Zg2. Since the antennas constitute a reciprocal 
network, Z1.= Zo. 

In the symmetric-antisymmetric method, both an- 
tennas are driven with equal voltages, first in phase 
(symmetric), then 180° out of phase (antisymmetric), - 
and the corresponding symmetric and antisymmetric 
input impedances are measured. Zm=Z12=0 


Symmetric case: 
Vio = Zsi1 10 =n Z 21 29 
Vio re Zi21 10 = Z il 2. 


It follows, 


s Vio 
ZS = as = Z51 + Zio. 
10 


960 


Antisymmetric case: 


Vio = Zgilio + Ziol 
% —Vio = Zielio + Zgil 20. 
It follows, 
: Vio 
2 LA a 753 — 212 
4 10 
Therefore, 
ms ZS + Za 
= PR a en me 
s 2 
i” ZS es Ze 
+a OE gs aces eas a 
A y) 
tJ 


i 

This procedure is in general quite difficult to carry 
through experimentally, especially if identical antennas 
are not driven by identical transmission lines. Probably 
‘one of the best methods that can be used to determine 
whether the antennas are driven 180° out of phase is to 
place a small auxiliary probe midway between them 
and then adjust the phase and amplitude of one of the 
driving voltages until a null is detected by this probe. 
This is an indication that the fields midway between 
the antennas cancel and therefore the input currents as 
well as the driving voltages are 180° out of phase. Since 
‘the null is very sharp, the antisymmetrical condition 
can be obtained quite precisely. The phase and ampli- 
tude of one of the driving voltages can be varied quite 
easily by inserting a line stretcher and variable atten- 
‘uator in one of the transmission lines driving the an- 
‘ttenna. A similar procedure cannot be used to obtain 
‘the symmetric condition since a broad maximum exists 
at the auxiliary probe. In order to obtain this condition, 
the line stretcher could be changed by 180° from the 
‘antisymmetric condition, but in general there are small 
‘reflections associated with line stretchers, thereby mak- 
‘ing it difficult to produce an accurate 180° shift with- 
out also slightly changing the amplitude. Naturally, it 
‘is not impossible to drive two antennas symmetrically. 
However, it is quite difficult to accomplish this pre- 
cisely. 

_ In the open circuit-short circuit method, one antenna 
‘is driven while the other remains parasitic; therefore, 
“Vn=0. For the driven antenna, Z11=9. 
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With the input terminals of the parasitic antenna 
open circuited, the input impedance of the driven an- 
tenna is measured. 


Vip = oO; Io) = O. 


10 
Thus, the self impedance of the antenna is simply the 
input impedance of the driven antenna with the para- 
site open circuited. The open circuit condition is ob- 
tained by locating a small current probe at the input of 
the parasite and then adjusting the position of a short 
circuit on the transmission line which feeds the para- 
site so that a null appears at the input. The input im- 
pedance of the driven antenna is then measured with 
the parasite short circuited. 

Vi0 Zi2° 


LG a eas 
Tio 81 


Z12 = 0. 


The short circuit condition is most easily obtained by 
simply short circuiting the input terminals of the para- 
site. Solving for Z12, 


Ze = + VZs1(Zs1 — Zs.c.). 


As can be seen, there exists an ambiguity in the sign of 
the mutual impedance. In order to determine the cor- 
rect sign, an additional measurement is required. The 
principal limitation of this method is that it cannot be 
used to measure small mutual impedances accurately. 
The determination of the mutual impedance depends 
on the difference of the input impedances of the driven 
antenna when the parasite is open and short circuited. 
This difference, which is usually quite small, is then 
multiplied by the self impedance, which is relatively 
large. In this process, the accuracy of the mutual im- 
pedance is reduced. 

By using a part of each of the previous methods, the 
self and mutual impedances can be measured quite 
easily and to a good degree of accuracy. As has been 
stated earlier, the self impedance of an antenna is 
simply the input impedance of the driven antenna when 
the parasite is open circuited. This measurement is 
quite exact. The antisymmetric impedance which can 
also be measured quite precisely is by definition: 


whe = Zs1 ee Z12. 


Therefore, the mutual impedance is simply the differ- 
ence between Z* and Zs1. 
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An Island as a Natural Very-Low-Frequency 
Transmitting Antenna* 


M. G. MORGANT, SENIOR MEMBER, IRE 


N thinking about how to generate whistlers arti- 
] ficially, it has occurred to me that an island of 

suitable size and shape, extending through the con- 
ducting sea, may constitute a naturally resonant, VLF 
slot antenna of high quality. Following this line of rea- 
soning, I thought first of the annular Pacific atolls, but 
knowing of the fresh-water lenses in them, rejected 
them as being too pervious to water to be satisfactory 
insulators. Also, of course, they are not found in suit- 
able latitudes for generating whistlers. The Pacific 
atolls are built upon submerged volcanic cones and this 
led me to think of Deception Island (Fig. 1, opposite) in 
the SubAntarctic, a remarkable, similarly shaped, vol- 
canic island in which the volcanic rock extends above 
the surface; and which is located in the South Shetland 
Islands where the rate of occurrence of natural whistlers 
has been found to be very great. Deception Island is a 
thin annulus surrounding a lagoon which has a single 
narrow entrance through the ring. The length of the 
island is about 30 km which is a half-wavelength at 5 kc, 
a frequency at which natural whistlers are very strong. 
To a rough approximation, the radiation pattern would 
be that of a vertical monopole, that is, proportional to 
the sine of the zenith angle. 

Taking the conductivity of sea water to be 4 mhos/ 
meter and the relative dielectric constant to be 80, one 
calculates the ratio of conduction current density to 
displacement current density (¢/we) at 5 kc to be 
1.8 X10°, so that the sea is certainly a good conductor. 
The “depth of penetration” (V2/wuc) is only 3.5 
meters so that depths in the lagoon, which run from 110 
meters near the shore to 165 in the middle, appear to be 
adequate. The entrance, called Neptune’s Bellows, is 
635 meters wide. For half of this width, the water is only 
3 to 4 meters deep, but there is a channel about 150 
meters wide which is 25 meters deep. The salinity in 
the lagoon does not differ significantly from that of the 
outside sea. 

The island appears to be constituted mainly of mas- 
sive tuff, a geological material and form which could be 
very favorable to its being a good insulator. The quality 
of the insulator most likely depends entirely upon the 
water content due to porosity, stratification, and frac- 
turing. It would be a relatively simple matter to make 


* Received by the PGAP, June 8, 1960. 
} Thayer School of Engrg., Dartmouth College, Hanover, N. H 


an im situ measurement of the driving point impedance 
of the proposed island antenna over a range of frequency 
and this would give a decisive answer to the critical 
question of the quality of the island as an insulator. 


The highest impedance across a half-wave resonant slot. 
is found at its center which, for the island antenna, 


would be between the outside sea and the lagoon at a 


( 


4 
| 


point opposite to the entrance to the lagoon, a point | 


aptly named Telefon Bay. 


At Telefon Bay, it is 2 km across the island to the | 


outside so that from the middle the connecting lines to 
the sea terminals would each be 1 km long. The induc- 
tive reactance of these leads will not be negligible. Inas- 


much as the current will not return in the ground di- | 


rectly under the wires if the quality of the island as an 
insulator is sufficiently good (the wires should be insu- 
lated from the moist surface), we use the expression 
(u/m) In (b/a) for the inductance per unit length of a 
parallel wire transmission line of conductor separation 
b and radius a. We take 6=2/ where / is the length 
from the center of the island to a point along either wire 
and obtain the total inductance as 


21 
i(1m — 1) 3 
r 


For /=1 km, r=1 mm (approximately No. 12 Amer- 
ican wire gauge), and M@=47 X10-’, the inductance is 
5.4 mh which leads to a reactance of 170 ohms at okt 

The inductance of the connecting wires can be re- 


21 
yr eee 


1 
Li 
=f. In ] 
Tv 0 oe T 


duced by running multiple wires in parallel. For uni- © 


formally spaced wires, the inductance is given by 


E + (n— 1)M 


nN 


=U xX 10-7] 


where L is the inductance without the added conduc- 
tors, M is the mutual inductance of any pair of adja- 
cent conductors, ” is the number of conductors on each 
side, and k is given by 

n 2 4 6 8 10 


k 0 0.621 We 


M is given by 


2 
21 X 10-7| n= os 1] 
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Fig. 1—Deception Island. There are fumaroles on the beach at Whalers Bay and De Mayo Bay. In 1930 an earthquake occurred 
at Whalers Bay, lowering the bottom nearly 5 metres. 
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where d is the separation of adjacent wires. The follow- 1 ohm for y=1/16m or a sphere whose area is but fifty 


ing calculations show that the inductive reactance can 
be lowered to 35 ohms by using 10 wires, 10 meters 
apart, which would certainly be practical. 


d=5 meters 


2 4 6 8 10 20 
EL (mh) 3.20) 2202 91.62) 212385) 1-20) 10790) 
X (ohmsat5kc) 101 64 51 43 38 28 


d=10 meters 


n 2 oa 6 8 10 
£ (mh) 3.13 1593 1.50 1.25 ilgili 
X (ohmsat5ke) 98 60 47 39 35 


The resistance of the wires will be negligible and the 
connections to the sea present no problems either. To 
investigate the latter point, we make use of the fact 
that if the capacitance between two electrodes is C, the 
resistance between them when immersed in a medium 
of conductivity o and dielectric constant e, is e/a C. The 
capacitance of an isolated sphere of radius r is 47er giv- 
ing the resistance to a surrounding, boundless, conduct- 
ing medium as 1/4mor. For o=4, the resistance is 
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square centimeters which is less than the surface area 
of one-meter of wire of 1 mm radius. Thus, though the 
configuration has been idealized, it is obvious that it 
will suffice to run each wire out a few meters from shore, 
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Notes 


1) For a general description of Deception Island, see “Sailing | 
Directions for Antarctica,” U. S. Navy Hydrographic Office, pp. 101— 
104, 1943, and supplement of 1959; or “The Antarctic Pilot,” U. K. | 
Admiralty Hydrographic Dept., pp. 164-166, 1948. : 

For a detailed scientific description, see “Contribucion al Cono- 
cimiento del Sector Anartico Argentino,’ Instituto Antartico Ar-— 
gentino, publ. no. 1, 1955; and “Contribucion a la Geologia de la 
aes Occidental,” Instituto Antartico Argentino, publ. no. 2, 
1 : 

2) MKS rationalized units are used. The formula for the reduc- 
tion of inductance by using multiple wires, is from Circular 74, U. S. 
National Bureau of Standards, p. 249; 1924 (reprinted 1937). 
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Canada. His work has been mostly in the 
field of antenna, propagation, and system 
engineering studies. Since 1958, he has been 
responsible for the theoretical aspects of the 
microwave subsystem in a wide band micro- 
wave communication equipment. 
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Robert C. Hansen (S’47—A’49-M’55- 
SM’S56) was born on August 3, 1926, in St. 
Louis, Mo. He attended the Missouri School 

- of Mines and Metal- 
lurgy, Rolla, from 
which he received the 
B.S.E.E. degree in 
1949. He received the 
M.S. and Ph.D. de- 
grees in 1950 and 
1955, respectively, 
from the University 
of Illinois, Urbana. 

While at the Uni- 
versity of Illinois, he 
worked at the Anten- 
na Laboratory on 
ferrite loops and streamlined airborne an- 
tennas. From 1955 to 1959, he was at the 
Microwave Laboratory of Hughes Aircraft 
Company, Culver City, Calif., where, as 
senior staff engineer, he worked on surface- 
wave antennas, slot arrays, and related 
fields. In 1959, he became senior staff 
engineer in the Telecommunications Lab- 
oratory of Space Technology Laboratories, 
Los Angeles Calif., where he is at present. 

Dr. Hansen is a member of the American 
Physical Society, Commission VI of URSI, 
Tau Beta Pi, Sigma Xi, Eta Kappa Nu, and 
Phi Kappa Phi. 
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Julius Kane was born on January 23, 
1935, in Vorozenh, Russia. He received the 
B.A. degree from Brooklyn College, Brook- 
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lyn, N. Y. in 1954, majoring in physics. 

Following graduation he was employed 
as a radiation physicist in cancer re- 
search. Subsequently, 
he spent a summer 
at Sylvania Electric 
Products, Bayside, 
at NG Verasjanres 
search physicist, work- 
ing on problems 
of electroluminescent 
displays. He was a 
lecturer in physics at 
the College of the 
City of New York, 
and has done post- 
graduate work at New 
York University, New York, N. Y., where 
he received the Ph.D. degree in applied 
mathematics. Currently, he is a staff con- 
sultant with Dorne and Margolin, Inc., 
Westbury, L. I., N. Y., engaged in research 
on antennas and propagation. In July, 1960, 
he became an associate professor of Elec- 
trical Engineering at the University of 
Rhode Island, Kingston. 

Dr. Kane is a member of the American 
Physical Society and the American Mathe- 
matical Society. 
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Aharon A. Ksienski (S’53-M’58) was 
born on June 24, 1924, in Warsaw, Poland. 
He received the M.S. degree in 1952 and the 
Ph.D. degree in 1958, 
both in electrical en- 
gineering from the 
University of South- 
ern California, Los 
Angeles. 

From 1948 to 
1951, he was in charge 
of an Aircraft Elec- 
trical School in the 
Israeli Air Force. 
From 1951 to 1952, 
he attended the Uni- 
versity of Southern 
California as a graduate student and a part 
time member of the teaching staff. Between 
1953 and 1955, he was with Wiancko Engi- 
gineering Company, Pasadena, Calif., en- 
gaged in developing instrumentation sys- 
tems for the measurement of various flight 
parameters and in the design of electro- 
mechanical transducers. In 1955, he returned 
to the University of Southern California as 
4 lecturer in mathematics and electrical en- 
gineering, and continued graduate research 
in the field of antennas. After graduating in 
1958, he joined the Antenna Research De- 
partment in the Microwave Laboratory of 
Hughes Aircraft Company, Culver City, 
Calif., where he is presently a senior staff 
engineer engaged in research in the field of 
antennas and antenna systems. 

Dr. Ksienski is a member of Phi Kappa 
Phi, Tau Beta Pi, Eta Kappa Nu, and Sigma 
Xi. 


A. A. KSIENSKI 
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Curt A. Levis (S’48—A’52—M’57-SM’59) 
was born in Berlin, Germany, on April 16, 
1926. During World War II he served in the 
U. S. Navy as an 
electronic technician 
and subsequently 
studied electrical en- 
gineering at Case 
Institute of Technol- 
ogy, Cleveland, Ohio, 
where he received the 
B.S. degree in 1949, 
During this time he 
was also employed 
as studio engineer 

C. A. Levis at Radio Station 

WSRS. 

He was a Gerard Swope Fellow at Hary- 
ard University and received the M.A. degree 
in engineering science and applied physics 
there in 1950, and the Ph.D. in electrical 
engineering at The Ohio State University, 
Columbus, in 1956. 

Since 1950, Dr. Levis has been employed 
at the Antenna Laboratory of The Ohio 
State University. His present position is 
associate supervisor of the Antenna Lab- 
oratory and associate professor of electrical 
engineering. 
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Donald. R. McCoy (A’54—M’58) was 
born in Burkeburnette, Tex., on August 30, 
1920. He attended The Ohio State Univer- 
sity, Columbus, from 
which he received the 
B.S.E.E, degree in 
1951 and the M.S. 
degree in 1956. 

He was an in- 
structor in the Signal 
Corps Radar schools 
from 1942 to 1945, 
and worked with the 
Signal Corps Engi- 
neering Laboratory 
Development De- 
tachment from 1945— 
1946. He was with the Antenna Laboratory 
of The Ohio State University Research 
Foundation from 1949-1956, He joined the 


D. R. McCoy 
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General Electric Company, Syracuse, N. Y., 
in 1956, and is presently engaged in Sys- 
tems Engineering. 

Mr. McCoy is a member of Sigma Xi. 
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Ralph J. Rainey (S’59) was born in AI- 
buquerque, N. M. on April 18, 1931. He 
received the B.S.I.A. degree in 1953 and 
the B.S.E.E. degree 
in 1959, both from 
the University of 
New Mexico, Albu- 
querque. He served 
as a rated officer in 
the Air Force for 
three years from 1953 
to 1956. 

Since September, 
1958, he has been em- 
ployed in the Elec- 
trical Engineering 
Division of the 
Engineering Experiment Station, Univer- 
sity of New Mexico, where he is engaged in 
tropospheric propagation studies. 

Mr. Rainey is a member of Sigma Tau 
and an associate member of Sigma Xi. 


R. J. Rainey 
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Sheldon S, Sandler, for a photograph and 
biography, please see p. 453 of the July, 
1960 issue of these TRANSACTIONS, 
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Howard E. Shanks was born in Los 
Angeles, Calif.,on January 28, 1933. He at- 
tended Pasadena City College, Pasadena, 
Calif., from 1950 to 1952, and then trans- 
ferred to the California Institute of Tech- 
nology, Pasadena, where he received the 
BiS.E.E. degree in 1954. In 1956, he re- 
ceived the M.S. degree in physics from the 
University of California at Los Angeles 
under the Hughes Fellowship program. 

He was with the Antenna Research De- 
partment of the Hughes Aircraft Com- 


Septembe: 


pany, Culver City, Calif., for six years an | 
contributed to research in the areas of UHF- 


VHF antennas, scattering cross sections, 


circular-aperture an- 


tennas, flush-mounted 
radar antennas, ar- 
bitrarily polarized ra- 
diators, advanced an- 
tenna techniques, and 
pattern synthesismeth- 
ods. He joined Amer- 


Inglewood, Calif., in 
March, 1960. He is 
the author of numer- 
ous technical reports 
and journal articles, 


H. E. SHANKS 


and has seven patent applications in the 


field of antennas and microwave devices. 
He is a lecturer in electrical engineering at 


the University of Southern California, Los_ 


Angeles. 
Mr. Shanks is a member of Tau Beta Pi 
and RESA. 
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Carlton H. Walter (S’48—A’51-M’56-— 
SM’59) was born at Willard, Ohio, on July 
22, 1924. He received the B.E.E. degree at 
The Ohio State Uni- 
versity, Columbus, 
in June, 1948, at 
which time he joined 
the staff of The Ohio 
State University An- 
tenna_ Laboratory. 
He did graduate work 
while employed at 
the Antenna Labora- 
tory and received the 
M.S. degree in phys- 
ics in March, 1951 
and the Ph.D. in 
electrical engineering in June, 1957. 

He is now an associate supervisor at The 
Ohio State University Antenna Laboratory 
and an assistant professor in the Electrical 
Engineering Department. 

Dr. Walter is a member of Sigma Xi, 


C. H. WALTER 


Tau Beta Pi, Eta Kappa Nu, and Sigma Pi | 


Sigma, and a Registered Professional Engi- 
neer in Ohio. 
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FOR INFORMATION CONCERNING ADVERTISING RATES 


Contact 
MR. DELMER C. PORTS 
Jansky and Bailey, Inc. 


1339 Wisconsin Avenue N.W. 
Washington 7, D.C. 


Telephone: Federal 3-4800 


FILED IN STACKS . 
INSTITUTIONAL LISTINGS 


The IRE Professional Group on Antennas and Propagation is grateful for 
the assistance given by the firms listed below, and invites application for 
Institutional Listing from other firms interested in the field of Antennas and 
Propagation. 


AERO GEO ASTRO CORP., 1200 Duke St., Alexandria, Va. ae q 
Research and Development; Antennas; Transponders; Command Receivers; Augmenters; Telemetry-Radar 


ANDREW CORPORATION, P.O. Box 807, Chicago 42, IIl. P 
Antennas, Antenna Systems, Transmission Lines, Development and Production. ’ 


ANTLAB, INC., 6330 Proprietors Rd., Worthington, Ohio 
Antenna Pattern Range Systems— Recorders & Mounts, & Telemetry Servo Pedestals. 


BLAINE ELECTRONETICS, INC., 14757 Keswick St., Van Nuys, Calif. 
Antennas, Paraboloids, Scale Models, Antenna Radiation Pattern Measurement Towers. 


. 
GABRIEL ELECTRONICS, Division of The Gabriel Company, Main & Pleasant Sts., Millis, Mass. — 
Research, Engineering and Manufacture of Antenna Equipment for Government and Industry. 


HUGHES AIRCRAFT COMPANY, Florence and Teale Sts., Culver City, Calif. 
Res., Dev., Mfg.: Radar Systems & Components; Antennas, Tubes, Radomes, Solid-State Devices. 


I-T-E CIRCUIT BREAKER CO., Special Products Div., 601 E. Erie Ave., Philadelphia 34, Pa. 
Design, Development and Manufacture of Antennas, and Related Equipment. 


ee 


JANSKY & BAILEY, INC., An Affiliate of Atlantic Research Corp., 
1339 Wisconsin Ave., N.W., Washington, D. C. , ; 
Complete Engineering Services for Antennas and Propagation Programs. 


MARK PRODUCTS CO., 5439 W. Fargo Ave., Skokie, II. } 
Antennas for Two-Way Communications, Grid Parabolas, Research & Development. 


TECHNICAL APPLIANCE CORP., 1 Taco St., Sherburne, N. Y. 
Des., Dev., & Mfg.: Antennas & Antenna Systems for Communications, Telemetry, & Tracking. 


WEINSCHEL ENGINEERING COMPANY, INC., Kensington, Md. 
Antenna Pattern Receivers; Bolometer Amplifiers; Modulated Microwave Sources; 
Insertion Loss Measuring Systems. 


WHEELER LABORATORIES, INC., Great Neck, N. Y.; Antenna Lab., Smithtown, N. Y. 
Consulting Services, Research and Development, Microwave Antennas and Waveguide Components. 


ate charge for Institutional Listing is $25 for one issue or $100 for six consecu- 
tive eaues (one year). Application may be made to the Technical Secretary, 
The Institute of Radio Engineers, 1 East 79th Street, New York 21, N. Y. 


